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Preface 



By intention, my project has two parts. The first one covers the classical 
electron theory. It is essentially self-contained and will be presented in the 
following chapters. 75 years after the discovery of quantum mechanics, to 
discuss only the classical version of the theory looks somewhat obsolete, in 
particular since many phenomena, like the stability of atoms, the existence of 
spectral lines and their life time, the binding of atoms, and many others, are 
described only by the quantized theory. Thus it is a necessity to discuss the 
quantized version of the classical models studied here. This is not quantum 
electrodynamics. It is the quantum theory of electrons, stable nuclei, and 
photons with no pair production allowed. Well said, but the quantum part 
turns out to be a difficult task. There is a lot of material with the mathe- 
matical physics side in flux and very active at present. Thus it remains to 
be seen whether the quantum part will be ever finished. In the meantime I 
invite the reader to comments, criticisms, and improvements on the classical 
part. 

In thank my collaborators, Sasha Komech and Markus Kunze, for their 
constant help and insistence. I am very grateful to Joel Lebowitz. He ini- 
tiated my interest in tracer particle problems and I always wanted to apply 
these ideas to electrons coupled to the Maxwell field. I am indebted to F. 
Rohrlich for discussions and important hints on the literature. I acknowl- 
edge instructive discussions with A. Arai, V. Bach, D. Bambusi, G. Bauer, 
J. Bellissard, F. Bonetto, D. Diirr, J. -P. Eckmann, L. Erdos, J. Frohlich, L. 
Galgani, G. Gallavotti, S. Goldstein, M. Hirokawa, F. Hiroshima, V. Jaksic, 
M. Kiessling, E. Lieb, M. Loss, D. Noja, C.-A. Pillet, M. Rauscher, L. Rey- 
Bellet, H.T. Yau. 

I am most thankful to Mrs. L. Vasquez de Rosswag for her cheerful and 
tireless effort in typing. 

Herbert Spohn Miinchen, August 1, 1999 
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Introduction 



"Classical Electron Theory" is an attempt of building a dynamical theory 
of electromagnetic fields coupled to well concentrated lumps of charges, like 
electrons and nuclei. As the name indicates, it is a theory of classical fields in 
interaction with classical particles. Two limiting cases have been confirmed 
experimentally and are the subject of any course on electrodynamics: either 
the currents, i.e. the motion of the charges, or the electromagnetic fields 
are given with the task to predict the behavior of the remaining piece. In 
stark contrast, there are only few examples where one is still in the classical 
domain and the full power of a coupled theory is required. In fact, I only 
know one example, namely the motion of a single electron in a Penning trap 
which will be discussed in Section 0. Thus the classical electron theory is a 
mostly theoretical enterprise, but with this taste it has intrigued physicists 
for almost a century. 

The motives have varied with time. The founding fathers, like Abraham 
and Lorentz, tried to develop a dynamic theory for the then newly discovered 
electron. In particular, depending on the model, they predicted its energy- 
momentum relation, cf. Section [| This enterprise came to a stand still with 
the advent of the theory of special relativity. Based on totally disjoint argu- 
ments it required for any massive particle the relativistically covariant link 
between energy and momentum. The classical electron theory flourished for 
a second time in the early days of quantum electrodynamics. The hope was 
that a refined understanding of the classical theory would give a hint on how 
to properly quantize and how to correctly handle the infinities. The notion of 
mass renormalization was taken from the classical theory. But as the proper 
quantum theory surfaced, the classical considerations had little value. In 
fact, there is no classical analogue of the renormalization in quantum field 
theory. There the interaction with the quantized field reduces the quantum 
fluctuations of the electron. The, yet unproven, construction is to counter- 
balance by increasing the bare electron fluctuations through letting its bare 
mass tend to zero in such a way that, upon removing the ultraviolet cut-off, 
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finite fluctuations remain which are adjusted to yield the experimental mass 
of the electron. 

At large, the classical electron theory has a poor reputation. One has to 
fight with infinities, the bare mass of the electron is supposedly negative and 
tends to — oo in the limit of a point charge, the theory has instabilities as 
demonstrated by the runaway solutions, physical solutions suffer from preac- 
celeration in contradiction to a causal description, and more. In contrast, 
our approach will be rather conservative with no need for such spectacular 
revisions of the conventional concepts. It is based on two corner stones 

• a well-defined dynamical theory of extended charges in interaction with 
the electromagnetic field, 

• a study of the effective dynamics of the charges under the condition 
that the external potentials vary slowly on the scale given by the size 
of the charge distribution. This is the adiabatic limit. 

Of course, the familiar objects reappear in disguise with one important dif- 
ference. In the adiabatic limit the ratio 'mass induced through the self- 
interaction with the electromagnetic field' to 'bare mass' remains finite and 
does not diverge as it is the case in the traditional point charge limit. 

Our approach reflects the great progress which has taken place in the 
theory of dynamical systems. After all, a charge coupled to its radiation 
field is just one particular case, however with some rather special features. 
Perhaps the most unusual one is the appearance of a central manifold in the 
effective dynamics if friction through radiation is included. 

A few words on the style and scope of the book are in order. First 
of all we systematically develop the theory, no review is intended. For a 
subject with a long history, such an attitude looks questionable. After all, 
what did the many physicists in that field contribute? To make up, we 
included one historical chapter, which as very often in physics is the history 
as viewed from our present understanding. Since there are already excellent 
historical studies, we hope to be excused. Further we added at the end of each 
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chapter Notes and References which are intended as a guide to all the material 
which has been left out. As prerequisites for reading, only a basic knowledge 
of Maxwell's theory of electromagnetism and of Hamiltonian mechanics is 
required. 

The reader will notice that we state several theorems and give their proof. 
In some cases the proof is complete. In other cases we only present the 
essential idea. The more technical steps can be found in the original literature 
and are therefore not reproduced. I am convinced that, in particularly in a 
subject where there is so little control through experience, one needs fixed 
points in terms of mathematical theorems. Over the years a lot of common 
knowledge has been accumulated and we made an effort to disentangle hard 
facts from truth by tradition. 

The introduction might give the impression as if all problems are resolved 
and the classical electron theory is in good shape. This would be too simplis- 
tic an attitude. What I hope is to convince the reader to view our dynamical 
problem from a particular perspective. Once this point has been reached, 
there are then many loose ends. To mention only two: our comparison be- 
tween the true and approximate particle dynamics could be sharper and a 
fully relativistic theory of an extended charge distribution exists only on a 
formal level. The greatest reward would be if my notes encourage further 
research. 
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1 A Charge Coupled to its Electromagnetic 
Field 



We plan to study the dynamics of a well localized charge, like an electron 
or a proton, when coupled to its own electromagnetic field. The case of sev- 
eral particles is reserved for Chapter || In a first attempt, one models the 
particle as a point charge with some definite mass. If its world line is pre- 
scribed, then the fields are determined through the inhomogeneous Maxwell 
equations. On the other hand, if the electromagnetic fields are given, then 
the motion of the point charge is governed by the Lorentz force equation. 
While it then seems obvious how to marry the two equations, so to have 
a coupled dynamics for the charge and its electromagnetic field, ambigui- 
ties and inconsistencies arise due to the infinite electrostatic energy of the 
Coulomb field of the point charge. Thus one is forced to introduce a slightly 
smeared charge distribution, i.e. an extended charge model. Mathematically 
this means that the interaction between particle and field is cut-off or reg- 
ularized at short distances, which seems to leave a lot of freedom. There 
are strong constraints however. In particular one has to satisfy local charge 
conservation, the theory should be of Lagrangian form, and it has to repro- 
duce the limiting cases mentioned above. In addition, as to be expected for 
any decent physical model, the theory should be well-defined and empiri- 
cally accurate in its domain of validity. In fact, only two models have been 
worked out in some detail: (i) the semi-relativistic Abraham model of a rigid 
charge distribution, and (ii) the Lorentz model of a relativistically covariant 
extended charge distribution. The aim of this chapter is to introduce both 
models at some length. On the way we recall some properties of the Maxwell 
equations for later use. 
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1.1 The Maxwell equations 

We prescribe a charge density p(x, t) and an associated current, j(x, t), which 
are linked through the conservation law 

d tP (x,t) + V-j(x,t) = 0. (1.1) 

Of course, x G M 3 , the physical space, and i e 1, the time. The Maxwell 
equations for the electric field E and the magnetic field B consist of the two 
evolution equations 

c^dtB^x.t) = -VxE(x,t), 

c" l d t E(x,t) = V xB(x,t)-c~ 1 j(x,t) (1.2) 
and the two constraints 

V • E(x, t) = p(x, t) , V ■ B(x,t) = 0. (1.3) 

We use the Heaviside-Lorentz units. The vacuum susceptibilities are 
Eq = 1 = fiQ, which fixes the unit of charge, c is the speed of light. Mostly we 
will set c = 1 for convenience, thereby linking the units of space and time. If 
needed, one can easily reintroduce these natural constants in the convential 
way. At some parts we will do this without notice, so to have the dimensions 
right and to better keep track of the order of magnitudes. 

We solve the Maxwell equations as a Cauchy problem, i.e. by prescribing 
the fields at time t — 0. If the constraints ( |1.3| ) are satisfied at t — 0, then by 
the continuity equation ( |1 . 1| ) they are satisfied at all times. Thus the initial 
data are 

E(x,0), B{x,0) (1.4) 

together with the constraints 

V • E(x, 0) = p(x, 0) , V-B(aj,0) = 0. (1.5) 
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The choice t = is merely a convention. In some cases it is preferable to 
prescribe the fields either in the distant past or the remote future. We will 
only consider physical situations where the fields decay at infinity and thus 
have the finite energy 

£ = \J d3x { E ( x , + B ( x ' t) 2 ) < 00 • (1-6) 

In a thermal state at non-zero temperature, one would be forced to consider 
infinite energy solutions. But this is outside the present scope. 

The Maxwell equations ( |1.2|) , ( |1.3|) are inhomogeneous wave equations 
and thus easy to solve. This will be done in Fourier space first, where Fourier 
transform is denoted by ^ and defined through 

f(k) = (2vr)- n/2 J d n xe~ lh - x f{x). (1.7) 

Then ( |1.2|) becomes 

d t B(k,t) = -ikxE(k,t), 

d t E(k,t) = ik x B(k,t) -j(k,t) (1.8) 

with the constraints 

ik ■ E(k, t) = p(k, t), ik- B(k, t) = (1.9) 

and the conservation law 

d t p(k,t)+ik-j{k,t) = 0. (1.10) 

To solve the inhomogeneous equations (|1.8|), as usual, we rely on the solution 
of the homogeneous equations, 

%^\) = U C os\k\t + (l-cos\k\t)^\k) { k\)(% k ^ 
B [0) (k,t) J L 1 1 V lJ { B(k,0) 

+(-L S m\k\t)ikx ( ^i fc ' 0) I . I II 

1*1 V ~E(k,0) 1 
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We insert ( |L11|) in the time-integrated version of ( |1.8| ). We impose the 
constraints and, with one partial integration and using the conservation law, 
arrive at 



E(k,t) = (cos \k\t)E(k, 0) + (|fc| _1 sin \k\t)ik x B(k,0) 

t 

+ J ds{ - (|fc| -1 sin|Jfe|(t- s))ikp(k, s) - (cos|fc|(i- s))j(k,s)) 
o 

= E (0) (k,t) + E Iet (k,t) , (1.12) 

B(k,t) = (cos|fe|t)B(fe,0) - (|fc| _1 sin|fc|t)iJfe x ^(fe,0) 

t 

+ J ds^k]' 1 sin \k\(t — s))ik x j(k, s) 
o 

= B {0) {k,t) + B vct {k,t). (1.13) 



The first terms are the initial fields propagated up to time t, while the second 
terms are the retarded fields. If we would solve the Maxwell equations into 
the past, then the retarded fields are to be replaced by the advanced fields. 

Let us introduce the fundamental propagator, G t (x), of the wave equation 
which satisfies 



d 2 t G - AG = S(x)S(t) 



;i.l4) 



and is defined as the Fourier transform of (2ir) 3 ^ 2 \k\ 1 sin \ k\t. This means 
Gt(aj) = (27T)- 1 5(\x\ 2 - t 2 ) and for t > 



G t {x) = —5{\x\-t). 



;i.l5) 



Then in physical space the solution ( 1.12 ), (PL. 13 ) of the inhomogeneous 
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Maxwell equations reads as 

E(t) = dtG t * E(0) +V x G t * B(0) - fd3(VGi-.*p(8)+ dtGt-.* j(s)) 

o 

= E {0) (t) + E Tet (t) , (1.16) 

t 

B(t) = d t G t * B(0) - VxG,* £(0) + J dsV x G t _ s * j(s) 

o 

= S (0) (t)+B ret (t). (1.17) 

Here * denotes convolution, i.e. /x * /^(a;) = / d n yfi(x — y)f<i{y)- 

The expressions ( |1.16|) , (|1.17|) remain meaningful when p,j are generated 
by the motion of a single point charge. Let us denote by q(t) the position 
and v(t) = q(t) the velocity of the particle carrying the charge e. Then 

p(x, t) = e6(x - q(t)) , j(x, t) = e5(x - q(t))v(t) . (1.18) 

It is assumed that the particle is relativistic and therefore \v(t)\ < 1. Insert- 
ing in ( 1.16Q , Q1.17]) one arrives at the Lienard-Wiechert fields. Since their 



derivation is handled in any textbook, we do not repeat the computation here 
and discuss only the result. We imagine that the world line, t i— > q(t), of the 
particle is given for all times. We prescribe the initial data for the fields at 
time t = to and take the limit to ~~ > ~°° m ( 1-16|) , ( |1 . 1 7|) . Then, at a fixed 



space-time point (t,x), the contribution from the initial fields vanishes and 
the retarded fields become the Lienard-Wiechert fields. To define them we 
introduce the retarded time £ ret , depending on t, x, as the unique solution of 

t Tet = t-\x-q(t iet )\. (1.19) 

t rct is the unique time when the world line crosses the backwards light cone 
with apex at (t, x). We also introduce the unit vector 

1-E q{ l ret)\ 
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Then the electric field generated by a moving point charge is given by 

(1 — v 2 )(n — v) nx [(n — v) x v] 



E(x,t) = (e/4n) 
and the corresponding magnetic field is 



^1 — v ■ n) 3 \x — q\ . 



\t = tr. 



1.21) 



B(x,t) = fix E(x,t) . 



(1.22) 



( |1.21| ) and ( |1.22| ) are less explicit than the notation suggests, since £ ret de- 
pends through (|1.19|) on the reference point (t, x) and the particle trajectory. 



The first contribution in ( 1.21 ) is proportional to \x — q\~ 2 and independent 
of the acceleration. This is the near field which in a certain sense remains 
attached to the particle. The second contribution is proportional to \x — q\ _1 
and to the acceleration. This is the far field, which carries the information 
on the radiation field escaping to infinity. If q(t) is smooth in t, then the 
Lienard-Wiechert fields are smooth functions except at x — q(t), where 
they diverge as \x — q(t)\~ 2 . The corresponding potentials have a Coulomb 
singularity at the world line of the particle. 



1.2 The Lorentz force equation 

We take now the point of view that the electromagnetic fields E, B are given. 
The motion of a charged particle, charge e, position q(t), velocity v(t), is then 
governed by the Lorentz force equation 

j t (m ol v(t)) = e(E(q(t),t) + c~ l v(t) x B(q(t),t)) , (1.23) 

which as ordinary differential equation has to be supplemented by the initial 
conditions q(0),v(0). Here "f(v) = l/y/T— (v/c) 2 . The particle is rela- 
tivistic with rest mass mo as measured experimentally. Once the particle is 
coupled to the Maxwell field, mo will attain a new meaning. 

The (J5, B) fields in (|1.23| ) are not completely arbitrary. They have to 
be solutions of the Maxwell equations with sources (p,j). So to speak, we 
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have separated all charges into a single charged particle whose motion is to 
be determined through (|1.23|) and a rest whose motion is taken to be known. 

( |1.23| ) is of Hamiltonian form. To see this we introduce the vector poten- 
tials 4>, A such that 

E(x, t) = -V<j>(x, t) - c~ l d t A(x, t) , B(x, t) = V x A(x, t) . (1.24) 

Then the Lagrangian associated to (|1.23| ) is 

L(q, q) = -m c 2 (l - c^q 2 ) 1 ! 2 - e(0(q, t) - c~ x q ■ A(q, t)) . (1.25) 



Introducing the canonical momentum 

g 

p = m l(q)q + -A(q,t) 

c 

the Hamiltonian is given by 

H(q,p,t) = ((cp-eA(q,t)) 2 + m2 c 4 ) 1 /2 + e 0(g )t ). 

In particular, if the fields are time-independent, then the energy 

S(q,v) = m 7(v) + e(f>(q) 



1.26) 



(1.27) 



(1.28) 



is conserved along the solution trajectories of (|1.23|) . 

It should be remarked that in general the solutions to the Lorentz force 
equation will have a complicated structure even for time-independent fields. 
This has been amply demonstrated for particular cases. Depending on the 
external fields the motion ranges from regular to fully chaotic with a mixed 
phase space as a rule. 



1.3 The coupled Maxwell— Lorentz equations 

While, up to a minute fraction, sufficient for all of electrodynamics, from a 
more fundamental point of view it is unsatisfactory that in the Maxwell and 
Lorentz force equations so to speak "one half" must be prescribed and one 
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would hope to have a coupled system of equations for the time evolution of 
the charged particles together with their electromagnetic field. If we restrict 
ourselves to a single particle, it is obvious how to proceed. From ( |1.2|) , (|1.3|) 
we have 

d t B(x,t) = -VxE(x,t), 

d t E(x,t) = V x B(x,t) - eS(x - q(t))v(t) (1.29) 
with the constraints 

V-E(x,t) = e6{x-q{t)), V-B(x,t) = (1.30) 
and from (11.231) we have 



-{m Ql v{t)) = e(E ex (q(t)) + E(q(t),t) + v(t) x (B ex (q(t)) + B(q(t), t))) . 

(1.31) 

We added in explicitely external electromagnetic fields E ex , B ex , which will 
play a prominent role later on. They are derived from potentials as 

E cx = -V0 CX , B cs = Vx A ex . (1.32) 

We assume the potentials to be time-independent for simplicity, although 
considerable parts of the theory to be developed will work also for time- 
dependent fields. As before, ( |1.29| ), ( |1.31| ) are to be solved as initial value 
problem. Thus E(x,0), B(x,0), q(0), v(0) are given. Note that the conti- 
nuity equation is satisfied by fiat. 

( |1.29p , (1.31) are the stationary points of a Lagrangian action, which 



strengthens our trust in these equations, since every microscopic classical 
evolution equation seems to be of that form. We introduce the electromag- 
netic potentials as in ( |1.24j ), ( |1.32| ). Then the action for ( |1.29j ), ( |1.31| ) reads 



S([q,<!>,A]) = ydt[-mo(l-g(t) 2 ) 1/2 - e (0 cx (9(t)) + 0(q(t),t) 

-q(t)-(A ex (q(t))+A(q(t),t)))] (1.33) 
+\J dt J d 3 % [(V0(aj, t) + d t A(x, t)) 2 - (V x A{x, t)) 2 ] . 
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The only difficulty is that ( |1.29| ), ( |1.31| ) make no sense mathematically. As 
explained, the solution to the Maxwell equations is singular at x = q(t) and 
in the Lorentz force equation we are asked to evaluate the fields exactly at 
that point. One might be tempted to put the blame on the mathematics 
which refuses to handle equations as singular as (|1.29|) , ( |1.31|) . However 
before such a drastic conclusion is drawn, the physics should be properly 
understood. The point charge carries along with it a potential which at 
short distances diverges as the Coulomb potential, cf. (|1.21|) , and which 
therefore has the electrostatic energy 



^ J d 3 xE(x,t) 2 ~ J drr 2 {r~ 2 ) 2 = J drr' 2 = oo . (1.34) 

{|iB-q(t)|<Jl} 

Taking literally such an object would have an infinite mass and cannot re- 
spond to external forces. It would maintain its velocity forever, which is not 
what is observed. 



Thus we are forced to regularize the Maxwell-Lorentz equation \1.29Q , 



( \1.31\) at short distances. 

To carry out such a program there are two in part complementary points 
of view. The first one, which we will not adopt here, regards the regular- 
ization as a mathematical device with the sole purpose to make sense out 
of a singular mathematical object through a suitable limiting procedure. To 
mention only one prominent mathematical physics example. The free scalar 
field, (f>(x), in Euclidean quantum field theory in 1+1 dimensions fluctu- 
ates so wildly at short distances that an interaction as J d 2 xV(<f)(x)) with 
V(4>) = <p 2 + A0 4 cannot be properly defined. One way, not necessarily op- 
timal, to regularize the theory is to introduce a spatial lattice with lattice 
spacing a. Such a lattice field theory is well defined in finite volume. One 
then carries out the limit spacing a — > at a simultaneous readjustment of 
the interaction potential, V(4>) = V a {4>) and obtains a Euclidean invariant, 
interacting quantum field theory. Ideally the limit theory should be indepen- 
dent of the regularization procedure. E.g. we could start with the free scalar 



18 



field in the continuum and also regularize (j>(x) as (fi*g(x) with g a test func- 
tion concentrated at 0. Then the regularized interaction is J d 2 xV((fi * g(x)) 
and in the limit g(y) i— > S(y) a quantum field theory should be obtained 
identical to the one of the lattice regularization. 

In the second point of view one argues that there is a physical cut-off 
coming from a more refined theory, which is then modelled in a phenomeno- 
logical way. While this is a standard procedure, it is worthwhile to illustrate 
it in a concrete example. We consider many (= 10 23 ) He 4 atoms in a con- 
tainer of adjustable size and we want to compute their free energy according 
to the rules of statistical mechanics. The more refined theory is here non- 
relativistic quantum mechanics which treats the electrons and nuclei as point 
particles carrying a spin |, resp. spin 0. As far as we can tell, this model 
approximately covers a temperature range T = °K to T = 10 5 °K, way 
beyond dissociation, and a density range p = to p = close packing. Beyond 
that relativistic effects must be taken into account. However there is a more 
limited range where we can get away with a model of classical point particles 
interacting through an effective potential of Lennard- Jones type. Once we 
specify this pair potential, classical statistical mechanics makes well-defined 
predictions at any T, p. There is no limitation in theory. Only outside a 
certain range the classical model looses the correspondence to the real world. 
Already from the way we describe the physical cut-off, there is a fair amount 
of vagueness. How much error do we allow in the free energy? What about 
more refined properties like density correlations? An effective potential can 
be defined quantum mechanically, but it is temperature-dependent and never 
strictly pair. Despite all these imprecisions and shortcomings, the equilib- 
rium theory of fluids relies heavily on the availability of a classical model. 

In the same spirit we modify the Maxwell-Lorentz equations by intro- 
ducing an extended charge distribution as a phenomenological model for the 
left out quantum electrodynamics. The charge distribution is stabilized by 
strong interactions which are outside the realm of electromagnetic forces. 
If the particle is an electron then on the classical level it looks like an ex- 
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tended charged object with a size roughly of the order of its Compton wave 
length, i.e. 4 x 10~ n cm. We impose the condition that the extended charge 
distribution has to be adjusted such that, over the range where classical elec- 
trodynamics is applicable, the coupled Maxwell-Lorentz equations correctly 
reproduce the empirical observations. 

Such general clauses seem to leave a lot of freedom. However, charge 
conservation and the derivability of the equations of motion from an action 
severely limits the possibilities. In fact, essentially only two models of an 
extended charge distribution have been investigated so far. 

(i) The semi-relativistic Abraham model of a rigid charge distribution. The 
^-function in ( |i.29| ), ( |1.3U| ), and ( |1.31| ) is replaced by a smooth function 
p(x) which is radial and vanishes for \x\ > R p . In Fourier space it means 
that couplings between the particle and Fourier modes with \k\ > R~ x are 
smoothly suppressed. The Abraham model will be studied in considerable 
detail . While defined for all velocities \v(t)\ < 1, it becomes empirically 
incorrect at velocities close to one. Despite this draw-back we hope that the 
Abraham model serves as a blue-print for more realistic cut-off prescriptions. 

(ii) The relativistic Lorentz model of a deformable charge distribution. The 
Abraham model violates relativistic invariance, since the choice of p singles 
out a specific reference frame. More in accord with relativity is to require 
that the charge distribution is radial in its momentary rest frame. We will 
discuss details in Section |1.5| and only remark already now that relativistic 
invariance forces the equations of motion to be nonlocal in time. 

We emphasize that for extended charge models the diameter R p of the 
charge distribution defines a length (and upon dividing by c also a time) 
scale, relative to which the approximate validity of effective theories, like the 
Lorentz-Dirac equation, can be addressed quantitatively. In fact, apart from 
the external forces, R p is the only length scale available. 
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1.4 The Abraham model 



Following Abraham, we model the charge as a spherically symmetric, rigid 
body to which the charge is permanently attached. To be specific the charge 
distribution, p, with total charge e is assumed to be smooth, radial, and 
supported in a ball of radius R p , i.e. 

Condition (C): 

p e C^(M. 3 ) , p(x) = p r (\x\) , p(x) = for \x\ > R p , J d 3 xp{x) = e . 

(1.35) 

To be definite, we require e > 0. With the obvious sign change, also e < is 
then covered. 

Equivalenty, we could introduce the form factor / such that 

p(x) = ef(x 2 ) , Jd 3 xf(x 2 ) = l. (1.36) 
The equations of motion for the Abraham model read then 
d t B(x,t) = -VxE(x,t), 

d t E(x,t) = V x B(x,t) - p(x - q(t))v(t) , (1.37) 
V-E(x,t) = p(x-q(t)), V-B(x,t) =0, (1.38) 

— (m b 7«(0) = / d3x P( x ~ <?W) [E cx (x) + v(t) x B cx (x) 

+E(x,t) +v(t) x B(x,t)] . (1.39) 

In contrast to the Lorentz force equation, for the Abraham model we denote 
the mechanical mass of the particle by m\, to emphasize that this bare mass 
will differ from the observed mass of the compound object "particle plus 
surrounding Coulomb field". For the external potentials <f) cx , A cx we have 
considerable freedom. We require them to be smooth and locally bounded, 
including their derivates to avoid too strong local oscillations. No condition 
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on the increase at infinity is needed, since \v(t)\ < 1. However, it is con- 
venient to have the energy (|1.42|) uniformly bounded from below. To keep 
matters simple we make the unnecessarily strong assumptions 

Condition (P): 

<f) cx G C°°(M 3 ) , A cx G C°°(M 3 ,M 3 ), CX > > -oo . (1.40) 

There exists a constant C such that |V fc ex | < C, |V fc A cx | < C component- 
wise for k = 1,2, 3. 

Physically, the Abraham model is not quite consistent. Besides the center 
of mass q(t) and its velocity v(t) a rigid body has also rotational degrees 
of freedom. As we will see in Chapter [5] even if initially non-rotating the 
external fields and the self-interaction necessarily induce a rotation. The 
translational and rotational degrees of freedom are coupled through the field. 
This makes the model considerably more intricate and it is an advisable 
strategy to understand the simplified version first. 

The Abraham model is derived from the Lagrangian 

L = -m b (l - q 2 ) 1 ' 2 - (0 CX + cj) — q ■ A cx — q ■ A) * p(q) 

+ l - Jd 3 x [( V0 + d t A) 2 - ( V x A) 2 ] . (1.41) 
Correspondingly the energy 

S(E, B, q, v) = m hl (v) + e0 ex * p (q) + ± J d 3 x (E(x) 2 + B{x) 2 ) (1.42) 
is conserved. 

As for any dynamical system, we first have to construct a suitable phase 
space. The dynamical variables are (E(x), B(x), q,v) — Y which is called a 
state of the system. We have q G M 3 , v G V = {v : | v\ < 1}. In addition the 
field energy (|1.42j ) should be bounded. Thus it is natural to introduce the 



(real) Hilbert space 

L 2 = L 2 (M 3 ,R 3 ) (1.43) 
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with norm ||JS|| = (/ d 3 x\E(x)\ 2 ) 1 ' 2 and to define L as the set of states 
satisfying 

\\Y \\ c = \\E\\ + + \q\ + \l{v)v\ < oo . (1.44) 

In particular the field energy, |(||.E|| 2 +||.B|| 2 ),is bounded. We equip C with 
the metric 

d(Y 1 ,Y 2 ) = \\E 1 - E 2 \\ + \\B 1 - B 2 \\ + | 9l - q 2 \ + \-f{v 1 )v 1 - -f{v 2 )v 2 \ . 

(1.45) 



In addition, one has to satisfy the constraints ( |1.3S| ). Thus the phase space 



A4, for the Abraham model is the nonlinear submanifold of £ defined through 

V ■ E(x) = p(x-q), V-B(a;) = 0. (1.46) 

M. inherits its metric from C. 

At several occasions we will need that the system forgets its initial field 
data. For this purpose it is helpful to have a little bit of smoothness and 
some decay at infinity. Formally we introduce the "good" subset M"cM, 
< a < 1, which consists of fields such that componentwise and outside 
some ball of radius Ro, \x\ > Ro, we have 

\E(x)\ + \B(x)\ + |aj|(|V£7(sc)| + |VB(aj)|) < C {x^ 1 -* . (1.47) 

The Lienard-Wiechert fields (|1.21|) , (|1.22|) are in M. 1 and M 1 is dense in M . 
However A4 a = for a > 1, since J d 3 xp(x) = e ^ 0. 

The evolution equations (|1.37f) to (|1.39| ) are of the general form 



jY{t) = F{Y(t)) (1.48) 
with Y(0) = Y° G M.. We have to turn to the question of the existence and 



uniqueness of solutions of the Abraham model (|L48|). 
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Theorem 1.1 Let conditions (C) and (P) hold and let Y° = {E°(x), B°(x), q°, 
v°) G M.. Then the integrated version of Equation (\1.4fy) , 

t 

Y(t)=Y° + J dsF(Y(s)), (1.49) 


has a unique solution Y(t) = (E(x,t), B(x,t), q(t),v(t)) G M., which is 
continuous in t and satisfies Y(0) = Y°. Along the solution trajectory 

£(Y(t)) = £(Y°), (1.50) 

i.e. the energy is conserved. 

For short times existence and uniqueness follows through the contraction 
mapping principle with constants depending only on the initial energy. For 
smooth initial data energy conservation is verified directly and by continuity 
it extends to all finite energy data. Thus we can construct iteratively the 
solution for all times. 

We first summarize some properties of the inhomogeneous Maxwell equa- 
tions. They follow directly from the Fourier and convolution representations 
(|ng), flTTJ), resp. QTTTBp, (|ET7p. 



Lemma 1.2 In the inhomogeneous Maxwell equations jjl-dj ) let p(x, t) 



p(x — q(t)), j(x,t) = p(x — q(t))v(t) with t 1— > q(t),v(t) continuous. Then 
has a unique solution in C{R,L 2 © L 2 ). The solution map 
(E°,B°) 1 — > (E(t), B(t)) depends continuously on q(t),v(t). 

Proof of Theorem 1.1: Let us fix some b > and choose initial data such 
that £(Y°) < b. 

(i) There exists a unique solution Y(t) G C([0, e), M) for e = e(b) sufficiently 
small. 

We write ( |1.42|) in the form 



j t {m hl v{t)) = F cx (t) + F (0) (t) + Fscif(t) (1.51) 



24 



by inserting E(x,t), B(x,t) from the Maxwell equations according to ( |1.16 ), 
flTTTp . Let 

W t (x) = [ d 3 k\p(k)\ 2 e ik - x -]- sin I felt 
7 |fe| 

= (2rr) 3 Jd 3 y Jd 3 y'p(y)p(y')^- t 5(\y + x-y'\-t). (1.52) 

Then 

F M (t) = E ex * pit) + v(t) x B cx * p(t) , (1.53) 

F (0) (t) = y d 3 zp(a; - * E°(x) + (V x G t ) * £ o (a0 

+«(*) x d t G t * B°(aj) - w(t) x(VxG t * E°(x))] ,(1.54) 

t 

Fseif(t) = J ds[-VW t - a {q{t) - q(s)) - v(s)d t W t - 8 (q(t) - q(s)) 



o 



+v(t) x (V x v(s)W t „ s (q(t) ~ q(s)))\ ■ (1-55) 

We now integrate both sides of (|1.51|) over the time interval [0, t}. The re- 
sulting expression is regarded as a map from the trajectory t i— > qf(t), < 
£ < 5, to the trajectory t i— ► u(t) and is defined by 

t 

= q° + J dsv(s), (1.56) 



o 



m b7 (v(t))v(t) = m b7 (t;> + y ds(F ex (s) + F {0) (s) + F self (s)) 



o 



where F fflt (s), F( )(s), F se if(s) are functionals of q(-), f (') according to (|1.53|) 
to ( |1.55| ). Since p, W, </> ex , A ex are smooth, this map is a contraction in C([0, t], 
R 3 x V), i.e. 

sup (\q 2 (s) - qx(s)\ + \v 2 (s) - «i(s)|) 



0<s<t 



<c(t,6) sup (|9i(s)-qf 2 (s)| + |«i(s)-« 2 (s)|), (1.57) 



0<s<t 
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with a constant c(t, b) depending on b and c(t, b) < 1 for sufficiently small t. 
Such a map has a unique fixed point which is the desired solution q(t), v(t). 
By the Maxwell equations also E(x,t),B(x,t) is uniquely determined. 

(ii) The solution map Y° *—>■ Y(t) is continuous in M.. 

This follows from Lemma |1.2| and the continuity of q(t),v(t) in depen- 
dence on the initial data. 

(iii) The energy is conserved. 

We choose smooth initial fields such that E, B G C°°(1R 3 ) and 

\V a E(x)\ + \V a B(x)\ < C{1 + |x|)- (2+|a|) . (1.58) 

Here a = (01,02,03) is a multi-index with aij = 0, 1,2,... . This set is 
dense in Ji4. By the convolution representation (|1.16|) , ( |1.17|) of the solution 
to the Maxwell equations we have E(x,t), B(x,t) G C 1 ([0,<5] x M 3 ) and 
\E(x,t)\ + \B(x,t)\ < C{1 + \x\Y 2 . Also v(t) G ^([0,5]). Thus we are 
allowed to differentiate, 

- E(Y(t)) = ^vv + v V0 ex * p(q) + / d 3 x(E -d t E + B- d t B) 

= Jd 3 x(E- (V x B) - B ■ (V x E)) = 0, (1.59) 

since the fields decay and hence the surface terms vanish. Thus S(Y(t)) = 
S(Y°) for < t < 5. By continuity this equality extends to all of A4. 

(iv) The global solution exists. 

^From (iii) we know that S(Y(S)) = S(Y°) < b. Thus we can repeat the 
previous argument for 5 < t < 25, etc.. Backwards in time we still have the 
solution ( 1.16 ), ( 1.17|) of the Maxwell equations, only the retarded fields have 
to be replaced by the advanced ones. Thereby we obtain the solution for all 
times. □ 



Theorem |T7T| ensures the existence and uniqueness of solutions for the 
Abraham model. For initial data Y° G M. the solution trajectory t — > Y(t) 
lies in the phase space M., is continuous in t, and its energy is conserved. We 
have thus established the basis for further investigations on the dynamics of 
the Abraham model. 
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1.5 Appendix: Long time asymptotics 

For dynamical systems one of the first qualitative question is to understand 
whether there are general patterns governing the long time behavior. For 
the Abraham model the long time asymptotics is dominated through the 
loss of energy radiated to infinity, which is proportional to v{t) 2 according to 
Larmor's formula. Since the energy is bounded from below, we expect that 

lim v{t) = (1.60) 

t— *oo 

under rather general initial conditions. In fact, one would also expect that 

lim v(t) = G V , (1.61) 

t— »OC 

where = for bounded motion and ^ for a scattering solution. 
In this section we will prove ( 1.60|) under the extra hypothesis 

Wiener Condition (W): 

p(k)>0. (1.62) 

The proof follows rather closely the physical intuition and leads to an equa- 
tion of convolution type which has a definite long time limit only under (W). 
(W) means that all modes of the charge distribution couple to electromag- 
netic field. According to (|1.60|) the Abraham model does not admit then any 
periodic solution. Since p has compact support, in general, p(k) may vanish 
for a discrete set of shell radii |fe|. At present, it remains as an open problem, 
whether periodic solution become then possible. 

Let us consider a ball of radius R centered at the origin. At time t the 
field energy in this ball and the mechanical energy of the charge is given by 

e R (t) = e{0)-± j d 3 x(E(x,t) 2 + B(x,t) 2 ) (1.63) 

{\ X \>R} 

for R sufficiently large, using the conservation of total energy. Sr changes in 
time as 

A g R (t) = _i?2 J d 2 U0J . [E(Ru, t) x B(Ru>, t)} , (1.64) 
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where a; is a vector on the unit sphere, cPu the surface measure normalized 
to Att, and ExB the Poynting vector for the flux in energy at the surface of 
the ball under consideration. Since the total energy is bounded from below, 
we conclude that 

R+t 

/d 
ds — 8 R {s)<C (1.65) 

R 

with the constant C = 8(0) — (f> independent of R, t. 

In (|1.65|) we first take the limit R — > oo, which yields the energy radiated 
to infinity during the time interval [0, t] through a large sphere centered at 
the origin. Subsequently we take the limit t — > oo to obtain the total radiated 
energy. To state the result we define 

1 r„ 



Eoofat) = -— J d 6 y p(y - q(t + u> ■ y)) (1.66) 

[(1 - U? ■ V)' 1 ^ + (1 - OJ ■ v)~ 2 (u> ■ V)(V - Cl>)] \t+u.y 

which is a functional of the actual trajectory of the particle. Whatever its 
motion we must have 

oo 

J dt J d 2 uu\E^(u,t) | 2 < C < oo. (1.67) 



Note that the integrand in Ql.67 ) is proportional to v{t) 2 , which therefore 
has to decay to zero for large t. 



To establish ( 1.67 ) is somewhat tedious with pieces of the argument 



explained in Sections £0| and |6.4| , |6.5| . One imagines that the trajectory 



t i — > q(t) is given and solves the inhomogeneous Maxwell equations accord- 
ing to ([L. 16 ), (|1.17 ). If the initial fields are in M a , \ < a < 1, then they 



decay and make no contribution to (|1.65|) in the limit R — > oo, cf. our 
treatment of the initial time slip in Section [4.3| . Next one has to study the 
asymptotics of the retarded fields, which is carried out in Sections |6.4| , |6.5| . 
There e is fixed and for our purpose we may set e — 1. In addition in ( |6.49 ) 
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the sphere is centered at q £ {t), rather than at the origin. This means we can 
use ( |6.52| , ( |6.53|) with q £ (t) replaced by for our case. 



The real task is to extract from ( |1.67|) that the acceleration vanishes for 
long times. 

Theorem 1.3 For the Abraham model satisfying {C), (P), and the Wiener 
condition (W) let the initial data Y(0) = (E°, B°, q°, v°) e M u with \ < 
a < 1. Then 

lim v(t) = 0. (1.68) 

t— >oo 



Proof. By energy conservation \v(t)\ < v < 1. Inserting in ( 1.39|) and using 



(P) we conclude that \v(t)\ < C. Differentiating (|1.39|) and using again (P) 
also \v{t)\ < C uniformly in t. Therefore E^ioj^t) is Lipschitz continuous 
jointly in u>,t. Since the energy dissipation ( |1.67| ) is bounded, this implies 

lim E^uj.t) = (1.69) 

uniformly in u). 

We analyze the structure of the integrand in (|1.66|) . In the retarded 
argument only y» — u; • y appears. Therefore the integration over y ± = 
y — y\\uj can be carried out and we are left with a one-dimensional integral 
of convolution type. We set p a {xz) = J dx\dx-i p(x). Then 

Eoo(w,t) = ±Jdnto(M-m(t+w)) ( L7 °) 

[(1 - w ■ v)- 2 u x ((w - v) x v)] | t+2/|| 
= ^ J dsp a (t - (s - q n (s)))[(l - u ■ v)~ 2 u x ((w - «) x v)]\ s . 

Since < 1, we can substitute 9 = s — q\\(s) and obtain the convolution 

representation 

E 00 (u,t) = [ d9 Pa (t-9)g w (6)=p a *g w (t), (1.71) 
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where 



9M = i- [(1 - w • v)- 2 u x ((w - v) x «)] | s(e) . (1.72) 



From ( |1.69|) we know that lim p a * g w {t) = 0. If p(fc ) = for some fc , 



and hence the Wiener condition would not be satisfied, then the convolution 
integral admits a periodic solution and no further progress seems to be pos- 
sible. However with (W) and the smoothness of g w (9) already established, 
Pitt's extension to the Tauberian theorem of Wiener ensures us that 



lim g„(e)=0, (1.73) 

which, since 9(t) — > oo as t — > oo, implies 

lim«x((w-«(i))x«(t)) = (1.74) 

for every u) in the unit sphere. Replacing u> by —uj and summing both 
expressions yields u> x (w x — > as t — ► oo. Since this is true for every 
w, the claim follows. □ 

Next we study the comoving electromagnetic fields for large times. As 
used already, under our assumptions the initial fields decay as t — ► oo. Thus 
we only have to consider the true retarded fields E iet (x + q(t),t), B rct (x + 
q(t),t) centered at the position of the particle. Since v(t) — > 0, these fields 
become almost stationary. In Section |3.1| we show that there is a unique 
comoving field with velocity v. These charge soliton fields are denoted by 
E V ,B V when centered at the origin, compare with (|3.21 ), (|3.22| ). Thus the 



true retarded fields E ret (x + q(t),t), B ret (x + q(t),t) have to compared with 
the soliton fields E v (t),B v (t)- F° r this purpose we use the representations 
(HOI , (|3~22|) for the charge soliton and ([i~T6|) , ( |i~T7|) for the retarded fields. 
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We insert the explicit form ( |1 . 1 5| ) of the propagator. This yields 

E v (x) = j d?y{kK\x-y\)~ l (\x-y\~ l p{y - v\x - y\)n 

+v ■ Vp(y — v\x — y\)(v — n)) , (1-75) 

B v (x) = J <i 3 ?/(47r|a3 — y|) _1 n x ( — \x — y\~ x p{y — \x — y\v)v 

+v ■ Vp(y — \x — y\v)v) , (1-76) 

where n = (x — y)/\x — y\. Similarly for the retarded fields 

E Tet {x + q{t),t) = J d z y{A-K\x-y\)- l (\x-y\~ l p{y + q{t)-q{T))n 

+v(t) ■ Vp{y + q(t) - q(r))(v(r) - n) 

- p (y + q(t)-q(r))v(r)), (1.77) 

B Ict (x + q(t),t) = J <Py{An\x - y\)- l n x (- \x - y\~ x p{y + q{t) - q{r)) 

v(t) + v(t) ■ Vp(y + q(t) - q(r))v(r) 

-p( y + q (t)-q(r))v(T)), (1.78) 

where r — t — \x — y\ and t >t p = 2R P /(1 — v). 

We compare the fields locally and use that lim v(t) = 0. Then, for any 
fixed R > 0, 



lim J d 3 x((E{x + q{t),t)-E v{t) {x))'' 



(B(x + q(t),t)-B v(t) (x)) 2 ) =0. (1.79) 



t— >oo 

{\ X \<R} 



Armed with this information we can sketch qualitatively the long time 
behavior of the Abraham model. We distinguish bounded and scattering 
trajectories. If ex and A ex decay sufficiently fast for large \x\ and if the 
particle escapes into the essentially force free region, then lim v(t) = ^ 
and the fields are well approximated by the soliton fields, compare with 



( 1.7S ). In particular, if there are no external fields the charge travels with 
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some definite velocity in the long time limit. This point is discussed in more 



detail by a different method in Section [O 
If the motion is bounded, 

\q(t)\<q (1.80) 

for all t, then v(t) — > for t —>■ oo implies that 

limv(*) = 0. (1.81) 

Inserting in the Lorentz force equation ( |1.39| ) and using that the fields become 
soliton-like, one infers that 

lim V0 ex *p(q(t)) = 0, (1.82) 

i.e. q(t) approaches the set of critical points A = {q, V0 ex * p(q) = 0}. 

( |1.81| ). ( |1.82| ) still leave a lot of freedom. Generically, one expects that 
q(t) approaches a definite limit. This is indeed the case if A happens to be 
a discrete set. By continuity of solutions in t, q(t) has to converge then to 
some definite q* G A. In particular if ex * p is strictly convex, the charge 
will come to rest at the minimum of ex * p- HA is not discrete, the long 
time behavior depends on be the specific situation. E.g. let A ex = and 
CX be strictly convex outside a ball of radius one and let ex = inside this 
ball. Each time the particle is reflected by the confining potential it looses in 
energy. Thus Q1.81 ) holds, but q(t) has no limit as t — > oo. A more realistic 



example is a constant magnetic field and ex = 0. Then A = M 3 , but using 
the sharper estimates of Section |4.3j it can be shown that the particle spirals 
inwards to come to rest at its center of gyration. 

So far ( |1.80| ) was an assumption and it would be nice to have some suffi- 



cient criteria. It seems that the only one available is conservation of energy. 
Bounded motion requires then sufficiently deep local minima of ex . Thus we 
see that bounded energy dissipation yields reasonable results for electrostatic 
forces but is rather weak for motion in magnetic fields, except for the general 
fact that the acceleration has to vanish for large t. 
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In general, bounded and unbounded motion coexist. The prime example 
is a charge bound by an infinitely heavy nucleus, which is modelled by an 
attractive, locally smoothened Coulomb potential. If the initial field supplies 
sufficient energy, then the "atom" becomes ionized and the charge will travel 
freely in the long time limit. Conversely a charge may loose energy through 
radiation and become trapped by the external potential. 

1.6 The Lorentz model 

To improve on the semi-relativistic Abraham model, with Lorentz it is natu- 
ral to assume that the charge distribution is rigid in its own rest frame. The 
actual construction of such a charge distribution requires some effort. For 
obvious reasons we will switch now to a relativistic notation. 

We fix a laboratory frame, S, with coordinates (t, x) = x M , /i — 0, 1, 2, 3, 
and invariant length t 2 —x 2 = t 2 —XiXi = x^x^ = g [iV x ll x v ', where we adopt the 
standard summation convention and use latin indices for three-vectors, greek 
indices for four-vectors. In S we prescribe the world line of a charged particle 
t i-> Zi(t),i = 1,2,3, with velocity Vi(t) = ^Zi{t), v(t) 2 < 1. We will also 
parametrize the world line through its eigentime as s i— > z^(s), /i — 0, 1, 2, 3. 

To construct the extended charge distribution we consider some time 
t — to when the particle has velocity Vi = Vi(to). The rest frame for the 
particle at time to is denoted by S' with coordinates x ,fl . It is connected to 
S by the Lorentz transformation 

x ,lx = \{vy v x v . (1.83) 

We require the origins of S and S' to coincide and their spatial axes to be 
parallel to each other, i.e. A(v) is a Lorentz boost and reads explicitly 

x'i = Xi + (j 2 / {I + l))viV k x k -<yvit, 
t > = ^t-v^). (1.84) 

The Lorentz model assumes that in S' the four-current j' M is given by 

J" = e(/((4 - WWi - 4(0)), 0)|^ (1-85) 
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with the relativistic form factor /, where ef(x 2 ) = p r (\x\) as before. In 
particular, / vanishes for \x\ > R p . 

We have to transform j 1 ^ back to the laboratory frame S. Using (|1.84 ), 
(|i~85| ) yields 



f = e( 7 ,7^)/((4 - ztfWj ~ **(0))lf=* • (L86) 
The condition t' = t' means j(t — v^x^) = 7(^0 — VkZk{to)), i.e. 

*o = t - v k (x k - Zk{t Q )) . (1.87) 

Therefore 

f(t,x) = e(7, 7 -y)/((t - t ) 2 - - ^(t ))(^ - ^-(t ))) (1-88) 

with t = t (t, x) defined through ( p..87| ). 

An example is shown in Figure 1. x) vanishes outside a tube around 
the world line. This tube is fibered into cross sections corresponding to 
{f = t' } in the momentary rest frame. Along each cross section the charge 
is smeared according to the form factor /. From the figure we notice that, 
when the acceleration becomes too large, the cross sections overlap, which 
means that for given (t, x) Equation ( |1.87p has several solutions. ( |1.88|) 



seems to indicate that one should add the contribution from each solution. 
This however would violate charge conservation and the proper prescription 
is to reverse the sign of a charge element when it moves backwards in time. 
Taking into account multiple solutions to ( 1.87| ) and their proper sign leads 
to the four-current 



f(x)= J dsv^l-v^x- zY]f((x- zy)5(v x (x- zf) (1.89) 

—00 

with x = (x°, . . . , x 3 ). Here z^ = z**(s) is parametrized by its eigentime s and 
r M (s) = 4-z^(s). The current (|1.89|) satisfies the charge conservation 

d^f = 0. (1.90) 
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There is an instructive way to rewrite the four-current by using the 
Thomas precession, which we recall first. The elementary observation is 



that the Lorentz boosts, cf. ( |1.84|) , do not form a subgroup of the Lorentz 
group. Let us consider the inertial frame S' with velocity v' relative to S 
and the inertial frame S" with velocity v" relative to S' and thus velocity 
w, the relativistic sum of v and v', relative to S. If we denote by A(v) the 
Lorentz boost with velocity v, then 

A(v')A(v) = R(v, v')A(w) , (1.91) 

where R(v, v') is a suitable pure spatial rotation. 

To apply this observation to a relativistically rigid charge distribution we 
imagine that it has a body fixed frame. We want to define an inertial frame 
/C(s) such that relative to this frame the body axes maintain their orientation 
throughout time. We agree that the time axis of JC(s) is parallel to the four- 
velocity t> M (s) and that the origin of S and JC(s) coincide. Let A(s) be the 
Lorentz transformation from S to K,(s). To determine it we subdivide the 
time axis into little intervals of length As and require that JC((m + l)As) is 
related to /C(raAs) by a Lorentz boost with the properly adjusted velocities, 
i.e. the space axes of JC((m + l)As) and /C(mAs) are parallel to each other. 
Taking the limit As — > one finds that A(s) is determined through 

A^ = V \A\, A^(0) = A">(0)) , (1.92) 

where 



1.93) 



Using this definition and noting that v v £ u 8(v\£ x ) = the four-current ( |1.89|) 
equals 



fix) 



ds[v» - rf„(x - zY\f{{x - zf) S(v x (x - z) x ) 



1.94) 



35 



in which form it has a transparent physical meaning. The current is the 
sum of two contributions. There is a translational part proportional to v M 
as expected from the nonrelativistic limit. In addition there is a rotational 
component due to the Thomas precession. The minus sign is a convention. 
With our definition a vector time-independent in JC(s) has components 
in S and they change as 

x il = -rf v x v . (1.95) 

Before proceeding to the action for the dynamics, we should understand 
whether the current ( 1.89|) conforms with the naive physical intuition. An 



instructive example is a uniformly accelerated charge, the so-called hyper- 
bolic motion. We assume that the particle is accelerated along the positive 
1-axis starting from rest at the origin. In the orthogonal direction the current 
traces out a tube of diameter 2R P and it suffices to treat the two-dimensional 
space—time problem. The center, C, of the charge moves along the orbit 

C = (t,g- 1 (y/T+g~W-l)), t>0, (1.96) 

where g > is the acceleration. The curves traced by the right and left ends, 



C + and C_, are determined from ( 1.84|) and are given in parameter form as 



C ± = ((1 ± R p g)t, g- l {{l + R p g)y/l+gH*- 1)), t > . (1.97) 

The equal-time distance between the center and C + is t~~ 1 ( y (R p g) 2 +2R p g) / {2g 2 
(1 + R p g)) for large t and thus well localized. However the left end motion 
depends crucially on the magnitude of R p g. If R p g < 1, then the distance 
to the center is t~ l ((R p g) 2 — 2R p g)/{2g 2 (l — R p g)) for large t. On the other 
hand for R p g > 1, the left end moves into the past, cf. ( |1.97| ), and the 



current density looks strangely distorted. To have a feeling for the order of 
magnitudes involved we insert the classical electron radius. Then 

c 2 

g > — = 10 31 [raj sec 2 ] , (1.98) 
R p 



36 



which is huge and way beyond the domain of the validity of the theory. Of 
course, one would hope that for reasonable initial data such accelerations 
can never be reached. But the mere fact that charge elements may move 
backwards in time is an extra difficulty 

Armed with the four-current ( |1.94| ) we can write down the relativistically 
covariant action, S, as a functional of the the four-potential A 11 and the 
particle coordinates z^ . The action of the Lorentz model has four pieces, 

S = S + S f + S int + S cx . (1.99) 

So is the mechanical action, 

*2 



S = m b / ds. (1.100) 



Sl 



Sf is the field action, 



S f = -~ J d*xF^(x)F^(x) , (1.101) 



where the field tensor F^ v is defined through 

= d tM A" -d u A». (1.102) 
The interaction is bilinear in the current and the field, 

Stat = J d*xA»{x) 3 ^x) , (1.103) 
h 

and correspondingly for the interaction with the external potentials 

d 4 xA^(x)j^x). (1.104) 



•"ex 

n 



The world line z^{s) is specified for s\ < s < S2- Consequently Q is the 
volume in S bounded by the two hyperplanes v^i^s^^x — z(si))^ = and 
u"(s 2 )(x-z(s 2 )V = 0. 
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The actual dynamical trajectory is a stationary point of the action S at 
fixed endpoints. The variation with respect to the potentials leads to the 
inhomogeneous Maxwell equations 

d v F v,l = f. (1.105) 

The constraints ( |1.3| ) are automatically satisfied, since the field tensor is 
derived from potentials according to ( |1.102| ). Next we vary the world line z^. 
From So we obtain the mechanical acceleration m^v 11 . The variation of S; nt is 
somewhat lengthy and deferred to Section |1.7| . The final evolution equation 
reads 

ro b ^(«) = eJd 4 U(e)S(v x (s)e)[F^(z(s)+0 + F^(z(s) + 0) 

m<o - T^oor) ■ (i-loe) 

At first glance the equations of motion ( |1.10(jp look rather similar to its 
semi-relativistic sister. It seems natural to specify then z(0),z(0) = v(0) 
and the field tensor on the hyperplane determined by (x — z(0)) x v\(0) = 0. 
If we assume that R p \v\ < 1, then the future is decoupled from the past 
and one would hope to have a unique solution. For short times one can 
presumably copy the proof given for the Abraham model. We are not aware 
of mathematical results which ensure the existence global in time. If the 
acceleration becomes large and R p \v\ > 1 for some time span, it is not even 
clear how to properly specify the initial data. Since such questions remain 
largely unexplored, the only option is to proceed as if a solution is well 
defined. 

As seen from the laboratory frame S the rigid charge distribution Thomas 
precesses. To be physically consistent we have to allow then also for a rotation 
of the body fixed frame relative to JC(s). For sure, such a rotation will 
be induced through the back reaction of the electromagnetic field onto the 
relativistically rigid charge. Let us denote by 

u)» v (s) (1.107) 
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the angular velocities of the motion of the body fixed frame relative to /C(s). 
From the point of view of the laboratory system they simply have to be added 
to rj fl u , i.e. in the current Ql.94|) Tf" v has to be substituted by 77^ + u^ u . The 
variation of the action is now with respect to A^, z M , and the Euler angles 
of the body fixed frame relative to JC(s). This results in the inhomogeneous 
Maxwell equations, as before, and in a coupled set of translational and rota- 
tional equations of motion which are nonlocal in time. 



1.7 Appendix: Variation of the action 



We carry out the variation <5Si nt of ( |1.103| ). Let us first list the necessary 
identities. Since {ds) 2 = dz^ dz^ we have 

dz a d 
5(ds) = —Sdz a = v a ds — Sz a (1.108) 
ds ds 



and 



dz» dzK + 5dz» dz» d d 

8v — 8 — - — = — : — — = — bz^ -v**v a — 5z a . (1.109) 

ds ds + ods ds ds ds 



Similarly 



5v»= 5 v »-v» Va — 5z a . (1.110) 
ds ds 



In the variation of r] fl u we have to maintain the frames at z^(s) and z^(s) + 
Sz^^s) without relative spatial rotation. This leads to 

Srf = v »5v u - v v 8v„ =v»£- 8v v -v v ^- 8v^ - if v v a 5z° . (1.111) 

ds as as 

If in ( Llll) we would also vary v^, then the constraint of zero spatial rotation 
is no longer satisfied. 

In the following we will assume that ( |1.87 ) has a unique solution within 
the tube traced out by the form factor. For this we need v u (x — z) v < 1 inside 
the tube, equivalently R p \v\ < 1, i.e. the acceleration has to be sufficiently 
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small on the scale of 1/R p . Using this condition, S int can be written more 
explicitly as 



S2 

Sint, = e 



3 int 

si 



JdsJ ^ /(e 2 )^* W«) + W«) - v WO ■ (i-ii2) 



Then 



si 

+ejdsj d 4 U(e)5(v x e)A^z + 



si 



S2 

+ejds j d'U\e) [6(v x e)Sz a d^A^z + 

si 

+iW>^^ + 0]("/.-Vf)- (1-113) 



We define 
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Then 



si Si 
S2 

+ejdsj d 4 U(envxe)A»(z + 0^6z„ 

Si 
S2 

si 
S2 

+e j ds5z,j d^myivx^F^iz + Oiv.-Vuan 

Si 
S2 

+ e J dsdz.J rf 4 e/(e^K A )^(2+e)K-ve ff ). (1.115) 



Si 



The last term can be rewritten as, using xS(x) = 0, 



d 4 i f{e) Kvxi X ) d" A»{z + (v v - vu,n 
= - I d 4 £A»(z + ST [fit 2 ) 6(v x £ x )(v v - v v v a n\ 

= I d 4 tA»{z+i)f{e)5{v x e)v„v v 
- j d 4 i jxe) s\v x e) a»(z + - van 

= - j d 4 xf((x-z) 2 )5'(v x (x-z) x )A»(x)(l-v a (x-zy) 
= ^ J 'd 4 xf((x-z) 2 )6(v x (x-z) x )A»(x) 

j d 4 U(e)S(v x e)A^z + (1.116) 



ds 
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which cancels against the second term of ( 1.11 5| ) . Therefore 



S2 S 2 



SS int = - Jds-^(Q» SvJ + J ds 8z» (-£(Q»- v»v v Q v ) 



SI Si 



Under variation of 5z fl (s) the term in the big brackets has to vanish. 

The term containing F^ u is already in its final form. The remaining 
summands are the time derivative of /, where 



(1.118) 

We have, using partial integration with respect to £, x8(x) = 0, and xS'(x) = 
-8(x), 



= ejd 4 U(e)6(vxe)(A v (z + 0v u ^ + d*A»(z + Z)v a v v ^) 

+e f d^f{e)5'{v x e)A v {z + 0^ a v^ 
= e J d^f{e)6{v x e)A v {z + 0{v v ^-v»v v ) 

-e j d^f{e)8'{v x e)A v {z + 0v v ^{l-v a ?), (1.119) 
VuQT = -eJd 4 U(e)S(v x e)A v (z + 0v v v^, (1.120) 

q^-i/w = e j d 4 u(e)t(v x e)A»(z+t)v u e 

-e j d^f{e)8{v x e)A v {z + ti)v v v»{l-v a (i°) 
-e f ^/(0%fV"(2 + eKf(l-^)- (1.121) 
42 



We turn to the second term in ( |1.118| ) and use x8'(x) = —S(x), x 2 5'(x) = 0, 

+e j d^f{e)5{v x e)A v {z + t){v v v»-v v v»v^-v v e)- (1-122) 

Thus 1 = 0. 

We conclude that the variation 

5(S + S int )=0 (1.123) 

at fixed endpoints leads to the Lorentz force equation 

mh Ts V " = 6 J d^fi^i^F^iz + O^-v^n- (1-124) 



Notes and References 



ad ad l.l: The material discussed can be found in any textbook. Par- 



ticularly useful I find Jackson (1999) and Scharf (1994). 

ad \1 . 3j : In our history chapter, Section [|, we discuss other approaches which 
cannot be subsumed under short distance regularization. In the literature the 
size of a classical electron, r c , is usually determined through equating the rest 
mass with the Coulomb energy, m e c 2 = e 2 /r c , which gives r c = 3 x 10 _13 cm. 
This is really a lower bound in the sense that an even smaller radius would 
be in contradiction to the experimentally observed mass of the electron (as- 



suming a positive bare mass, cf. the discussion in Section [4.5|) . Milonni 
(1994) argues that due to quantum fluctuations the electron appears to have 
a classical extent, which is given by its Compton wave length A c = r c /a with 
a the fine structure constant. Renormalization in Euclidean quantum field 
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theory is explained by Glimm, Jaffe (1987). Effective potentials for classical 
fluids are discussed, e.g., in Huang (1987). 



ad \1.4\ : The Abraham model was very popular at the beginning of this cen- 
tury and was studied by Lorentz (1892,1915), Abraham (1903,1905), Som- 
merfeld (1904,1905), and Schott (1912), amongst others. Apparently a rotat- 
ing rigid charge had been considered only much later Frenkel (1926), Bhabha, 
Corben (1941). The proof on the existence and uniqueness of the dynamics 
is taken from Komech, Spohn (1999), where a much wider class of external 
potentials are allowed. A somewhat different technique is used by Bauer, 
Diirr (1999). They cover also the case of a negative bare mass and discuss 
the smoothness of solutions in terms of smoothness of initial data. 



ad [77||.- The long time asymptotics is treated in Komech, Spohn (1999), 
where the details of the proof can be found. Pitt's version of the Wiener the- 
orem is proved in Rudin (1977), Theorem 9.7(b). We remark that Theorem 
1.3 gives no rate of convergence. Thus to investigate the asymptotics of the 
velocity and position always requires some extra considerations. Komech, 
Spohn (1998) study the long time asymptotics for zero external potentials 
in the case of a scalar wave field. In particular, by these methods one can 
handle the asymptotics of scattering trajectories. Presumably these results 
extend to the Maxwell field. They require the Wiener condition. In Section 
[4.3| we give a different proof without Wiener condition (W) but for a suffi- 
ciently small charge. This indicates that (W) is an artifact of the proof. One 
might wonder, how bounded energy dissipation works in the case of several 
particles. The result is somewhat disappointing. One concludes only that 
the center of mass acceleration vanishes in the long time limit. To show that 
also the relative motions come to rest requires novel techniques. 

In the literature, Bohm, Weinstein (1948), in particular the review by 
Pearle (1982), periodic solutions of the Abraham model have been reported 
repeatedly for the case of a charged sphere, i.e. p(x) = e(47ra 2 ) _1 <5(|a;| — a), 
which is not covered by Theorem [0| since (W) is violated. These compu- 
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tations invoke certain approximations and it is not clear whether the full 
model, as defined by (|1.37|) to Q1.39Q , has periodic solutions. Kunze (1998) 



excludes periodic solutions in some small part of phase space for the scalar 
field without the Wiener condition (W). 



ad \1 . 6 : This section is based on the monumental work of Nodvik (1964). 



Nodvik includes the rotation of the body fixed frame. A complementary dis- 
cussion is given by Rohrlich (1990), Chapter 7-4. The relativistic Thomas 
precession is discussed in M0ller (1952). 

Of course relativistic theories have been studied much before, e.g. Born 
(1909). One difficulty is to write down the proper variant of the Lorentz force 
equation, which may be circumvented by considering the Lienard-Wiechert 
fields generated by a point charge and to infer through the balance of mo- 
mentum the actual motion of the charge, e.g. Rohrlich (1990), Teitelbom 
et.al. (1980). One then looses the notion of a true trajectory and its ap- 
proximation through a comparison dynamics. In particular, the selection of 
physical solutions can come only through an additional postulate. 

The current generated by a point charge can be written as 



f(x) = j dsv fl (s)5(x - z(s)) . (1.125) 

— oo 

McManus (1948) proposed to smear out the 5-function as 

f{x) = [ dsv»(s)f((x - z(s)) 2 ) , (1.126) 



which is to be inserted in the action ( |1.103| ), ( |1.104j ). The resulting equa- 



tions of motion are nonlocal in space-time. McManus did not identify the 
conserved four-momentum. The case of rectilinear motion only is discussed 
by Schwinger (1983). It seems to us that Nodvik (1964) is the only worked 
out example of an extended relativistic charge model. As to be explained, 
for slowly varying external potentials it agrees with experience over the full 
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range of allowed velocities. 

ad \1. T[ : The variation of Nodvik (1964), which includes the Euler angles of 
the body fixed frame, is adapted here to the restricted variation over the 
world line of the particle only. 
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2 Historical Notes 



2.1 Extended charge models (1897-1912) 

When in 1897 J.J. Thomson identified the cathode rays as consisting of parti- 
cles with charge — e, he had not only discovered the first elementary particle, 
but he also challenged the theoretical physicists to compute the energy- 
momentum relation of the electron. To put it slightly differently, we write 
the equations of motion in crossed E and B fields as 

m(v)v = e(E + c~ l v x B) (2.1) 

with m(v) the velocity dependent mass as a 3 x 3 matrix. The challenge was 
then to determine the ratio m(v) /e. In a long series is of experiments starting 
in 1897 the velocity dependent mass was measured by W. Kaufmann (1901). 
With improving technology the experiments were repeated and extended 
to a larger range of velocities. Latest by 1914 the relativistic dependence 
m(v) = mo7(l + 7 2 c~ 2 |v)(f|),7 = (1 — t> 2 /c 2 ) -1 / 2 , and m the rest mass of 
the electron, was regarded as well confirmed. 

So which theory could be used to determine m(v)7 In fact, there was 
little choice. Since the phenomenon under consideration is clearly electro- 
magnetic, one had to use Maxwell equations, and since the trajectory of a 
single charge was measured, one had to couple to the Lorentz force equation. 
Thus the electron was pictured as a tiny sphere charged with electricity. In 
the inhomogeneous Maxwell equations one had to insert the current gener- 
ated by that moving sphere. On the other hand the electromagnetic fields 
back react on the charge distribution through the Lorentz force equation. 
Thereby one has introduced what is called an extended charge model. Abra- 
ham (1905) adopted first a in the absolute laboratory frame rigid charge 
distribution. The corresponding energy-momentum relation is discussed at 
length in the second volume of his book on electromagnetism, compare with 



Section EO. For the Abraham model, Sommerfeld (1904,1905) obtained an 



exact equation of motion for the electron, which as a complicating and un- 
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familiar feature contained memory terms as a result of integrating over the 
retarded fields. Lorentz (1904a,b) advertised a charge distribution which is 
rigid in its momentary rest frame and therefore, as seen from the labora- 
tory frame, is contracted parallel to its momentary velocity Of course, it 
was left completely open by which forces this charge distribution is kept in 
place. Poincare (1906) developed non-electromagnetic models where addi- 
tional stresses counteracted the Coulomb repulsion. 

In all extended charge models the velocity dependent mass has the addi- 
tive structure m(v) = m^l + rrif(v), where is the bare mechanical mass 
of the particle, in accordance with Newtonian mechanics taken to be velocity 
independent, and where m^{v) is the mass due to the coupling to the field, 
which had to be computed from the model charge distribution. Clearly, in 
the experiment only the sum + rrif can be measured. Lorentz apparently 
favored to set = 0. His model yields then the usual relativistic velocity 
dependence, cf. Notes to Section |3.1| . However, in these times a fully rela- 
tivistic model was out of reach. Rather one relied on the semi-relativistic 
Abraham model of a rigid sphere and substituted at the appropriate places 
a relativistically contracting charge distribution. 

By 1904 the theoretical predictions were worked out with the experiments 
not yet precise enough to clearly distinguish between them. Nevertheless the 
whole enterprise came to a rather sudden end, since Einstein (1905a,b) force- 
fully argued that just like electromagnetism in vacuum also the mechanical 
laws had to be Lorentz invariant. But if Einstein was right, then the energy- 
momentum relation of the electron had to be the relativistic one, as empha- 
sized independently by Poincare (1906). Thus the only free parameter was 
the rest mass of the electron which anyhow could not be deduced from the- 
ory, since the actual charge distribution was not known. There was simply 
nothing left to compute. By latest in 1913 with the atomic model of Bohr, 
it was obvious that a theory based on classical electromagnetism could not 
account for the observed stability of atoms nor for the sharp spectral lines. 
As a tool to explain properties of atoms, electrons, and nuclei the classical 
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electron theory was abandoned. 

The effective equation of motion for the electron as given by Equation 
(2.1) could not possibly have been the full story. Through the work of Larmor 



it was already understood that a charge looses energy through radiation 
roughly proportional to v 2 . Lorentz observed that in the approximation of 
small velocities this loss could be written as the friction or radiation reaction 
force 

e 2 

F "- = ^— 3*, (2-2) 



which had to be added to the effective Lorentz force equation fl2.1| ). In 1904 
Abraham obtained this friction force for arbitrary velocities as 

e 2 

F rr = [ 7 4 C - 2 (v ■ v)v + 3 7 V 4 (v • vfv + 37 4 c" 2 (v • v)v + 7 2 -i)l , 

one 5 

(2.3) 

7 = (1 — v 2 1 c 2 ) -1 / 2 . He argued that energy and momentum is transported 
to infinity through the far field. On that scale the charge distribution is 
like a point charge and the electromagnetic fields can be computed from the 
Lienard-Wiechert potentials. Using conservation of energy and momentum 
for the total system he showed that the loss at infinity could be accounted for 
by the friction like force (|2.3| ). Von Laue (1909) recognized that the radiation 
reaction is relativistically invariant and can be written as, ' denoting now 
differentiation with respect to the eigentime, 

Frr = S~A^- C- 2 V X V X Vn . (2.4) 
OTIC 6 

This is how the radiation reaction appears in the famous 1921 review arti- 
cle of Pauli on relativity. But apparently, there was no incentive to study 
properties of the effective Lorentz force equation Q2.1|) including the full ra- 



diation reaction correction (|2.3| ). According to Schott (1912) after studying 
the motion in a uniform electric field: Hence the effect of the reaction due 
to radiation is quite inappreciable in this and probably in all practical cases. 
For applications, simpler phenomenological approaches sufficed. 
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The first chapter on the dynamics of classical electrons closes around 1912 
as compiled and worked out in great detail by Schott (1912). In essence there 
were two results: (i) a relativistically invariant expression for the radiation 
reaction and (ii) energy-momentum relations for the charged particle which 
were depending on the particular model charge distribution. However all 
models were inconsistent with Einstein's theory of special relativity. In par- 
ticular, the rest mass came out to be different from the electrostatic energy 
of the charge distribution. 

2.2 The point charge 

Our second chapter consists of a single paper: "Classical theory of radiat- 
ing electrons" submitted by P.A.M. Dirac on March 15, 1938. But before we 
have to follow up the intermission during which some research on the classical 
electron theory continued. We mention only the studies by Fermi (1922) and 
Frenkel (1925). Fermi argues that Abraham and Lorentz had not used the 
relativistically proper definition of energy and momentum which explained 
their disagreement with Einstein's theory. Frenkel, apparently influenced by 
Ehrenfest, proposes to consider electrons as undivisable, therefore without 
any extent, and studies the general structure of equations of motion. Of 
course, the most important event during the intermission was the develop- 
ment of quantum mechanics, which almost immediately after its discovery 
was applied to quantizing the electromagnetic field. Thereby the line shape 
and life time for excited states of atoms could be determined. Quantum 
mechanics gave a strong push to the classical theory. One had to quantize 
in Hamiltonian form. Thus the Lagrangian and Hamiltonian structure of 
the coupled Maxwell-Lorentz equations had to be explored, the role of the 
constraints and of the gauge freedom had to be understood. 

It became apparent fairly soon that the newly born quantum electrody- 
namics yields infinities when one tries to remove the ultraviolet cutoff, i.e. 
in the limit of a point charge distribution. Thus a problem which had been 
dropped over 15 years before reappeared in a different guise. In the '30 and 
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early '40 it was a fairly widespread believe that one way to overcome the diffi- 
culties of quantum electrodynamics is a better understanding of the classical 
theory of point charges coupled to their radiation field. Of course, this was 
only a vehicle to the final goal, namely a consistent quantized theory. We 
do not describe the various attempts, since the proper formulation of quan- 
tum electrodynamics eventually went a very different route. Dirac's paper 
was equally motivated by quantum electrodynamics. However, as such it is 
concerned only with the classical electron theory. 

We have to report the findings of Dirac in fair detail, since most further 
activities start from there. The formal argument in the original paper can 
be well followed and alternative versions can be found in Rohrlich (1990), 
Teitelbom et. al. (1980), Thirring (1997). So there is no need for repetition 
and we concentrate on the conclusions. At first reading it is mandatory to 
disregard all philosophical claims and to concentrate on the equations. But 
before, let us see how Dirac himself viewed the 1897-1912 period: The Lorentz 
model of the electron as a small sphere charged with electricity, possessing 
mass account of the energy of the electric field around it, has proved very 
valuable in accounting for the motion and radiation of electrons in a certain 
domain of problems, in which the electromagnetic field does not vary too 
rapidly and the accelerations of the electrons are not too great. Dirac wanted 
to construct quantum electrodynamics. There the electron is regarded as an 
elementary particle with, almost by definition, no internal structure. Thus 
Dirac had to dispense with model charges and to develop a theory of point 
electrons. 

What did Dirac really accomplish? Of course, he assumes the validity of 
the inhomogeneous Maxwell equations. The current is generated by a point 
charge whose motion is yet to be determined. Mechanically this point charge 
is relativistic with bare mass m b . There is no explicit back reaction of the 
field onto the charge, since at no stage Dirac would invoke the Lorentz force 
equation. Rather conservation of energy and momentum should suffice to fix 
the true trajectory of the point charge. Note that this is very different from 
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the extended charge models where the starting point is a closed systems of 
equations for the particle and the Maxwell field. Dirac studies the flow of 
energy and momentum through a thin tube of radius R around the world 
line of the particle. The computation simplifies by writing the retarded fields 
generated by the motion of the point charge as 

-Fret = r (-^ret + -F adv ) + - (F ret — -F adv ) (2.5) 

in all of space-time. The difference term turns out to be finite on the world 
line of the charge and yields in the limit R — > 0, through a balancing of 
energy and momentum, the relativistic radiation reaction fl2.4j ). Thus in 
retrospect one can understand why in the semi-relativistic Abraham model 
the radiation reaction is nevertheless of relativistic form. 

The more delicate term in ( |2.5| ) is the sum, which is divergent on the 
world line of the particle. At the expense of ignoring other divergent terms, 
cf. Thirring (1997), Equation (8.4.16), Dirac obtains the expected result, 
namely 

e 2 

-——v^ = -m^. (2.6) 
vkRc a 

Adding the radiation reaction ( |2.4| ) and equating with the mechanical four- 
momentum, the final result is an equation of motion which determines the 
trajectory of the particle, 

2 

(m b + m f K = m exp ^ = -F£ v v + [v» - c~ 2 v x v x v»] + O(R) (2.7) 

c bnc 6 

with an error of the size of the tube, where we have added in the prescribed 
tensor F£? of the external fields. 

To complete the argument Dirac has to take the limit R — ► 0. Since 
ffif — > oo, this amounts to 

mb — * — oo , m.f — > oo , m cxp = + rrif fixed , (2.8) 

where m exp is adjusted such that it agrees with the experimentally determined 
mass of the charged particle. ( |2.8|) is the classical charge renormalization. 
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Dirac admits that such a model is hardly a plausible one according to 
current physical ideas but this is not an objection to the theory provided we 
have a reasonable mathematical scheme. 

Equation (|2.7|), dropping the terms O(R), is the Lorentz-Dirac equa- 
tion. Within the framework of Dirac it makes no sense to ask whether the 
Lorentz-Dirac equation is "exact" , since there is nothing to compare to. The 
Lorentz-Dirac equation comes as one package, so to speak. One only could 
compare with real experiments, which is difficult since the radiation reaction 
is so small, or one could compare with higher level theories as quantum elec- 
trodynamics. But this has never been seriously attempted, since it would 
require to have a well defined relativistic quantum field theory which is a 
difficult task to begin with. 

The Lorentz-Dirac equation is identical to the effective equations of mo- 
tion obtained from extended charge models, if we ignore for a moment that 
the kinetic energy might come out differently depending on which model 
charge is used. In this sense Dirac has recovered the classical results through 
a novel approach. However there is an important distinction. For extended 
charge models one has a true solution for the position of the charged particle, 
say q{t). One can compare then q(t) with a solution of the Lorentz-Dirac 
equation and hope for agreement in asymptotic regimes, like slowly varying 
potentials. In addition for an extended charge model a negative bare mass 
might have drastic consequences which cannot be ignored. 

Dirac continues with an observation which shattered the naive trust in 
the classical electron theory. He observes that even for zero external fields 
(|2.7f ) has solutions where \v(t)/c\ — > 1 as t — > oo and \v(t)\ increases without 
bound. Such unphysical solutions he called runaway. If one inserts numbers, 
then runaways grow very fast. E.g. for an electron one has v(t) = v(0)e^ T 
with r = f 0~ 23 sec. If the Lorentz-Dirac equation (|2.7|) is a valid approxima- 
tion in an extended charge model, which after all was the main consensus of 
the f897-f9f2 period, then also there one encounters runaway solutions. It 
is somewhat surprising that apparently runaways went completely unnoticed 
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before, wich only proves that no attempt was made to apply the Lorentz- 
Dirac equation to a concrete physical problem. 

Dirac proposed to eliminate the runaway solutions by requiring the asymp- 
totic condition 

lim v^t) = . (2.9) 

t— >oo 

As additional bonus the problem of the missing initial condition is resolved: 



Since in (2.7) the third derivative appears, one has to know z M (0), i M (0), as 
in any mechanical problem, and in addition z M (0). If one accepts ( |2.9| ), the 
initial condition ?) M (0) is replaced by the asymptotic condition (|2.9| ). Dirac 
checked that for zero external forces and for a spatially constant but time- 
dependent force the asymptotic condition singles out physically meaningful 
solutions. 

By the end of 1938 the classical electron theory was in an awkward shape, 
in fact in a much worse shape than by the end of the 1912. Formal, but even 
by strict standards careful, derivations yielded an equation with unphysical 
solutions. How did they come into existence? While Dirac's asymptotic 
condition seemed to be physically sensible, it was very much ad hoc and 
imposed post festum to get rid of unwanted guests. Even those physicists 
willing to accept the asymptotic condition as a new principle, like Haag 
(1955), could not be too happy. Solutions satisfying the asymptotic condition 
are acausal in the sense that the charge starts moving even before any force 
is acting. To be sure the causality violation is on the time scale of r = 10~ 23 
sec for an electron, and even shorter for a proton, and thus has no observable 
consequences. But acausality remains as a dark spot in a relativistic theory. 

The clear recognition of runaway solutions generated a sort of consensus 
that the coupled Maxwell-Lorentz equations have internal difficulties. 

To quote from the preface of the book by Rohrlich: Most applications 
treat electrons as point particles. At the same time, there was the widespread 
belief that the theory of point particles is beset with various difficulties such as 
infinite electrostatic self-energy, a rather doubtful equation of motion which 
admits physically meaningless solutions, violation of causality, and others. It 
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is not surprising, therefore, that the very existence of a consistent classical 
theory of charge particles is often questioned. 

To quote from Chapter 28 of the Feynman Lectures: Classical mechanics 
is a mathematically consistent theory; it just doesn't agree with experience. 
It is interesting, though, that the classical theory of electromagnetism is an 
unsatisfactory theory all by itself. The electromagnetic theory predicts the 
existence of an electromagnetic mass, but it also falls on its face in doing so, 
because it does not produce a consistent theory. 

To quote from the textbook on mathematical physics by Thirring: Not 
all solutions to ( \2. Tj ) are crazy. Attempts, have been made to separate sense 
from nonsense by imposing special initial conditions. It is to be hoped that 
some day the real solution of the problem of the charge-field interaction will 
look differently, and the equations describing nature will not be so highly un- 
stable that the balancing act can only succeed by having the system correctly 
prepared ahead of time by a convenient coincidence. 

To be sure, these issues were of concern only to theoretical physicists 
in search for a secure foundation. Synchroton radiation sources were built 
anyhow. The loss in energy of an electron during one revolution can be ac- 
counted for by Larmor's formula. This is then the amount of energy which 
has to be supplied in order to maintain a stationary electron current. The 
radiation emitted from the synchroton source is computed from the inho- 
mogeneous Maxwell equations with a point charge source, i.e. from the 
Lienard-Wiechert potentials. No problem. 

2.3 Wheeler— Feynman electrodynamics 

To avoid the infinities of self-interaction Wheeler and Feynman (1945,1949) 
designed a radical solution, at least on the classical level since the quantized 
version of their theory was never accomplished. The Wheeler-Feynman the- 
ory departs in two essential aspects from standard electrodynamics, 
(i) The only dynamical variables are the trajectories of the charges. As such 
there are no electromagnetic fields, even though one uses them as a familiar 
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and convenient notational device. 

(ii) To achieve agreement with observation, the theory requires to have many 
particles. For example, in the two-body scattering problem there is no ra- 
diation damping. Such friction forces are understood as the result of the 
interaction with the charged particles in the surrounding matter. 

The starting point of the Wheeler-Feynman electrodynamics is an action 
first written down by Fokker (1929). Let us consider iV particles with mass 
rrii, charge e^, and motion given by the world line zrfifc), i = 1, . . . , N. The 
world line is parametrized by its eigentime and ' denotes differentiation 
with respect to that eigentime. The action functional has the form 

S = -^m,c 2 / yjiz^dn (2.10) 

i=l J 

+ e * e i / / 6 ( z d) - %))%) • Z{j)dTidTj . 



>,j=i 



A formal variation of S leads to the equations of motion 



N 

_ £i y> 1 

'(9 c 2 



Here F^-y F^ti) are ^ ne retarded and advanced Lienard-Wiechert fields 
generated by the charge at zu\ and evaluated at zuy They are derived from 
the retarded and advanced potentials 

A ret(j)( x ) = e i iJ 1 j) (r jrct )/(a; CT -2; (j)CT (r jrct ))i[ r i) (r irct ), (2.12) 

A adv(j)( X ) = e i^)( T iadv)/(^-2(i)a(r ja dv))i(;-)(r ja dv) (2.13) 

with rj ret , resp. r jadv , the eigentime when the trajectory zu^ crosses the 
backward, resp. the forward, light cone with apex at x. 

To transform ( |2.11| ) into a familiar form, we use the decomposition (|2.5|) 



and Dirac's observation that (F Tet — F a d v )/2 at the trajectory of the particle 
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yields the radiation reaction. Then 

N 



3=1 



N 1 

+7 EzKw^o) - • (2-14) 



Of course, being symmetric in time, we could have equally transformed to 
the advanced fields for the force and a radiation reaction with reversed sign. 
We note that in (|2J^) the mass of the particle is not renormalized. The 



retarded force is of the usual form. The radiation reaction has runaways. 
So one must either impose the asymptotic condition ( |2.9| ) or have the good 
faith that (|2.11 ) does not posess such unphysical solutions. The last term in 



( 2.14| ) is unwanted and Wheeler and Feynman spend a considerable effort to 



argue that for a sum over a large number of charges in disordered motion this 
last term vanishes. If it is exactly zero, the condition of a perfect absorber is 
satisfied and the standard equations of motion for charged particles result. 

As with Dirac, one can accept only the whole Wheeler-Feynman package. 
Consequently, there has been little further work on the theory. In particular, 
it has never been checked how well the assumption of a perfect absorber is 
satisfied. 



Notes and References 

ad [|: A more detailed account on the history of the classical electron theory 
can be found in Pais (1972,1982), Rohrlich (1973), and in the introductory 
chapters of Rohrlich (1990). The interconnection with quantum electrody- 
namics before the 1947 Shelter Island conference is well described in Schweber 
(1994). 



ad \2.2j . Kramers (1948) investigations on the mass renormalization in the 



classical theory were instrumental for a correct computation of the Lamb 
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shift. We refer to Dresden (1987 ) and Schweber (1994). 

ad \2.3[ : The two-body problem in Wheeler- Feynman electrodynamics is dis- 
cussed by Schild (1963). The existence and classification of solutions is stud- 
ied by Bauer (1997). A few explicit solutions are listed in Stephas (1992). 
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3 Energy— Momentum Relation 



For the Abraham model we established already that its energy £ is conserved, 
S given by Equation (|1.42| ). If the external fields vanish, then the dynamics 
is invariant under spatial translations. Thus the total momentum, denoted 
by V, must also be conserved. The minimum of S at fixed V defines the 
energy-momentum relation. 

If the external forces vanish, the simplest solution to the equations of 
motion has the particle travelling at constant velocity v in company with its 
electromagnetic fields. There seems to be no accepted terminology for this 
object. Since it will be used as a basic building block later on, we need a 
short descriptive name and we call this particular solution a charge soliton, 
or simply soliton, at velocity v, in analogy to solitons of nonlinear wave 
equations. The soliton has an energy and a momentum which are linked 
through the energy-momentum relation. 

In the following two sections we compute the conserved energy and mo- 
mentum, the charge solitons, and the energy-momentum relation for both 
the Abraham and the Lorentz model. We will assume ex = 0, A ex = 
throughout. 

3.1 Abraham model 

The mechanical momentum of the particle is given by 



rrib'yv 



(3.1) 



and the momentum of the field by 




(3.2) 



Thus we set the total momentum 



V = m h ^v + Vf 



(3.3) 
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as a functional on M. . One easily checks that V is conserved by the Maxwell- 
Lorentz equations ( |1.37| ) - ( |1.39|) . The corresponding Lagrangian, compare 
with (|1.41|) , is invariant under spatial translations and V is the conserved 
quantity which, by Nother's theorem, corresponds to this symmetry. 

We want to minimize the energy at fixed total momentum. We eliminate 
v between ( p..42|) and (|3.3j ) and thus have to minimize 



K + P-J d '< E x B )f) 1/2 + \J d "< E2 + B " 



(3.4) 



at fixed V and subject to the constraints V ■ E = p ,V ■ B = 0. By translation 
invariance we can center p at an arbitrary q 6 R 3 . For q = 0, say, the 
minimizer is unique and given by 

E v (x) = -V(j) v (x)+v(v-V(j) v (x)), 

B v (x) = -v x V(j) v {x) , (3.5) 

where 

f v {k) = [k 2 - (v ■ k) 2 }- 1 ^) (3.6) 

and in physical space 

(j) v {x) = J d 3 y(47T^-*(x - yf + (v-(x- y)f)~ 1 p{v) . (3.7) 

Here \v\ < 1, i.e. v e V, and v has to be adjusted such that V = P s (v) with 

P s {v) = m h -fv + J d 3 k\p(k)\ 2 ([k 2 - (k- v) 2 ]-^ 

- 7 - 2 [fc 2 - (k ■ v) 2 ]- 2 (k ■ v)k) 
-1 + V 2 , 1 + \v\ 



V 



{m b 7 + m e \v\ 2 [ log- — - 1 1 } , (3.i 

1 L 2m 1 — \v\ 



where m e is the electrostatic energy of the charge distribution p, 

m e = I J d 3 xd 3 x'p(x) p(x')(4:7i\x - x'^ 1 . (3.9) 
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The map Vbuh P s (v) G M 3 is one to one and therefore V = P s { v ) has a 
unique solution. The minimizing energy is given by 

E s {v) = m b7 + i J(Pk\p(k)\ 2 [k 2 -(k-v) 2 ]- 2 ((l + v 2 )k 2 -(3-v 2 )(v-k) 2 ) 

r 1 1 + \V\ n 

= m b 7 + m e — log- — r - 1 . (3.10) 

L \v\ 1 — \v\ J 

Eliminating now v between E s and P s yields the energy-momentum relation 

E eS (p) = E s {v(p)) , (3.11) 

where v(P s ) is the function inverse to P s (v). As to be underlined, E e ^ 
depends on the charge distribution only through its electrostatic energy. 
We note that 

P s {v) = V v T(v) (3.12) 

with 

T{v) = -m hl ~ l + \r 2 j d'klpik^lk'-ik-v) 2 ]- 1 (3.13) 
and that 

E s (v) = P s (v)-v-T(v). (3.14) 

This suggests that T will play the role of the inertial term in an effective 
Lagrangian and E s the role of an effective Hamiltonian as our notation in 
(|3.11| ) indicates already. In particular, 



v = V p E cS (p) (3.15) 

and, equivalent ly, 

v ^l = VvEs (v) (3.16) 
which implies that v is to be interpreted as velocity. 
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For a relativistic theory one expects that 



E s {v) = (m b + m e )7, P s (v) = (m b + m e )jv . (3.17) 

Since the Abraham model is semi-relativistic, there is no reason for such 
a property to be satisfied. Still we found that, as in the relativistic case, 
the energy-momentum relation depends on the charge distribution p only 
through m e . 

For small v we have 

E s (v) -E a (0) = ^{m h + -m e )v 2 , P s (v) = (m h + ^m e )v . (3.18) 

Thus the effective mass in the nonrelativistic approximation is 

4 

m cS = m h + -m e . (3.19) 
o 

In Figure 2 we plot E s (v), P s (v) for the extreme case m b = and compare 
with the relativistic dispersion of mass |m e . Clearly at speeds \v\ > 0.3, the 
Abraham model looses its empirical validity. One could partially save the 
Abraham model by declaring the Compton wave length as the characteristic 
size of the charge distribution. Then m e /m b = 0.01 and the relativistic 
dispersion is violated only for speeds close to one. 

The energy minimizer has a simple dynamical interpretation. We look for 
a solution, S q>v , of ( |1.37| ) , ( |1.39| ) travelling at constant velocity v and find 



S q>v {t) = (E v (x -q-vt), B v (x - q-vt),q + vt, v) (3.20) 

with v G V, q G M 3 , and E V ,B V from ( |3.5|) . S q>v is the charge soliton 
labeled by its center q and velocity v. It has the energy S(S qtV ) = E s {y) and 
momentum V(S q ^ v ) = P s (v). 

There is an instructive alternate way to represent the charged soliton. 
We consider the inhomogeneous Maxwell equations ( |1.37| ) and prescribe the 
initial data at time r. We require that the particle travels along the straight 
line q = vt. If we let r — > — oo and consider the solution at time t — 0, then 
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in ( |1.16 ), ( 1.17 ) the initial fields will have escaped to infinity and only the 



retarded fields survive. Using ( |1.16|) , (|1.17|) this leads to 

o 

E v (x) = - J dt J d 3 y (VG_t (x - y) p(y - vt) 

— oo 

+d t G- t (x-y)vp(y-vt)), (3.21) 

o 

: J dt J d 3 yVx G- t (x - y) vp{y - vt) , (3.22) 



BJx) 



which can be checked either directly in Fourier space or as being a solution 
of the Maxwell equations travelling at constant velocity v. 

3.2 Lorentz model 

We look for a solution travelling at constant velocity v. Since the model is 
relativistic, we first determine the four-potential of the charge soliton in its 
rest frame which yields 

AV) = (9(x' 2 ), 0) (3.23) 
with g(|a?| 2 ) = — (A~ 1 p)(x). Then in the laboratory frame S 

A»(x) = v»g(x 2 - (v x x x ) 2 ) (3.24) 
and the electromagnetic field tensor has the form 

F" v (x) = d»A v tx)-d v A»{x) (3.25) 
= 2 [{a? - {v x x x )v»)v v - {x v - Kx A KK] g'{x 2 - (v x x x ) 2 ) . 



F^ v indeed satisfies the Maxwell equations (|1.105|) with the current 

f(x)=v»ef(x 2 -(v x x x ) 2 ), (3.26) 
where e/(|a;| 2 ) = p(x), in accordance with (|1.36|) , (|1.85|) . 

63 



Expressed in terms of electric and magnetic fields we have 

E(x,t) = —V<fi v (x — q — vt)+v(v-V<fi v (x — q — vt)), 

B(x,t) = -v x <p v (x - q - vt) , (3.27) 

where 

$ v (k) = [k 2 - {v • fc) 2 ]- 1 p( 7 - 1 fc,| + k ± ) (3.28) 

with fey parallel and fej_ orthogonal to v. In contrast to the nonrelativistic 
coupling, the charge distribution is now Lorentz contracted as seen from the 
laboratory frame, compare with ( p.6|) , where we note that k 2 — {y ■ k) 2 = 
7 - 2 fcJ + k\ and p{T l h + k ± ) = p r ([k 2 - (v ■ k) 2 } 1 ' 2 ). 

To determine energy and momentum of the relativistic soliton we first 
have to find out how these quantities are even defined. We start from the 
energy-momentum tensor of the electromagnetic field 

T^ u = F uX F x " - g" u F a/3 F al3 . (3.29) 

^From the Maxwell equations it satisfies the local balance 

d v T^ = F^j u . (3.30) 

We now claim that 

^(m b ^ + J d A xT^{x)v u 5{v x {x - z) x )) =0. (3.31) 

Thus it is natural to regard 

pi" = m h v^ + J Ar(i)^«5(?; A (j; - z) x ) (3.32) 

as the conserved four-momentum. 

To derive (|3.31| ) we multiply (|3.30|) by [1 — v a (x — z) a ] 8(v\(x — z) x ) and 



integrate over all space-time. Then 

d A xF^{x) 3u {x) [1 - v a {x - zf] 6(v x (x - z) x ) (3.33) 

d*xd v T^{x) [1 - v a {x - zY] 6{v x {x - z) x ) 

- J d A x T^(x) [ - v u 5(v x (x - z) x ) + (1 - v a (x - zY)v v 5'(v x (x - z) x )} . 
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We have 
d 



-f / d^xT^ix^div^x- z) x ) 
ds 



d 4 x T^(x) (v v 5{v x {x - z) x ) - v v {l - v a (x - zY) 8'{v x {x - ^) A )) 
d 4 x F^{x) ]u {x){l - v a (x - zf) 8{v x {x - z) x ) 
d A xF» v {x){l - v a {x - zf) 5{v x {x - z) x ) 

x J ds' [(«„ - v v v a {x - z)°)f({x - zf) 8{v x (x - zf)](s') . (3.34) 

Under our assumption of not too large an acceleration, compare with ( |1.90| ) 
below, the hyperplane {x : V\(s)(x — z(s)) x = 0} intersects the hyperball 
{x : vx(s') (x - z(s')) x = 0, f((x - z(s')) 2 ) > 0} only if s = s' and 

(1 - v a {s){x - z(s)f) 8(v x (s)(x - z(s)f) 8(v x (s')(x - z(s')) x ) 

= 5{v x {s){x - z{s)) x )5{s - s') . (3.35) 

Thus 



(3.34) 



d A x F^{x) 5{v x {x - z) x ) (v u - r) ua {x - zY) f{{x - zf) 



(3.36) 



where we used the equations of motion (|1.106|) . This proves (|3.31 



The expression ( |3.32j ) for P^ is covariant. Thus we are allowed to work 
out the integral in the frame S' traveling with velocity v relative to S. In 
this frame 

p'° = 7 ( mb + I Jd 3 x'E(x'f), P'» = for fx = 1,2,3 (3.37) 



and thus 



pv = (m h + m e ) 7 u A 



(3.38) 



which shows that the Lorentz model has the physically correct relativistic 
four-momentum. Of course, experimentally only the sum, + m e , can be 
observed. 
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Notes and References 



ad WA , \3. %\ : Abraham (1905) computed the energy-momentum relation in 
essence along the same route as outlined here (except for the variational 
characterization). Sommerfeld (1904,1905) used the expansion of the exact 
self-force, as will be explained in Chapter |5|. Lorentz (1904a) proposed a 
model charge which relativistically contracts along its momentary velocity. 
Thus provisionally we replace the charge distribution p{x) by its Lorentz 
contracted version 

p L (x) = 7Pr([* 2 + 7 2 (*^) 2 ] 1/2 ), (3-39) 
p L (k) = p r ([k 2 - (v • k) 2 ] 1 ' 2 ) • (3.40) 



This expression is substituted in ( |3.6[ ) and gives the electromagnetic fields 
comoving with the charge at velocity v. Their energy and momentum is 
computed as before with the result 

4 

P L (v) = v(m h ^(v) + -m e j(v)) , (3.41) 

E L (v) = m hl (v)+m cl (v)(l + ^v 2 ) . (3.42) 

The momentum has the anticipated form, except for the factor 4/3 which 
should be 1. The energy has an unwanted v 2 /3. In particular the relation 
( 3.16p does not hold, which implies that the power equation ^El(v) differs 
from the force equation v ■ |Pi(u). We refer to Yaghjian (1992) for a 
thorough discussion. 

Schott (1912) employed as a model charge a deformable elastic medium. 
To compute the velocity dependent mass he used in principle the same 
method as Sommerfeld, an exact self-force and an expansion in the charge 
diameter. Schott considered also electron models different from those of 
Abraham and Lorentz. 

There have been various attempts to improve on the oversimplistic version 
( |3.39| ) of the Lorentz model. Fermi (1922) argues that in a relativistic theory 
energy and momentum have to be redefined. His argument has been redis- 
covered several times and is explained in Rohrlich (1990). Poincare (1906) 
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takes the elastic stresses into account. His theory is excellently presented 
in Yaghjian (1992). The material of Section |3.2j is adapted from of Nodvik 
(1964). 
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4 Adiabatic Limit 



If we assume that the mass of an electron is purely electromagnetic, then 
by equating its rest energy and electrostatic Coulomb energy the charge 
distribution p must be concentrated in a ball of radius 

e 2 

R p = — - = 3x 1(T 13 cm (4.1) 

mc 2 

which is the so called classical electron radius. Quantum mechanically one 
argues that through fluctuations the electron appears to have an effective 
size of the order of the Compton wave length A c = hm/c = (e 2 /hc)~ l R p = 
137 R p . Electromagnetic fields which can be manipulated in the laboratory 
vary little over that length scale. R p defines a time scale through the time it 
takes light to cross the diameter of the charge distribution, 

t p = Rp/c = 1CT 23 sec , equivalently as a frequency, uj p = 10 23 Hz . (4.2) 

Again, manufactured frequencies are much smaller than u p . Space—time vari- 
ations as fast as (|4.1|) and ( |4.2| ) lead us deeply into the quantum regime. Thus 
it is a natural and physically a mandatory problem to study the dynamics of 
a charged particle under external potentials which vary slowly on the scale of 
R p . This means we have to introduce a scale of potentials an enquire about 
an approximately autonomous particle dynamics with an error depending on 
the scale under consideration. We will introduce such a scheme formally in 
the following section. The resulting problem has many similarities with the 
derivation of hydrodynamics from Newtonian particle dynamics - with the 
most welcome addition that it is simpler mathematically by many order of 
magnitudes. Still, the comparison is instructive. 

4.1 Scaling limit 

We assumed that in the Lorentz force equation there are in addition to dy- 
namical fields E(x,t), B(x,t) also prescribed external fields acting on the 
particle, which are the gradients of the external potentials <j) ex (x), A ex (x), 
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compare with Equation ( |1.39| ). We want to impose that ex and A ex are 
slowly varying on the scale of R p . Formally we introduce a small dimension- 
less parameter e and consider the potentials 

4> cx (ex) , A cx (ex) , (4.3) 

which are slowly varying in the limit e —>■ 0. Most of our results extend 
to potentials which vary also slowly in time. But for simplicity we restrict 
ourselves to time-independent potentials here. Clearly, e appears as a pa- 
rameter of the potential, just like ujq is a parameter of the harmonic potential 
| mujQX 2 . But one really should think of e as a book keeping device which or- 
ders the magnitude of the various terms and the space-time scales in powers 
of e. Such a scheme is familiar in very diverse contexts and appears whenever 
one has to deal with a problem involving scale separation. 

So how small is e ? From the discussion above one might infer that 
if (f) ex , A cx vary over a scale of 1 mm, then e = 1CT 12 . This is a strictly 
meaningless statement, because e<f) ex , - c A ex have the dimension of an energy 
and thus the variation depends on the adopted energy scale. In ( fOD we fix 
the energy scale and merely stretch the spatial axes by a factor e^ 1 . Since 
from experience this point is likely to be confusing, let us consider the specific 
example of a charge circling in the uniform magnetic field (0,0, B ). Since 
the corresponding vector potential is linear in x, to introduce e as in fl4.3| ) 
just means that the magnetic field strength equals eBq and the limit e — > 
is a limit of small magnetic field relative to some reference field B . Thus to 
obtain e we first have to determine the reference field and compare it with 
the magnetic field of interest. This shows that in order to fix e we have to 
specify the physical situation concretely, in particular the external potentials, 
the mass and charge of the particle, 7(v), and the time span of interest. 

The scaling scheme ( |4.3j ) has the enormous advantage that the analysis 
can be carried out in generality. In a second step one has to figure out e 
for a concrete situation, which leads to a quantitative estimate on the error 
terms. E.g. if in the case above we consider an electron with velocities such 
that 7 < 10, then, by comparing the Hamiltonian term and the friction term, 
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the reference field turns out to be B = 10 17 Gauss. Laboratory magnetic 
fields are less than 10 5 Gauss and thus e < 10~ 12 . In practice, e is always 
very small, less than 10~ 10 . This means that, firstly, all corrections beyond 
radiation reaction are negligible. Secondly, we do not have to go each time 
through the scheme indicated above and may as well set e — 1 thereby 
returning to the conventional units. Still on an theoretical level the use of 
the scale parameter e is very convenient. In Section |4.3| we will work out 
the example of a constant magnetic field more explicitly. If the reader feels 
uneasy about the scaling limit, (s)he should consult this example first. 
Adopting (|4.3|), the Lorentz force equation reads now 

— (m h jv(t)) = / d 3 x p{x - q(t)) [eE ex (ex) + E(x, t) 

+v(t) x (sB ex (ex) + B(x, t))] , (4.4) 

where 

E ex = - V0 ex , Bex = Vx A cx . (4.5) 



It has to be supplemented with the Maxwell equations ( |1.37|) , (|1.38|) . Our 



goal is to understand the structure of the solution for small e and as a first 
qualitative step one should discuss the rough order of magnitudes in powers 
of e. But before we have to specify the initial data. We give ourselves q°, v° 
as initial position and velocity of the charge. The initial fields are assumed 
to be Coulombic, centered at q° with velocity v°, i.e. 

Condition (I): 

Y(0) = S q o, v o(0), (4.6) 

compare with ( |3.20| ). Equivalently, according to ( |3.21| ), ( |3.22| ), we may 
say that the particle has travelled freely with velocity v for the infinite 
time span (— oo,0]. At time t = the external potentials are turned on. 
More geometrically we define the six-dimensional charge soliton manifold 
S = {S qjV , q G IR 3 , v G V} as a submanifold of the phase space M.. Then 
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our initial data are exactly on S. If there are no external forces, the solution 
remains on S and moves along a straight line. For slowly varying external 



potentials as in (|4.3|) we will show that the solution remains e-close to S in 
the local energy distance. 

On general grounds one may wonder whether such specific initial data 
are really required. In analogy to hydrodynamics, we call this the initial 
slip problem. In times of order t p , the fields close to the charge acquire their 
Coulombic form. However, during that period the particle might gain or loose 
in momentum and energy and the data at time t p close to the particle are 
approximately of the form S q ^, where q and v are to be computed from the 
full solution. Of course, at a distance ct away from the charge, the field still 
remembers its t — data. Thus we see that the initial slip problem translates 
into the long time asymptotics of a charge at zero external potentials but with 



general initial field data. We study this point in more detail in Section [4~3 . 
At the moment we just circumvent the initial slip by fiat. 

Let us discuss the three relevant time scales, where we recall that t p = 

R P /c 

(i) Microscopic scale, t = 0(t p ), q = 0(R P ). On that scale the particle 
moves along an essentially straight line. The electromagnetic fields adjust 
themselves to their comoving Coulombic form. As we will see, they do this 
with a precision 0(e) in the energy norm. 

(ii) Macroscopic scale, t = 0(e~H p ), q = 0(e~ l R p ). This scale is de- 
fined by the variation of the potentials, i.e. on that scale the potentials are 
4> ex (x), A cx (x). The particle follows the external forces. Since it is in com- 
pany with the almost Coulombic fields, the particle responds to the forces 
according to the effective energy-momentum relation, which we determined 
in the previous section. On the macroscopic scale the motion is Hamiltonian 
up to errors of order e. There is no dissipation of energy and momentum. 

(iii) Friction scale. Accelerated charges loose energy through radiation, which 
means that there must be friction corrections to the effective Hamiltonian 
motion. According to Larmor's formula the radiation losses are proportional 
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to v(t) 2 . Since the external forces are of the order e, these losses are pro- 
portional to e 2 when measured in microscopic units. Integrated over a time 
span e~H p the friction results in an effect of order e. Thus we expect order e 
dissipative corrections to the conservative motion on the macroscopic scale. 
Followed over the even longer time scale e~ 2 t p , the radiation reaction results 
in 0(1) deviations from the Hamiltonian trajectory. 

On the friction time scale the motion either comes to a stand still or stays 
uniform. In addition, as to be shown, the dissipative effective equation has 
the same long time behavior as the true solution. Thus we expect no further 
qualitatively distinct time scale beyond the friction scale. 

/,From our description, in a certain sense, the most natural scale is the 
macroscopic scale and we transform the Maxwell-Lorentz equations to this 
new scale by setting 

t' = et , x = ex. (4.7) 

We have the freedom of how to scale the amplitudes of the dynamic part of 
the electromagnetic fields. We require that their energy is independent of e. 
Then 

E'(x', t') = e~ 3/2 E(x, t) , B'(x', t') = e~ 3/2 B(x, t) . (4.8) 

Finally the new position and velocity are 

q'(t') = eq(t) , v\t') = v(t) , (4.9) 

so that Jp- q' = v'. There is little risk in omitting the prime. We denote then 

q £ {t) = eq{e-H) , v £ {t) = v{e~H) , p e (x) = e~ 3 p^x) , (4.10) 

which means that J d 3 x p £ (x) = e independent of e and p £ is supported in 
a ball of radius eR p . In the macroscopic coordinates the Maxwell-Lorentz 
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equations read 

d t B(x,t) = -VxE(x,t), 

d t E(x,t) = V x B(x,t) -V^Pe(x-q £ (t))v £ (t), (4.11) 
j {m hl v £ (t)) = E cx * p £ (q £ (t)) + v £ (t) x B cx * p £ (q £ (t)) 

+ J d z x^fsp e {x-q £ {t))[E(x,t) + v £ (t) x B(x,t)} 

together with the constraints 

V ■E = y/ep(--tf{t)) , V-B = 0. (4.12) 
On the macroscopic scale the conserved energy is 

£mac = m hl {v) + CX * p £ {q) + ^Jd 3 x (E(x) 2 + B{x) 2 ) . (4.13) 

Also the initial data have to be transformed and become now 
Condition (I e ): 

Y £ (0) = S £ qQjV o = (E v o(x - q°), B v o(* ~ <7°), 9°, «°) (4.14) 

with 

E v = -V4>% + v(v • V0 £ J , B v = -v x V<j>% , (4.15) 

where now 

#(fc) = Vip(ek)/[k 2 -(vk) 2 }. (4.16) 



On the macroscopic scale, the scaling parameter e can be absorbed into 
the "effective" charge distribution \fep £ . Its electrostatic energy, 



\ J d 3 ke\p £ (k)\ 2 ^ = \ J d 3 k\p(k)\ 2 p , (4.17) 
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is independent of e, whereas its charge 

J d 3 xy/ep £ (x) = yfee (4.18) 

vanishes as y/e. Recall that e is a book keeping device. 

We argued that on the macroscopic scale the response to the external 
potentials in the motion of the charges is of order one. We thus expect that 
q e (t) tends to a nondegenerate limit as e —* 0, i.e. 

lim q £ (t) = r(t) , lim v £ (t) = u(t) . (4.19) 

The position r(t) and velocity u(t) on the macroscopic scale should be gov- 



erned by an effective Lagrangian. In Section 3.1 we determined already the 



effective inertial term. If the potentials add in as usual, we have 

L cS (q, q) = T(q) - e(0 cx (q) - q ■ A cx (q)) , (4.20) 

which results in the equations of motion 

r = u , m{u)u = e(E ex (r) + u x B cx {r)) . (4-21) 

The velocity dependent mass m(u) has a bare and a field contribution. ^From 
( |3.12j ) we conclude that 

m{v) = ^1 (4.22) 
dv 

as a 3 x 3 matrix. If instead of the velocity we introduce the canonical 
momentum, p, then the effective Hamiltonian reads 

H cS (r,p) = E cS (p - eA cx (r)) + e0 ex (r) (4.23) 

with Hamilton's equations of motion 

r = V p H cS , p = -V r H cS . (4.24) 



Our plan is to establish the limit ( f4.19|) and to investigate the corrections 
due to radiation losses. 
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4.2 Comparison with the hydrodynamic limit 



In hydrodynamics one assumes that a small droplet of fluid with center r 
has its intrinsic velocity, u(r), and that relative to the moving frame the 
particles are distributed according to thermal equilibrium with density p(r) 
and temperature T(r). For such notions to be reasonably well defined, the 
hydrodynamic fields p, u, T must be slowly varying on the scale of the typ- 
ical interparticle distance. This is how the analogy to the Maxwell-Lorentz 
equations arises. As for them we have three characteristic space-time scales. 

(i) Microscopic scale. The microscopic scale is measured in units of a col- 
lision time, resp. interatomic distance. On that scale the hydrodynamics 
fields are frozen. Possible deviations from local equilibrium relax through 
collisions. To prove such a behavior one has to establish a sufficiently fast 
relaxation to equilibrium. For Newtonian particles no method is available. 
For the Maxwell field the situation is much simpler. Local deviations from 
the Coulomb field are transported off to infinity and are no longer seen. 

(ii) Macroscopic scale. The macroscopic space-time scale is defined by the 
variation of the hydrodynamic fields. If, as before, we introduce the dimen- 
sionless scaling parameter e, then space-time is 0(e~ r ) in microscopic units. 
On the macroscopic scale the time between collisions is 0(e), the interpar- 
ticle distance 0(e), and the pair potential between the particles at positions 
q { , qj is V(e~ 1 (q i — qj))- On the macroscopic scale the hydrodynamic fields 
evolve according to the Euler equations. These are first order equations, 
which must be so, since space and time are scaled in the same way. The 
Euler equations are formally of Hamiltonian form. There is no dissipation, 
no entropy is produced. In fact, there is a slight complication here. Even for 
smooth initial data the Euler equations develop shock discontinuities. There 
the assumption of slow variation fails and shocks are a source of entropy. 

(iii) Friction scale. In a real fluid there are frictional forces which are re- 
sponsible for the relaxation to global equilibrium. One adds to the Euler 
equations diffusive like, second order in spatial derivatives, terms and obtains 
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the compressible Navier-Stokes equations incorporating the shear and vol- 
ume viscosity for friction in momentum transport and thermal conductivity 
for friction in energy transport. On the macroscopic scale these corrections 
are of order e. In the same spirit, based on the full Maxwell-Lorentz equa- 
tions, there will be dissipative terms of order e emerging which have to be 
added to ( 4.2ip . Of course, in this context one only has to deal with ordinary 
differential equations as effective dynamics. 



4.3 Initial slip 

We adjusted the initial data for the electromagnetic field to be exactly on 
the charge soliton manifold, which physically means that without external 
forces the particle would travel forever at constant velocity accompanied 
by its comoving Coulombic fields. One may wonder whether such a rigid 
assumption is really necessary. Let us consider then times on the microscopic 
scale of the order e~ 1+s with some small 5 > 0. The external forces, which 
are 0(e), are still negligible, but the microscopic time span diverges. Thus 
we are led to investigate the long time limit of the Abraham model for zero 
external fields. Roughly one has the following picture: initially there is some 
exchange of momentum and energy between particle and field, but, since 
the total energy is bounded, eventually the particle relaxes to some definite 
velocity and the field builds up its comoving Coulombic shape. Thus after a 
very short macroscopic time the state is already close to the charge soliton 
manifold from whereon the adiabatic dynamics applies. 

To prove such a behavior we need little bit of preparation. Firstly we 
must have some decay and smoothness of the initial fields at infinity. We 
introduced already such a set of "good" initial data, M° ', in Section [L5] and 
therefore require Y(0) G Ai' 7 , < a < 1. Secondly, we need a notion of 
closeness of the fields. At a given time and far away from the particle the 
fields are determined by their initial data. Only close to the particle they are 
Coulombic. Therefore it is natural to measure closeness in the local energy 
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norm defined by 

\\(E,B)\\l=± J d 3 x(E(x) 2 + B{xf) (4.25) 

{\ X \<R} 

for given radius R. 

The true solution is Y(t) = (E(x,t), B(x,t), q(t),v(t)) which is to be 
compared with the charge soliton approximation (E v m(x — q(t)), B v (a(x — 
q(t)),q(t),v(t)),d. (pD. WesetZi(a:,t) = E(x,t)-E v{t) (x-q(t)), Z 2 (x,t) 
B(x,t) — B v ( t ){x — q(t)), Z = (Z\,Z<i) and want to establish that, for fixed 
R, \\Z(- + q(t), t)\\n— > for large times t — > oo. 

Proposition 4.1 For the Abraham model with zero external potentials and 



satisfying (C) let \e\ < e according to Theorem \5. i| and let the initial data 



Y(0) G Ai a for some a G (0, 1]. Then for every R> we have 

\\Z{. + q{t) 1 t)\\R<C R {l + \t\)- l - (T . (4.26) 
In particular, the acceleration is bounded as 

\v{t)\ < C(l+ \t\)- l - a (4.27) 
and there exists a G V such that 

lim v(t) = Voo • (4.28) 
Proof. From the Lorentz force equation and since \v(t)\ <v< 1 we have 

\v(t)\<Ce\\Z(- + q(t),t)\\ Rp . (4.29) 



Therefore (fOT|) follows from ( ^26|) . Then v(t) = v(0) + J dsv(s) and 



o 



\v(t)- Voo \ <C{l + \t\)-°. 

It remains to establish (|4.2(j|), which uses the method described in Ap- 



pendix Since Z(0) ^ 0, ( |5.35| ) reads now 



t 



Z(t) = e M Z(0) - J ds e A{t ~ s) g(s) . (4.30) 
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For the integrand in the second term we have the bound, compare with fl5.36| ) , 

\\e A(t ~ s) g(s)\\ Rp < C(v)e 2 (l + (t - s) 2 y l \\Z{. + q(s), s)\\ Rp . (4.31) 

For the first term of ( ^30|) we note that Z x (x, 0) = E°(x) - E v o(x) G M a 
by assumption. Using the solution of the inhomogeneous Maxwell equations 
in position space and the bound ( |1.47[ ) we have 

\Z 1 (x,t)\ + \Z 2 (x,t)\ < CT 2 J d i y5{\x-y\-t){\Z l {yM + \Z2{yM) 

+CT 1 J d 3 y5(\x-y\-t)(\VZ 1 (y,0)\ + \VZ 2 (y,0)\) 

< Cr 2 J d 3 yS(\x-y\-t)(l+ lyl)- 1 - 
+Cr 1 J d 3 yS(\x -y\- t)(l + \y\)- 2 -° . 

< C(l + ty l - u . (4.32) 
We choose R > R p . Then from flOID and ( |Og ) 



W + 9(t),t)IU < C(l + t)- 1 -' + C(7J)e 2 Jds(l + (t- s) 2 )- 1 ||Z(- + q(s), s) 

o 

(4.33) 

Let « = sup(l + t) 1+CT ||Z(- Then 

k < C + C(*;)e 2 (J ds (1 + (t - s) 2 )" 1 ^ + s)- 1 ^)/? , (4.34) 
o 

which implies k < oo provided C(?J) e 2 is sufficiently small. □ 
4.4 Appendix: How small is e ? 

We consider an electron moving in an external magnetic field oriented along 
the z-axis, B cx — (0, 0, Bq). The corresponding vector potential is A ex [x) — 
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~ Bq(— X2, xi, 0). According to our convention the slowly varying vector po- 
tential is given by A cx (ex) = ^eB (—X2,Xi,0). Thus B is a reference field 
strength, which we will determine, and B = sBq is the physical field strength 
in the laboratory. The motion of the electron is assumed to be in the 1-2 
plane and we set v = (it, 0). According to Section |7.2| , Example (ii), within 
a good approximation the motion of the electron is governed by 

711 = uj^u 1 - — (3uj c u) . (4.35) 

Here u 1 - = (— u 2 ,Ui), u c = eB/m c is the cyclotron frequency, and (3 = 
e 2 /67rc 3 mo. The first term is the Lorentz force and the second term accounts 
for the radiation reaction. 

We choose now the reference field B Q such that both terms balance, i.e. 

B = (Pe/moc)- 1 . (4.36) 

For electrons 

B = 1.1 x 10 17 Gauss (4.37) 

and even larger by a factor (1836) 2 for protons. For a laboratory field of 10 5 
Gauss this yields 

£ = 10- 12 . (4.38) 
Written in units of B (|4.35|) becomes 



7 it = elu^u 1 - epu° c u) (4.39) 

with f3u® = 1, lj® = eB^/m^c = 1.6 x 10 28 /sec. Thus friction is of relative 
order e and higher order corrections are then of relative order e 2 . As to be 
demonstrated, the dimensionless scaling parameter e merely serves as a book 
keeping device to keep track of the relative order of the various contributions. 
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4.5 Appendix: Point charge limit, negative bare mass 

The convential point charge limit is to let the diameter of the total charge 
distribution R p — > such that the charge remains fixed. Physically, this 
means that the charge diameter is small in units of the variation of the 
external potential, since this is the only other length scale available. At first 
sight, one just seems to say that the potentials vary slowly on the scale set 
by the charge diameter and that hence point charge limit and adiabatic limit 
coincide. 

To see the difference let us regard R p as a small parameter, relative to 
some reference scale. As before we require that the total charge 



d i xp(x)=e (4.40) 
is independent of R p . The electrostatic energy diverges then as 



i J d 3 k\p(k)\ 2 ^ = R 



p 



m e (4.41) 



for small R p , where m e is the electrostatic energy of the charge distribution 
at the reference scale which is independent of R p . In particular, the ratio 
field mass to bare mass grows as i?" 1 in the point charge limit and remains 
constant in the adiabatic limit. 

To display the order of magnitude of the various dynamical contributions 
we resort again to our standard example of an electron in a uniform magnetic 
field B = Bn, n = (0, 0, 1) with B of the order of 1 Tesla = 10 4 Gauss, say. 
It suffices to consider small velocities. In the adiabatic limit we set B = eB 
where the reference field B = 1.1 x 10 17 Gauss, compare with Section fL2 . 
Up to higher order corrections, the motion of the electron is then governed 
by 

(m b + % m e ) i) = - eB (v x n) + v + 0(e 3 ) (4.42) 
6 c one 6 

on the microscopic scale. Going over to the macroscopic time scale, t' = e~ 1 t, 
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qoq ) b 



ecomes 



4 e 6^ 

(m b + - m e ) v = - 5 (« x n) + - — - ev + 0(e 2 ^ 

6 C D7TC 



(4.43) 



Setting mo = + |m e , u;° = eB /m c, j3 = e 2 /6irc 3 m , and restricting 
to the critical manifold, as will be explained in Chapter [7|, Equation (|4.43|) 
becomes 



v = uj®(v x n + (v x n) x n) + C(£ 2 ) , 



(4.44) 



equivalently, on the microscopic time scale 



v = u c (v x n + (3u c (v x n) x n) + 



(4.45) 



with the cyclotron frequency u c = eeB /m c = eB/m c. 

For the point charge limit we rely on the Taylor expansion of Section |5T2 . 
Then, for small velocities, 



(m b + R' 1 - m e ) v = -B{v x n) 



p 3 c Giro 3 



v + 0(iL 



(4.46) 



Since based on the same expansion, as long as no limit is taken, of course, we 
can switch back and forth between (|4.46|) and (|4.42| ), resp. ( |4.43| ), provided 
the appropriate units are used. This can be seen more easily if we accept 
momentarily the differential-difference equation 



m h v{t) = e(E ex (q(t)) + c~ l v{t) x B ex (q(t))) 



UttcR* 



(v(t - 2R p /c) - v(t)) 
(4.47) 



cf. ( |5.14j ), as an approximate equation for the motion of the charge. If we 
expand in the charge diameter R p , then 

2 2 

(m h + g^j) v = e(E ex + c~ 1 vx B cx ) + ^ v + 0{R P ) , (4.48) 



Que 3 
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which is the analogue of (4.46). On the other hand, if we assume that the 



external fields are slowly varying, as discussed in Section flj.ip , then on the 
macroscopic scale 



em h v(t) = ee(E cx (q(t)) + c- 1 v(t)x B cx (q(t))) 

(v(t-2eR p /c)-v(jt)), (4.49) 



e 2 



where R p is now regarded as fixed. Taylor expansion in e yields 

2 2 

(m h + G J R c2 )v = e{E cx + c- l vx B cx ) + e v + 0{e 2 ) (4.50) 



which is the analogue of Q4.43|) . 



As can be seen from ( (4.46Q , in the point charge limit the total mass 
becomes so large that the particle hardly responds to the magnetic field. 
The only way out seems to formally compensate the diverging R~ x (4/3)m e 
by setting m b = —Rp 1 (4/3)m e + m exp . But this is asking for trouble, since 
the energy ( 1.42|) is no longer bounded from below and potential energy 



can be transferred to kinetic mechanical energy without limit. To see this 
mechanism in detail we consider the Abraham model with B cx = and 
4> cx varying only along the 1-axis. The bare mass of the particle is now 
— mb with mb > 0, as before. We set q(t) = (qt, 0,0), v(t) = (v t ,0, 0), 
E cx = (— (p'{q)i 0, 0). <p is assumed to be strictly convex with minimum at 
q = 0. Initially the particle is at rest at the minimum of the potential. Thus 
E(x, 0) = E (x) from (|3.5|) and B(x,0) = 0. We now give the particle a 
slight push to the right, which means q — 0, v > 0. By conservation of 
energy 

-m b c 2 7 (u t ) + e <f>{q t ) + ~ J d 3 x (E(x, tf + B(x, tf) 
= -m h c 2 j{v ) + e<p{q ) + 1 d 3 xE{x,0) 2 . (4.51) 



2 . 

We split E into longitudinal and transverse components, E = Eu+Ej_, E\ 
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\k\ 2 k(k ■ E). Clearly f d 3 x E\\E± = and therefore 

d 3 xE(x,t) 2 > j ' d 3 xE ll (x,t) 2 = j ' d 3 k\k\- 2 (k-E(k,t)) 2 

= J d 3 k\k\- 2 \p(k)\ 2 = J d 3 xE(x,0) 2 , (4.52) 
since the initial field has zero transverse component. Inserting in ( |4.51|) yields 



q 2 t > 1 " [7faO + (e/m b c 2 )(0( % ) - 0( 9o ))]~ 2 ■ (4.53) 

Since j(vq) > 1, q t > for short times. As the particle moves to the right 
(4>(qt) — 0(<7o)) is increasing and therefore q t — > 1 and q t — > 00 as t — > 00. Note 
that t>o an d m b can be arbitrarily small. Not so surprisingly, the Abraham 
model with a negative bare mass behaves rather unphysically. A tiny initial 
kick suffices to generate a runaway solution. 

The point charge limit is honored through a long tradition, which however 
seems to have constantly overlooked that physically it is more appropriate to 
have the external potentials slowly varying on the scale of a fixed size charge 
distribution. Then there is no need to introduce a negative bare mass and 
there are no runaway solutions. 



Notes and References 

ad \4 ■ 4 ' The importance of slowly varying external potentials has been em- 
phasized repeatedly. It is somewhat surprising then, that this notion was 
apparently never transcribed to the equations of motion. In the context 
of charges and the Maxwell field the adiabatic limit was first introduced in 
Komech, Kunze, Spohn (1999) and in Kunze, Spohn (1999). 

ad \4-Q : A more detailed discussion of the hydrodynamic limit can be found 
in Spohn (1991). 

ad |^. j| : The initial slip as discussed here is a side-remark in Komech, Kunze, 
Spohn (1999), where the adiabatic limit for a scalar wave field is studied. 
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Komech, Spohn (1998) prove the long-time asymptotics without the restric- 
tion |e| < e but imposing the Wiener condition instead. Orbital stability is 
established by Bambusi, Galgani (1993). 

ad \4-4\ : In the early work on the classical electron theory, one simply ex- 
panded in R p . R p was considered to be small, but finite, say, of the order 
of the classical electron radius. Schott (1912) pushed the expansion to in- 
clude the radiation reaction which he concluded to be "quite inappreciable 
in this and probably in all practical cases". According to Frenkel (1925) 
the electrodynamics of point, rather than extended, charges is an idea of P. 
Ehrenfest. The point charge limit is at the core of the famous Dirac (1938) 
paper, cf. Section (272"! Since then the limit m\, — ► — oo is a standard piece of 
the theory, reproduced in textbooks and survey articles. The negative bare 
mass was soon recognized as a source of instability. We refer to the review 
by Erber (1961). On a linearized level stability is studied by Wildermuth 
(1955) and by Moniz, Sharp (1977). Bambusi, Noja (1993) discuss the point 
charge limit in the dipole approximation and show that in the limit the true 
solution is well-approximated by the linear Lorentz-Dirac equation with the 
full solution manifold, physical and unphysical, explored. The bound Q4.53 ) 
is taken from Bauer, Diirr (1999), which is the only quantitative handling of 
the instability for the full nonlinear problem. 
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5 Self - Force 



The inhomogeneous Maxwell equations have been solved in ( 1.16Q , ( |1 . 1 7|) . 



Thus it is natural to insert them in the Lorentz force equation in order to 
obtain a closed, albeit memory equation for the position of the particle. 

According to ( |1.16| ), ( |1 . 1 7| ) the Maxwell fields are a sum of initial and 
retarded terms. We discuss first the contribution from the initial fields. By 
our specific choice of initial conditions they have the representation, for t > 0, 

o 

E (0) (x,t) = - J ds J d 3 y(VG^ s (x-y)p(y-q -v°s) (5.1) 

— oo 

+d t G t „ s (x-y)v p(y-q -v°s)), 

o 

B (0) (x,t) = J ds Jd 3 yX7xG t _ s (x-y)v p(y-q -v°s), (5.2) 



compare with (|3.21|) , ( |3.22| ). ( |3.12j ) and ( |1.16| ), can be checked by going 



to Fourier space and using (|1.17|) as initial condition in (|Q| ), (|5.2| ). Since 
G t is concentrated on the light cone, we conclude from (|5.1|), ( p.2| ) that 
E(q)(x, t) = 0, B( )(x, t) = for |q° — x\ <t — R p . If we would have allowed 
for more general initial data, such a property would hold only asymptotically 
for large t. 

Next we note that by energy conservation the particle cannot travel too 
far. Using the bound on the potential, we can find a v < 1 such that 

sup \v £ {t)\ <v<l, (5.3) 

cf. Equation (|5.26| ). The charge distribution vanishes for \x — q{t)\ > R p . 
Since \q{t) \ < v, the initial fields and the charge distribution have zero overlap 
once 

t>t p = 2R p /(l-v). (5.4) 

Thus for t > t p the initial fields make no contribution to the self-force and 
it remains to discuss the effect of the retarded fields. 
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We insert ( |1.12| ), ( 1.13 ) into the Lorentz force equation setting the exter- 
nal potentials equal to zero for a while. Then on the macroscopic scale, for 
t > et p , 



dt 



(m hl rf(t))=Fl u (t) 



(5.5) 



with the self-force 



KeU (t) 



12 p -ik-{qc{t)-q-{s)) 



J dse J d 3 k\p(ek)\ 2 e- 
o 

((|fc| _1 sin \k\(t - s))ik - (cos \k\(t - s))v e (s) 
-({k^ 1 sm\k\(t - s)) v e (t) x (ijfe x v e (s))) . 



(5.6) 



(|5.5|) is exact under the stated conditions on the initial fields. No infor- 
mation has been discarded. The interaction with the field has been merely 
transcribed into a memory term. To make further progress we have to use 
a suitable approximation which exploits that the external forces are slowly 
varying. Since this corresponds to small e, we just have to Taylor expand 
-Fg elf (t), which is carried out in Section ^| with the proper justification left 
for Appendix |5.3| . But before, also to make contact with previous work, we 
have a closer look at the memory term. 
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5.1 Memory equation 



Equation ( |5.4j ) can be further simplified where we set £ — 1 in this subsection. 
By partial integration 

ds J d 3 k\p(k)\ 2 e - iHq{t) - q(s)) v(s) ^{k^ 1 sin\k\(t - s^ 

- J d 3 k \p(k)\ 2 e -<*-(«C*)-*(0)) V (Q) j fc | -1 sin \ k \ t 
t 

-J ds J rf 3 fc|p(fc)| 2 e^ fc - (9(i) - q(s)) (|fc|- 1 sin|fc|(t-s))(i;(s) 



Since t > t p , the boundary term vanishes. Inserting ( |5.7|) into (|5.6|) , returning 
to physical space, and setting t — s = r, we have for t >t p 

oo 

itself (*) = -/ dT [*(* - T ) + (! - V ( f ) ■ V (t - T )) V - ( 5 - 8 ) 


+v(t-r)(v(t) -v(t-r)) ■ V x ]W t (x)\ x=qit y q{t _ T) , 

where, as in ( pL.52|) , 



W t {x) = J d 3 k\p(k)\ 2 e - tk - x \k\- 1 sm\k\t. (5.9) 

In fl5.8|) we have extended the integration to oo, since the integrand vanishes 
for r > t p . Carrying out the angle integrations in ( |5.9| ) we obtain 

W t (x) = |aj| -1 (A(|x|+t) - h(\x\-t)) , (5.10) 

oo 

h(w) = 2tc / dkg(k)coskw (5-H) 



o 

12 



with ( | fe | ) = \p(k)\ . Since p vanishes for \x\ > R p , h(w) = for \w\ > 
2R P . Note that \q(t) — q{t — r)| < vr. Thus for t > t p we indeed have 
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Wt(q(t) — q(t — t)) = 0, as claimed before. -F se if(t) has a finite memory 
extending backwards in time up to t — t p . 

To go beyond (|5.10|) one has use a specific p. Two, at the time popular, 
choices are p s (x) = e(47r-R 2 ) -1 S(\x\ — R p ) and Pb(x) = e (47r_R^/3) _1 for 
\x\ < Rp, Pb{x) = for \x\ > R p . For the uniformly charged sphere one 
finds 

h(RpW)= S eH8«R p )-Hl-H/2) for M<2, 
V 9 ' \ for \w\ > 2, K ' 

and for the uniformly charged ball 

f e>(S«R e rfh*hH for H<2, 
I for \w\ > 2 , 

with /i(u>) = (1 - w 2 )!!^^!} . 

For the charged sphere W t (x) is piecewise linear and, by first taking 
the gradient of W, the time integrations simplify. In the approximation of 
small velocities the motion of the charged particle is then governed by the 
differential-difference equation 

e 2 

m h v(t) = e(E cx (q(t)) + v(t) x B cx (q(t))) + ^-^ (v(t - 2R p ) - v(t)) , 

(5.14) 

where we have reintroduced the external fields. 

The memory equation ( |5.14| ) is of suggestive simplicity. However, to have 



a well defined dynamics one has to prescribe q(0) and v(t) for —2R p < t < 0. 
No instruction for that choice is provided by the supporters of differential- 
difference equations. More importantly, R p is a small parameter and we 
might allow on the top of small velocities a further error of 0(R P ) by Taylor 
expanding in (|5.14j) to obtain 

2 2 

m h v = e(E cx + vx B cx ) - v + |- v + 0{R P ) . (5.15) 

6nR p 6ir 



( 5 . 1 5| ) is a differential equation and only q, v, v are needed as initial data. As 
to be discussed in Section |6|, in fact v(0) is determined by q(0),v(0), since 
the physical solution has to lie on the critical manifold of fl5.15|) . 



5.2 Taylor expansion 



We return to Equation ( |5.5| ). As will be explained in Section |5]^ we know 
that 



sup \q £ (t)\ < C, sup | q (t)\ < C, sup | q £ (t)\ < C (5.16) 



uniformly in e, provided the total charge e is sufficiently small. This smallness 
condition only reflects that at present we do not know how to do better 
mathematically. Physically we expect (|5.16|) to hold no matter how large e. 

Because of (|5.16|) we are allowed to Taylor expand in ( |5.6| ). To simplify 
notation we set v £ {t) = v and t — s = r. Then 



V £ ( S ) = V £ (t-T)=V-VT+-VT 2 + 0(T 3 ), (5.17) 
e -i Htfi t)- q e(s)) = e -*(««(t)-«-(t-r)) = e -i(*.„)T ^1+1 r 2, (fc . _ 1 r 3 - (fe . 

-\ (1 r\k ■ v) - i r 3 (fc ^)) 2 + 0((|fc|r 2 ) 3 )) . (5.18) 
89 



Inserting in ( |5.6| ) and substituting s' = e 1 s, k' = ek yields 



Ft el{ (t) = J dre- 1 J d s k\p(k)\ 2 e- i(k - v)T 
o 

| sin \k\r)ik — (cos \ k\r)(v — erv + - e 2 r 2 v) 



— (\k\ 1 sin \k\r) (v x (ik x v) — v x (ik x erv) + - v x (ik x e 2 r 2 v)) 

+— er 2 i(k ■ v) sin \k\r)ik — (cos |fc|r)(v — erv) 

— (|fc| _1 sin \k\r)(v x (ik x v) — v x (ik x erv))) 

+ ( - ~ £ 2 r 3 «(fc • «) - ~ eV(Jfe • -i)) 2 ) ((Ifc]- 1 sin |jfe|r)iJfe 

-(cos \k\r)v - (Ik]' 1 sin \k\r)(v x (ife x «))) ) + C(e 2 ) . (5.19) 



The terms proportional to e 1 cancel by symmetry. We sort all other 
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terms, 

73k I 5/1,-112 



n*(t) = Jd 3 k\p(k)\ : 

e-H 

{(-{vv)V v + v(v-V v )) J dre-^ 7 \\k\- 1 sin \k\r) 



e~H 



+ (v + ^v(v-V v )) J drre' i(k - v)T (cos \k\r) 
o 

+e(^ [ - {v 2 - l)(v ■ V V )V V + v(v • V„)(t> • V„) + (v ■ v)W v 

-v(v ■ V v )} + i [ - (1 - v 2 )(v ■ V V )V V - v(v ■ W v )(v ■ W v ) 

+3(v ■ v)(v ■ V V )V V - 3v(v ■ V v )(v ■ V w )] + ^ [(v 2 - l)(v ■ V V ) 2 V V 



e~H 



-v(v-V v )(v-V v ) 2 ]j J drre-^ T (\k\- 1 sin \k\r) 

o 

+ e (-*-\ ■ v -) + 3 *>0 • v -)]) 

J dTT 2 e- l(k - v)T cos \k\r} +C(e 2 ). (5.20) 
o 

To take the limit e — > we go back to position space and use the funda- 
mental solution of the wave equation. Then, for p = 0, 1, 

lim J dr J d 3 k\p{k) | 2 e~ %{k - v)T {\k\~ l sin \k\r) r p (5.21) 
o 

oo 

J dt J d 3 x J d 3 yp(x)p(y)^5(\x + vt-y\-t)t p 
o 

Jd 3 k\p(k)\ 2 [k 2 - (k-v) 2 }- 1 for p = 0, 
j d 3 x p(x) J d 3 yp(y) {^ 2 /A-k) for p = 1 . 
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By the same method 



Urn / dr [ d 3 k\p(k)\ 2 e- l{k - v) W 1+p -^-(\k\- 1 sin \k\r) 



(5.22) 



oo 

(l+p+{v- V„)) y J d i k\p{k)\ 2 e- i{k - v) \\k\- 1 sin\k\t)t v 
o 

Jrf 3 A;|p(fc)| 2 (fc 2 + (fc-v) 2 )[fc 2 -(fc-i;) 2 ]- 2 for p = 0, 
/ d 3 x p(x) J d 3 yp(y) for p = 1 . 

Collecting all terms the final result reads 

—rri{(v)v + e(e 2 /6ir) [^(v ■ v)v + 3j 6 (v ■ v) 2 v 

+37 4 (v ■ v)v + 7 2 v] + 0(e 2 ) (5.23) 



with 
mAv) 



v 



-4„,2 



7 2 (3-i; 2 ) - (2\v\ 5 )~ 1 (3 + v 2 ) log 

1 + \v 



1 + 


\v\ 


1 - 


\v\ 



+ (-\v\- 2 + (2\v\ 3 y 1 (l + v 2 ) log 



\v 



)\v)(v\ 
(5.24) 



Note that mi{v) = d(P s — m^v) / as a 3 x 3 matrix. 

Up to order e, F £ sclf (t) consists of two parts with a rather different char- 
acter. The term — rrif(v)v is the contribution from the electromagnetic field 
to the change in total momentum. We computed this term already in Section 



3J]via a completely different route. As emphasized there, since the Abraham 
model is semi-relativistic, the velocity dependence of rrif has no reason to be 
of relativistic form and indeed it is not. The term proportional to e in (|5.23 ) 
is the radiation reaction. Again there is no a priori reason to expect it to be 
relativistic, but in fact it is. Using the four-vector notation of Section [L 
the radiation reaction can be rewritten as 



e(e 2 /6vr)[i; /i - v x vxv>*\ . 



(5.25) 
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5.3 Appendix: How to bound the acceleration? 

We return to the microscopic time scale. From the conservation of energy 
together with condition (P), we have 

E s (v°) + (<f> * p)(eq°) = S(E ,B ,q ,v°)=£(E(t),B(t),q(t),v(t)) 

> m h ^(v(t)) + e min (5.26) 

and therefore 

sup \v(t)\ <v< 1. (5.27) 

t 

The external forces are of order e. Superficially the self-force is of order 
one. However for a Coulombic field the self-force vanishes. Thus if we could 
show that the deviations from the appropriate local soliton field are of order 
e, then altogether 

sup |i>(t)| < (7e (5.28) 

t 

with C a suitable constant. This is what we want to prove. We will not keep 
track of the constants and the value of C changes from equation to equation. 
We make sure however that C depends only on v and is thus determined by 
the initial conditions. Of course, to justify the Taylor expansion of Section 



5T2| , we also need analoguous estimates on higher derivatives, which can be 
obtained with more effort through the same scheme. Here we want to explain 
how to get ( |5.28|) and why we need e to be sufficiently small, at least at 



present. 

^From the equations of motion we have 



v = m (v) 1 e J d 3 xp(x — q) (E ex (ex) + v x B cx (ex)) 

+ J d 3 xp(x-q)(E(x) + vx B(x)) , (5.29) 

where Wq 1 ^) = (mb7) _1 (l — |v|~ 2 |i>) (v\) is the matrix inverse of m (v). 
Clearly by ( |5.27|) we have ||mo(^) — 1 1| < C and, by condition (P), the first 
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term is bounded as 

e\ J d i xp(x - q)(E cx (ex) + v x B cx (ex)) \ < Ce . (5.30) 

On the other hand the self-force looks like order one. To reduce it we have 
to exploit that E, B deviate only little from E v , B v close to the charge 
distribution, i.e. we rewrite the self-force as 

J d 3 x p{x - q) [E(x) - E v (x) +vx (B(x) - B v (x))] (5.31) 

and have to show that the term in the square bracket is of order e. 
Let us define then 

Z(xt)=( E(x,t) - E v(t) (x - q(t)) \ 
\ B(x,t)-B v{t) (x-q(t)) J ' 

Using Maxwell equations and the relations (v ■ V) E v = — V x B v + pv, (v ■ 
V)B V = V x E v we obtain 

Z(t)=AZ(t)-g(t) } (5.33) 

where 

A= ( Vx ) { )= ( (*M-V.)*(. -,«)) V (5 .34) 
V -Vx j'"' 1 \(i,(t)-V v )B v (x-q(t)) ) ' 

Therefore ( |5.33| ) has again the structure of the inhomogeneous Maxwell equa- 
tions. Since by our assumption on the initial data Z(0) = 0, we have 

t 

Z(t) = - J dse A{t - s) g(s). (5.35) 
o 

We set W(t, s) = e A ^~ s ^g(s). Below we prove that 

|Wi(t,s,g(t) + x)\ + \W 2 {t,s,q{t) + x)| < eC|i>(s)|(l + (t - s) 2 )- 1 

(5.36) 
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for | a? | < R p . Therefore inserting in ( |5.29| ) we obtain 

t 

\i){t)\ < eC(e + e J ds{l + {t- s) 2 )" 1 \v(s)\) . (5.37) 
o 

Let k, = sup |v(£)|. Then ( |5.37| ) reads 



t>o 



oo 

k < eC(e + en J ds (1 + s 2 ) -1 ) 



o 

eC 

« < J-Tc^ (5-38) 

^From the computation below we will see that C depends on v (and on model 
parameters like the form factor /), but not on e. Thus taking e sufficiently 
small we can ensure e 2 C < 1 and k < Ce as claimed. 

We still have to establish ( |5.36| ). e At is given in Equation ( |1.12| ), ( |1.13|) . 



Since V • gf 1 (s) = = V • g 2 (s), the term proportional to |fe)(fe| drops out. 
In real space sin \k\t becomes G t from ( |1.15| ) and cos \k\t becomes d t G t - 
Therefore 

W 1 (t,s,x)= A . 1 .„ / d 3 y5(\x-y\-(t-s)) 



4vr(t - s) 2 j 

[(t - s)V x # 2 (y, s) + ^(j/, s) - (a: - y) ■ Vg^y, s)] , 
W 2 (t,s,x)= \ J d 3 y5(\x-y\-(t-s)) (5.39) 

[-(t - s)V x c/^y, s) + g 2 (y, s) - (x - y) ■ Vg 2 (y, s)] . 



We insert g from ([5.34 ). E v and B v are first order derivatives of the function 



<f) v which according to (|3.7f) is given by 

0,(8!) = J d 3 yp(x - 2 /)(47r)- 1 [((1 - v 2 )y 2 + (v y) 2 )]~ 1/2 . (5.40) 
Using (|3.5| ) we have component-wise 



iV^aOI + lV^aOl < C(|V0„(a;)| + |VV„^(a;)|), (5.41) 
\VV v E v (x)\ + \VV v B v (x)\ < C(|VV^(a:)| + iVVV^(cc)|). 
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Taking now successive derivatives in ( |5.40 ) we obtain the bounds 



|V0„(aO| + |VV„<k,(aO| < eC{l + \x\Y\ 
|VV0„(a:)| + |VVV„0„(a;)| < e C (1 + | | ) 3 , (5.42) 



which implies 



\ gi (x,s)\ + \g 2 (x,s)\ < eC\v(s)\(l + \x - q^s)] 2 )- 1 , 
\V gi (x,s)\ + \Vg 2 (x,s)\ < eCl^Kl + la:-^)! 3 )- 1 . (5.43) 



We insert the bound ( |5.43|) in ( |5.39|) which results in a bound on W(t, s, q{t) + 



x). We use that \x\ < R p and \q(t) — q(s)\ < v\t — s\, which finally yields 

(EH)- 

We summarize our findings as 

Theorem 5.1 For the Abraham model satisfying the conditions (C),(P), 
and (I) there exist constants e and C , depending only on the initial conditions 
through v, such that on the microscopic time scale we have 

sup \v{t)\ < v < 1, sup \(^-) n v(t)\ < Ce n } n= 1,2,3, (5.44) 
t t dt 

provided the charge is sufficiently small, i.e. e < e. 

If we would keep track of the constant C, we would get a bound of the 
admissable charge in Theorem |5j]. Since we believe this restriction to be an 



artifact of the method anyhow, there is no point in the effort. 



Notes and References 

ad |5. i| : Sommerfeld (1904,1905) first used systematically memory equations. 
In fact he considered the Abraham model with the kinetic energy m\yV 2 /2 
for the particle and wanted to understand what happens when v(0) > c. He 
argued that the particle rapidly looses its energy to become slower than c by 
emitting what we call now Cherenkov radiation. The differential-difference 
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equation was stated first by Page (1918) with its relativistic generalization 
Caldirola (1956). For reviews we refer to Erber (1961) and Pearle (1982). 
Moniz, Sharp (1974,1977) supplied a linear stability analysis and showed 
that the solutions to ( |5.14j ) are stable provided R p is not too small. For that 
reason Rohrlich (1997) regards ( |5.14| ) and its relativistic sister as the fun- 
damental starting point for the classical dynamics of extended charges. We 
take the Abraham model as the basic dynamical theory. Memory equations 
useful tool in analyzing its properties. 



ad \5. <\ : Taylor expansion is taken from Kunze, Spohn (1999). Such an expan- 
sion was first used by Sommerfeld (1904,1905) and then repeated in various 
disguises. The traditional expansion parameter is the size of the charge dis- 
tribution, which in our context is replaced by the scaling parameter e 



ad \5. 4 ' The bound on v £ {t) comes from Kunze, Spohn (1999) where also 
higher derivatives are discussed. The contraction argument first appeared in 
Komech, Kunze, Spohn (1999). 
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6 Comparison Dynamics 

If in ( [5.23| ) we simply ignore the error of order e 2 , then we obtain the following 
approximate equation for the motion of the charge, 

q = v, m(v)v = e(E cx (q) + v x B cx (q)) (6.1) 
+e(e 2 /6tc) [y 4 (v ■ v)v + 37 6 (t> ■ v) 2 v + 37 4 (t> ■ v)v + 7 2 -i)] . 

Here m(v) is the effective velocity dependent mass. It is the sum of the bare 
mass and the mass ( |5.24|) induced by the field, 

m{v) = m b (7ll + 7 3 |t>) (v\) + m f (f ) . (6.2) 



As anticipated in Section IO|, via a distinct route, the leading contribution 



to ( |6.1|) is derived from the effective Lagrangian 

L eS (q, q) = T(q) - e(0 cx (q) - q ■ A CK (q)) , (6.3) 

equivalently from the Hamiltonian 

E eS (p - eA cx (q)) + e0 cx (q) . (6.4) 

For later purposes it is more convenient to work with the energy function 

H{q,v) = E s {v) + e<f) cx {q) , (6.5) 

which is conserved by the solutions to (|6.1| ) with e = 0, compare with ( |3.14|) . 

The term of order e in ( |6.1| ) describes the radiation reaction. Globally its 
effect can be deduced from the energy balance. We add to (|6.5| ) the Schott 
energy, 

G £ (q, v, v) = H(q, v) - e (e 2 /6tt) ^{v ■ v) . (6.6) 
Then, along the solution trajectories of (|6.1|), 



j t G £ (q, v, v) = -e (e 2 /6vr) [^v 2 + ^{v • v) 2 } . (6.7) 
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Thus G £ is decreasing in time. Integrating both sides of ( p.7| ) we have 

-G e {q(t), v(t),v(t)) + G £ (q(0), v{0), v(0)) 
t 

= e(e 2 /Qn) [ ds [ 7 4 v(s) 2 + 7 6 (v(s) ■ v(s)) 2 ] . (6.8) 



The mechanical energy is bounded from below, but the Schott energy does 
not have a definite sign. If (!) the Schott energy remains bounded in the 



course of time, then 



dt [ 7 4 v(t) 2 + 7 6 (v(t) ■ v(t)) 2 } < oo (6.9) 

o 

which implies 

lim v(t) = 0. (6.10) 

t— »CXD 



Equation ( |6 . 9|) corresponds to the finite energy dissipation ( |1.67|) in Sec- 



tion and we can repeat verbatim the discussion there. In essence the limit 
(|6.10| ) allows only two scenarios. 

(i) lim v(t) = 7^ . This corresponds to a scattering situation where 

t— >oo 

the particle escapes into a region with E cx = = B cx and then travels with 
velocity along a straight line. 

(ii) lim v(t) = 0. This corresponds to a bounded motion where the particle 

t— >oo 

eventually comes to rest. At such a rest point, (q*,0), we have 

V0 cx (q*) = O (6.11) 

by©. 

As noted already in Section |1.5| , in general, ( |6.10| ) and ( |6.11 ) carry too 



little information for determining the true long-time behavior, as can be seen 
from the case of the motion in a uniform magnetic field. 

Unfortunately the energy balance does not tell the full story. As noticed 
apparently first by Dirac (1938), Equation ( |6.1|) has solutions which run away 



99 



exponentially fast. There is no contradiction to (|6.7| ). Since the Schott energy 
does not have a definite sign, in (|6.8|) both G £ (t) and the time-integral diverge 
as t — > oo. The occurence of runaway solutions can be seen most easily in 
the approximation of small velocities, setting B cx = 0, and linearizing ex 



around a stable minimum, say at q = 0. Then ( |6.1| ) becomes 



mv = —muj q + e kmv (6-12) 



with km = e 2 /6n. The three components of the linear equation ( |6.12 ) decou- 



ple and for each component there are three modes of the form e . The char- 
acteristic equation is z 2 = —uj 2 + e kz 3 and to leading order the eigenvalues 
are z± = ±iuo — e (Axg>q/2), z 3 = (1/ek) + 0(1). Thus in the 9-dimensional 
phase space for ( |6.12|) there is a stable 6-dimensional hyperplane, C £ . On C £ 
the motion is weakly damped, friction coefficient e(kuj 2 /2), and relaxes as 
t — > oo to rest at q = 0. Transverse to C £ the solution runs away as e" E *'. 

Clearly such runaway solutions violate the stability estimates ( p,16| ) . Thus 
the full Maxwell-Lorentz equations do not have runaways. They somehow 
appear as an artifact of the Taylor expansion in ( |5.6| ). Dirac simply postu- 
lated that physical solutions must satisfy the asymptotic condition 

lim v(t) = 0. (6.13) 



In the linearized version (|6.12| ) this means that the initial conditions have to 



lie on C e . In Theorem |1.3| we proved the asymptotic condition to hold for the 
Abraham model. Thus only those solutions to ( |6.1| ) satisfying the asymptotic 
condition can serve as a comparison dynamics to the true solution. We then 
have to understand how the asymptotic conditions arises, even better the 
global structure of the solution flow to Equation ( |6.1| ) . 

We note that in ( |6.1| ) the highest derivative is multiplied by a small pref- 
actor. Such equations have been studied in great detail under the heading 
of (geometric) singular perturbation theory. The main conclusion is that the 
structure found for the linear equation ( |6.12j ) persists for the nonlinear equa- 
tion ( |6.1|) . Of course the hyperplane C £ is now deformed to some manifold, 



100 



the critical (or center) manifold. We plan to explain the standard example 
in the following section and then to apply the theory to ( |6.1|) . 



6.1 An example for singular perturbation theory 

As purely mathematical example we consider the coupled system 

x = f(x,y), ey = y-h(x). (6.14) 

h and / are bounded, smooth functions. The phase space is M. 2 . One wants 
to understand how the solutions to (|6.14 ) behave for small e. If we just set 



e — 0, then y = h(x) and we obtain the autonomous equation 

x = f(x,h(x)). (6.15) 

Geometrically this means that the two-dimensional phase space has been 
squeezed to the line y = h(x) and the base point, x(t), is governed by ( |6.15|) . 
{y = h(x),x G R} = Co is the critical manifold to zero-th order in e. 

To see some motion in the phase space ambient to Cq we change from t 
to the slow time scale r = e t. Denoting differentiation with respect to r 
by ', ( |6.14|) goes over to 



x' = e f(x, y) , y ' = y- h(x) . (6.16) 

In the limit e —>■ we now have x' = 0, i.e. x(t) = xq and y' = y — h(xo) with 
solution y{t) = (y — h(x ))e t + h(x ). Thus on that time scale, C consist 
exclusively of repelling fixed points. This is why Co is called critical. The 
linearization at Co has the eigenvalue one transverse and the eigenvalue zero 
tangential to Co- In the theory of dynamical systems zero eigenvalues in the 
linearization are linked to center manifolds and Co is also called the center 
manifold (at e = 0). 

The basic result of singular perturbation theory is that for small e the 
critical manifold deforms smoothly into C e . Thus C £ is invariant under the 



solution flow to ( 6.14 ). Its linearization at (x,y) G C e has an eigenvalue of 
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(9(1) with eigenvector tangential to C £ and an eigenvalue 1/e with eigenvector 
transverse to C £ . Thus for an initial condition slightly off C £ the solution very 
rapidly diverges to infinity. Since Co is deformed by order e, also C e is of the 
form {y = h £ (x), x G R}. According to ( |6.14| ) the base point evolves as 



x = f(x, h £ (x)) . 
Since h £ is smooth in e it can be Taylor expanded as 

0{e m+l ) . 



j=0 



X 



By ( |6.14| ) and ( |6.17| ) we have the identity 

ed x h £ (x) f(x, h £ (x)) = h £ (x) - h(x) . 



(6.17) 



(6.18) 



(6.19) 



Substituting (|6.18|) and comparing powers of e we can thus determine recur- 
sively hj(x). To lowest order we obtain 



ho(x) = h(x) , hi(x) = h'(x) f(x, h{x)) 
and to order e the base point is governed by 

x = f(x, h(x)) + ed y f(x, h(x)) h'(x) f(x, h(x)) . 



(6.20) 



(6.21) 



Given the geometric picture of the center manifold, the stable (not run- 
away) solutions to (|6.14j) can be determined up to a set precision. 



6.2 The critical manifold 



Our task is to cast (6J.) into the canonical form used in singular perturbation 
theory. We set (xi, x 2 ) = x = (q, v) G R 3 x V, y = v G R 3 , 



f{x,y) = (x 2 ,y) G V x 



(6.22) 



and 



9{x,y,e) 



T 2 ^)- 1 ((67r/e 2 ) [m(x 2 )y - F ex (x)\ 

-s [3 7 6 (a; 2 • yf x 2 + 3 7 4 (x 2 ■ y)y}) , (6.23) 
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where 7 = (1 — a^) -1 / 2 as before, F cx (x) = e(E cx (x 1 ) + x 2 x B ex (xi)), and 
is the 3x3 matrix k{v) = 1 + 7 2 |tr) with inverse matrix = 
1 — \v)(v\. With this notation Equation ( |6.1|) reads 

i = /(a;, y), ey = ^(a:, y, e) . (6.24) 

We set h(x) = m(x2)~ 1 F ex (x). Then for e = the critical manifold, Co, 
is given by 

C = {(x, h(x)), x G R 3 x V} = {(q, v, v) : m(v)v = F cx (q, v)} , (6.25) 

which means that the critical manifold for e = is spanned by the solutions 
of the leading Hamiltonian part of Equation ( |6.1| ). Linearizing at Cq the 
repelling eigenvalue is dominated by r )~ 2 K,(x2)~ l m(x2) which tends to zero 
as \x 2 \ — > 1. Therefore Cq is not uniformly hyperbolic, which is needed to 
use the results from Sakamoto (1990). 

To overcome this difficulty we modify g to g s , 5 small, which agrees with 
g on R 3 x {v, \v\ < 1 — 8} x M. 3 and which is constantly extended to values 
\v\ > 1 — 5. Thus for | a^ 2 (t) | < 1 — 5 the solution to x = f, ey = g s agrees 
with the solution to x = f, ey = g. For sufficiently small e the modified 
equation has then a critical manifold C £ with the properties as discussed in 
the example of Section |6.1| . We only have to make sure that the modification 
is never seen. Thus, for the initial condition |v(0)| < v, we have to find a 
5 = 5(v) such that \v(t)\ < 1 — 5 for all times. To do so we need the energy 
balance ( |6.7|) . 

We consider the modified evolution with vector field (/, g s ) and we choose 
the initial velocity such that |v(0)| < v < 1. For e small enough this dynamics 
has a critical manifold of the form v = h e (q,v) and \h e (q,v)\ < c\ — Ci(5). 
We start the dynamics on C e . According to ( |6.7D , for all t > 0, 

G £ {q{t),v(t),h £ {t)) < G £ (0) = H(q(0),v(0))-e(e 2 /6n)(v(0)-h £ (0)) 

< + e0 HC (qf(O))+eci. (6.26) 

We choose now 5 such that v < 1 — 25. Since the initial conditions are 
on C £ , the solution will stay for a while on C £ until the first time, r, when 
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|u(t)| = 1 — 5. After that time the modification becomes visible. At time r 
we have, using the lower bound on the energy and ( |6.26 ), 

E s (v(r))+e0 < H(q(r) 1 v(r)) = G £ (r)+e(e 2 /Q7r) 1 4 (v(r)-h £ (r)) 

< E s {v)+e<f) cx (q(0)) + 2ec l (6.27) 

and therefore 

E s (l -5)< E s {\ - 26) + e (0 ex (q(0)) - 0) + 2e Cl . (6.28) 
E s (l — 5) = 1/y/S for small 5, which implies 

^<c 2 + 4eci (6.29) 
V o 

with c 2 = 2e (0(q(O)) — 0). We choose now 5 so small that l/y/5 > c 2 + 1 
and then e so small that Aeci < 1. Then Q6.29Q is a contradiction to the 
assumption that \v(t)\ = 1 — 5. We conclude that r = oo and the solution 
trajectory stays on C e for all times. 

Equipped with this information we have for small e the critical manifold 

v = h £ (q,v). (6.30) 

On the critical manifold the Schott energy is bounded and from the argument 
leading to ( |6.10p we conclude that Dirac's asymptotic condition holds on C £ . 
On the other hand, slightly off C £ the solution diverges with a rate of order 
\je. Therefore the asymptotic condition singles out, for given q(0), v(0), the 
unique v(Q) on C £ . 

The motion on the critical manifold is governed by an effective equation 
which can be determined in approximation as in Section |6.1| . We define 

h{q, v) = m(v)- 1 e(E cx {q) + vx B cx {q)) . (6.31) 

Then, up to errors of order e 2 , 

m(v)v = e (E cx (q) + vx B cx (q)) (6.32) 
+e (e 2 /6n) [j 2 k(v) (v-V q + h- V v h + (3 7 6 (i; ■ hfv + 3 7 4 (v ■ h)h))] . 
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The physical solutions of (|6.1|), in the sense of the asymptotic condition, 
are governed by Equation (|6.32|) , which thus should be regarded as the true 
comparison dynamics to the microscopic equation ( [4.521) . Note that the error 
made in going from ( |6.1| ) to (|6.32| ) is of the same order as the error made in 
the derivation of Equation (|6.1|). 

On a formal level (|6.32| ) is easily deduced from (|6.1| ). We regard m{v)v = 
e(E cx (q) + v x B cx (q)) as the "unperturbed" equation and substitute for 
the terms inside the square bracket, which means to replace v by h and v by 
h = (v ■ V q )h+ (h ■ V v )h. While yielding the correct answer, one misses the 
geometrical picture of the motion in phase space and of the critical manifold. 

For a numerical integration of the comparison dynamics it is advantageous 
to use directly (|6.32j ). The only other practical option would be to solve (|6.1|) 
backwards in time. Then the trajectory is pushed rapidly towards the critical 
manifold. On C e one solves however the time-reversed dynamics which means 
a final rather than an initial value problem. Instead of weakly damped the 
motion is now slowly accelerating. 



6.3 Tracking of the true solution 

^From ( [4 . 1 1| ) we have the true solution q £ (t),v £ (t) with initial conditions 
q°, v° and correspondingly adapted field data. We face the problem of how 
well this solution is tracked by the comparison dynamics (|6.1|) . Let us first 
disregard the radiation reaction. From our a priori estimates we know that 

q £ = v £ , m(v £ )v £ = e (E cx (q £ ) + v £ x B cx (q £ )) + 0{e) (6.33) 

which should be compared to 

r = u , m{u)u = e (E cx (r) + u x B cx (u)) . (6.34) 

We switched to the variables r, u instead of q, v so to more clearly distinguish 
between the true and comparison dynamics. 
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Theorem 6.1 For the Abraham model satisfying the conditions (C), (P), 
and (I) let e < e and e < so be sufficiently small. Let r(t),u(t) be the 
solution to the comparison dynamics ( \6. 34 ) with initial conditions r(0) = 
q°, it(0) = v°. Then for every r > there exist constants c(r) such that 



\q £ (t) - r(t)\ < c{t)e , \v e {t) - u(t)\ < c{r)e (6.35) 



for < t < r. 



Proof. One converts ( |6.33j ), ( |6.34| ) into a first order equation in its integral 
form. The difference is then estimated by using Gronwall's lemma, which 
yields an error as ee ct . □ 



Theorem |6.1| states that, up to an error of order e, the true solution is 
well approximated by the Hamiltonian dynamics (|6.34j) . In the next order 
the comparison dynamics is 

r £ = u e , m(u £ )u £ = e (E ex (r £ ) + u £ x B cx (r £ )) (6.36) 
+e{e 2 /67c) [ 7 4 {u £ ■ u £ )u £ + 3 7 6 (n £ • ii £ f u £ + 3 7 4 (n e ■ ii £ )ii £ + 7 2 u £ ] 

restricted to its critical manifold C e and one might expect that 

\q £ (t) - r £ (t)\ + \v £ (t) - u £ (t)\ = 0(e 2 ) . (6.37) 

Because of the improved precision one has the possibility to resolve the ra- 
diation reaction correction to ( |6.36|) . 

An alternative option to keep track of the correction would be to consider 
longer times, of the order e~ 1 t on the macroscopic time scale. Then the 
radiative effects add up to deviations of order one from the Hamiltonian 
trajectory. Thus 

\q £ (t) -r £ (t)\ = 0(e) for 0<t<e -1 r. (6.38) 

One should be somewhat careful here. In a scattering situation the 
charged particle reaches after a finite macroscopic time the force free re- 
gion. According to Q6.37 ) the error in the velocity is then 0(e 2 ) which builds 
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up to an error in the position of order e over a time span e^ 1 t. Thus we 
cannot hope to do better than ( |6.38 ). On the other hand when the motion 
remains bounded, as e.g. in a uniform external magnetic field, the charge 
comes to rest at some point q* in the long time limit and the rest point q* is 
the same for the true and the comparison dynamics. At least for an external 
electrostatic potential with a discrete set of critical points we have already es- 
tablished such a behavior and presumably it holds in general. Thus for large 
t we have q e (e~ 1 r) = q* and also r £ (s~ l r) = q*. Therefore we conjecture 
that ( |6.38 ) holds for all times. 

Conjecture 6.2 For the Abraham model satisfying (C), (P), and (I) let 
q e (t) be bounded, i.e. \q £ (t)\ < C for all t > 0, e < 6q . Then there 
exists (r £ (0) ,u £ {0), 
u £ (0)) G C £ such that 

sup \q £ (t) - r £ (t) \ = 0(e) , (6.39) 

where r £ (t) is the solution to ( 6. 3(\ ) with said initial conditions. 

At present we are far from such strong results. The problem is that an 
error of order e 2 in (|6.36|) is generically amplified as e 2 e t ^ £ . Although such 
an increase violates the a priori bounds, it renders a proof of (|6.39|) difficult. 
We seem to be back to ( |5.35[ ) which carries no information on the radiation 
reaction. Fortunately the radiation correction in ( |6.36j ) can be seen in the 
energy balance. 

Theorem 6.3 Under the assumptions of Theorem |^J| we have 

\{E s {v £ {t)) + e<P^{q £ {t))} - [E s (u £ (t)) + e <J) ex (r £ (t))}\ < Cc(r)e 2 (6.40) 

for < t < t, where r £ (t),u £ (t) is the solution to \6. 3dj ) with initial data 
r £ (0) = q°, u £ (0) = v°, ii £ (0) = h £ (q°,v°). 

Proof. We use the estimate ( |5.23| ) on the self-force, where \f £ (t)\ < Ce 2 for 
ei p < t and \f £ (t)\ < Ce for < t < ei p . Then, as in ( |6T7[) , 



j f G e (q £ , v £ , v £ ) = f{t) -v £ -e (e 2 /6vr) [ 7 4 (^) 2 + 7 V • v £ ) 2 } (6.41) 
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and therefore 



\H(q £ ,v s )-H(r*,u°)\ 

< e (e 2 /6n)\-f(v £ )\v £ ■ v £ ) - 7 (u £ ) 4 (u £ • ii £ )\ 
t 

+ J ds (|/ £ -v £ \+e (e 2 /67r)| 7 (v £ ) 4 {v £ f + 7 (« £ ) 6 {v £ ■ v £ f 
o 

- 7 (n £ ) 4 (ii £ ) 2 - 7 (u £ ) 6 (u £ ■ ii £ ) 2 \) . (6.42) 

Since \v £ \, \u £ \ remain bounded away from 1, we can use the bound \v £ (t) — 

u (t)\ < C ( T ) £ from Theorem ETT[ Reinserting ( |6.35[ ) into ( |6.33| ) and (|6.34j ) we 

t 

obtain \v £ (t) —ii £ (t)\ < c{r)e. Furthermore J ds |/ e (s)| < Cte 2 . We conclude 

o 

that 

\H(q £ (t), v £ (t)) - H(r £ (t), u £ (t))\ < C(t + c(t))e 2 . (6.43) 

□ 



6.4 Electromagnetic fields in the adiabatic limit 

So far we have concentrated on the Lorentz equation with retarded fields and 
have obtained approximate evolution equations for the charged particle. Such 
an approximate solution can be reinserted into the inhomogeneous Maxwell 
equations in order to obtain the electromagnetic fields in the adiabatic limit. 

As before, let q £ (t),v £ (t)),t > 0, be the true solution. We extend it to 
<f (t) = q° + v °t, v £ {t) = v° for t < 0. According to (|372T1) , (|3T22D and using 
the scaled fields as in ([Of ), we have 



(6.44) 



t 

E(t) = - J ds Jd 3 y (VG t _ s * p £ (s) + d t G t - s * j £ (s)) 
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with p £ (x,t) = p £ (x — q £ (t)),j £ (x,t) = p £ (x — q £ {t))v £ {t). Inserting from 
( |1 . 1 5| ) and by partial integration 

t 

-j=E(x,t) = - J ds j d 3 y 5{\x -y\-(t- a)) Vp £ {y, s) 

— oo 

t 

- j ' ds J d 3 y s y 6(\x-y\-(t- s)) [{y - x) ■ Vj e {y, s) + j £ (y, s)} 

— oo 

= - / d ' 3 y f tt- — i v p^y - ftt ~\ x ~ -\ x ~v\) ( 6 - 45 ) 

J \47r|£c — y\ 

+ -TT^ i? ~\ x ~ 2/IX 1 + (y~x)-V) p £ (y - q £ (t -\x- 2/|))) . 



47t\x — y 
In the same fashion 

-1= B(x, t) = - / rf 3 y — , v £ (t - | a; - 2/|) x Vp £ {y - q £ {t -\ x -y\)). 

Je J An\x-y\ 

(6.46) 



In the limit e — > we have p e (£c) — > 5(cc) and, by Theorem |6.1] , <? £ (t) — > 
r(t), v e (t) -> u(t), where r(t) = <?° + w°t, u(t) = v° for t < 0. We 
substitute y' = y — q £ {t — \x — y\) with volume element det(dy/dy') = 
[1 — v £ (t — \x — y\) ■ (x — y)/\x — 2/|] _1 . Then S(y') leads to the constraint 
= y — r(t — | x — y\) which has the unique solution y = r(t TCt ), compare 
with ( |1.19|) . In particular the volume element det(dy / dy') becomes in the 
limit [1 - n • u^)]- 1 with n = n(x,t) = (x - r(t Ict ))/\x - r(t iet )\. 
We conclude that 

]jm -=E{x,t) = E(x,t), (6.47) 

e^O Je 



lim -^=B(x,t) = B(x,t), (6.48) 

e^O Je 



where E , B are the Lienard-Wiechert fields ( |1.21| ) , ( ]l.22 ) generated by a 



point charge moving along the trajectory t i— > r(t). The convergence in 
(|6~47p , is pointwise if one excludes the Coulomb singularity at x = **(£)■ 
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6.5 Larmor's formula 



We want to determine the energy per unit time radiated to infinity and 
consider, for this purpose, a ball of radius R centered at q £ {t). At time t + R 
the energy in this ball is 

(* + R ) = 8 (°) - \ J d 3 x(E(x,t + R) 2 + B(x,t + R) 2 ) 

{\ X -qe(t)\>R} 

(6.49) 

using conservation of total energy. The radiation emitted from the charge at 
time t reaches the surface of the ball at time t + R and the energy loss per 
unit time is given by 
d 

(n(x) ■ v £ (t)) + E(x, t + R)- [n(x) x B(x, t + R)] 
= \R 2 J d 2 uo ((E{q £ {t) + Ru,t + Rf + B(q £ (t) + Ru,t + R) 2 ) 

(u> ■ v £ (t)) + 2E{q £ {t) + Ruj,t + R)-[ux B{q £ {t) + Ru,t + R)]) , 

(6.50) 

where n(x) is the outer normal of the ball and |u?| = 1 with d 2 u the integra- 
tion over the unit sphere. (|6.50 ) holds for sufficiently large R, since we used 
that {x\ \x-q £ (t)\ > R}n{x\ \x-q £ (t + R)\ < eR p } = 0, which is the case 
for (l-v)R > eR p . 

( |6.50|) still contains the reversible energy transport between the consid- 
ered ball and its complement. To isolate that part of the energy which is 
irreversibly lost we have to take the limit R — > oo. For this purpose we first 
partially integrate in ( |6.45|) , (|6.46| ) by using the identity 

Vp = V yP - (1 + fozjOlg )- 1 (^ . S7 y )p (6.51) 

\y - x\ K \y - x\ J y 
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at the argument y — q £ (t — \y — x\). For large R the fields in ( 6. 50 ) become 
then 



RE(q £ (t) + RuJ,t + R) = ^ J d 3 y±-p £ (y-q £ ) 

[ _ (i _ a, . v eyi i,e _ (1 _ w . „ e) -2 (w . ^ £)(v£ _ w) j 1^.^^ (6.52) 

i?B(q £ (t) + Ru, t + R) S ^ y rf 3 y -L p e (y - q £ ) 

[ - (1 - u> • v 6 ) -1 (w x v £ ) - (1 - cj • v £ y 2 (u ■ i> £ )(u> x v £ )] \ t+w<y _ q e m 
= u> x RE(q £ (t) + Rw,t + R) , (6.53) 



where we used that t + R - \q £ {t) + Bw - y\ = t + uj ■ (y - q £ (t)) + 0(1/ R) 
for large R. Inserting in ( |6.50| ) yields 



Km I R>E (t) = I e (t) 

R— >oo 

= -e J d 2 u (1 - w • ^ £ (t)) (RE(q £ (t) + Ruj,t + R)) 2 (6.54) 

= - e y d 2 * (1 ~ " ■«'(*)) ([^ | Aftd/-g E )(l-^^ 2 (^« £ )] 2 
-[^ J d*yp e (y-q £ )(l-u J .v £ )- 1 v £ 

+(!-«■ • i;> £ ] 2 ) | t+ „ W(i)) . (6.55) 



I e (£) is the energy radiated per unit time at e fixed. As argued before it is 
of order e. The expression (|6.54 ) shows that I e (t) < 0. 



( J6.55 ) is not yet Larmor's formula. For this we have to go to the adiabatic 
limit e — > 0. Then q £ (t) — > r(t). Since p e (cc) — > e5(a:) we have = g £ (i) — 
r(t) in ( |6.55 ). From the d 3 y volume element we get an additional factor of 
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(1 -w« e ) -1 . Thus 
lim/ £ (t) = I(t) = -e 2 J d 2 u;(l-u;-u(t))(4ir(l-u>-u(t)y 3 ) 2 

((w ■ u(t)) 2 - [(1 - w ■ u(t))ii(t) + (u> ■ u(t))w(t)] 2 ) 
= -(e 2 /67r) [ 7 4 ^W 2 + 7 6 M*) • ™W) 2 ] 
= -(e 2 /67r)7 4 [tt(t) 2 - (u(t) x u(t)) 2 ] , (6.56) 

which is the standard textbook formula of Larmor. Note that the same 
energy loss per unit time was obtained already in (|6.7|) using only the energy 
balance for the comparison dynamics. 

Starting from (|6.50| ) we could alternatively first take the limit £ _1 /R i£ (t) — > 



lR,o{t), which is the change of energy in a ball of radius R centered at the 
particle's position r(t) in the adiabatic limit. As before we have to isolate 
the irreversible energy loss through 

lim I Rfi {t) = I(t) . (6.57) 

R— >oo 

The energy loss does not depend on the order of limits, as it should be. 
Notes and References 

ad The radiation damped harmonic oscillator is discussed in Jackson 
(1999) with a variety of physical applications. The asymptotic condition is 
first stated in Dirac (1938). It has been reemphasized by Haag (1955) in 
analogy to a similar condition in quantum field theory. 



ad \6. j : Singular, or geometric, perturbation theory is a standard tool in the 
theory of dynamical systems. We refer to Jones (1995) for a review with many 
applications. In the context of synergetics, Haken (1983), one talks of slow 
and fast variables and the slaving principle, which means that fast variables 
are slaved by the slow ones. Within our context this would correspond to 
an attractive critical manifold. The renormalization group flows in critical 
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phenomena have a structure similar to the one discovered here. The critical 
surface corresponds to critical couplings which flow then to some fixed point 
governing the universal critical behavior. The critical surface is repelling 
and slightly off that surface the trajectory moves towards either the high 
temperature or low temperature fixed points. 



ad \6. jj : Particular cases have been studied before, most extensively the one- 
dimensional potential of finite width and with linear interpolation, Haag 
(1955), Carati, Galgani (1993), Carati et al (1995), Blanco (1995), in addition 
head on collision in the two-body problem, Huschilt and Baylis (1976), and 
motion in a uniform magnetic field, Endres (1993). These authors emphasize 
that there can be several solutions to the asymptotic condition. From the 
point of view of singular perturbation theory such a behavior is generic. If 
e is increased, then the critical manifold is strongly deformed and no longer 
given as a graph of a function. For specified q(0), v(0) there are then several 
v(0) on C £ which means that the solution to the asymptotic condition is not 
unique. However these authors miss to underline that the nonuniquess in the 
examples occurs only at such high field strengths where a classical theory 
has long lost its empirical validity. At moderate field strengths the worked 
out examples confirm our findings. The applicability of singular perturbation 
theory for a general class of potentials is first recognized in Spohn (1998). 



ad \6. 3j , \6-A , vA : These results are adapted from Kunze, Spohn (1999). 
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7 The Lorentz— Dirac Equation 



In relativistic notation the Lorentz-Dirac equation reads 

m = (e/c)F^(z)v u + (e 2 /Qnc 3 ) [«" - c~ 2 v x v x v»] , (7.1) 

where we reintroduced the speed of light, c. mo is the experimental rest mass 
of the charged particle. is the electromagnetic field tensor of the external 
fields. In this section we omit the index "ex" for better readability. Formally, 
Equation (7.1) can be derived from the Lorentz model in the adiabatic limit. 



To conform with the usual notation we have set the adiabatic scale parameter 
e=l. But it should be kept in mind that the radiation reaction in (7.1) is 
a small correction to the Hamiltonian part. 

If we fix a frame of reference and go over to three-vectors, then the 
Lorentz-Dirac equation becomes 

mQ^K{v)v = e(E(q) + c~ x v x B(q)) 

+ (e 2 /6vrc 3 ) 7 2 k(«) [v + 3 7 2 c" 2 (v •«)«], (7.2) 

with the 3x3 matrix k{v) = 1 + c~ 2 7 2 |v) (u| and its inverse k{v)~ 1 = 1 — 
c~ 2 \v)(v\. The Lorentz-Dirac equation ( |7.1| ) differs from ( |6.1D only through 
a proper relativistic kinetic energy. Clearly, qualitative properties of the 
solution flow should not depend on such a detail. In the analysis of ( |6.1|) we 
only used the critical manifold to be a uniform repeller, except for \v\/c close 
to one, and the energy balance 

4- (m c 2 j(v) + e(p(q) - (e 2 /67rc 3 ) ^(v ■ v)) 
at v 

= -(e 2 /6vrc 3 ) 7 4 (v, k{v)v) , (7.3) 

which ensured that the solution stays for all times on the critical manifold, 
provided the radiation reaction term is sufficiently small. Thus we can follow 



the blueprint of Section 6.1 to obtain an effective second order equation for 



the motion on the critical manifold. 
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In Section 17^ we work out some examples of experimental interest. While 



at present an actual test is rather indirect at best, the examples should 
convince the reader that the effective second order equation can be handled 
with ease. As an extra bonus we will make some predictions on the motion 
of the charge which could not have been guessed on the basis of Larmor's 
formula. 



7.1 Critical manifold, the Landau— Lifshitz equation 

We write ( |7.2| ) in the standard form of singular perturbation theory, compare 
with Section |Q[ Then 

x = f(x,y), ey = g{x,y,e) (7.4) 



with 



f(x,y) = (x 2 ,y), (7.5) 
g{x,y,e) = (6vrc 3 /e 2 ) (m 7 _1 y - e^- 2 K(x 2 )~ 1 (E(xi) + c -1 ^ x B{x x ))) 
-3e-f 2 c~ 2 (x 2 ■ y)y . (7.6) 



To conform with ( |6.1| ) we reintroduced the small parameter e. At zeroth order 
the critical manifold is {y = h(x)} with h(q, v) = (e/m )7 _1 k(w) -1 (E(q) + 
cr x v x B(q)). Linearizing (|7.5| ), (|7.6| ) at y = h(x) the repelling eigenvalue is 
(67rc 3 /e 2 ) mo7 _1 + 0(e), which vanishes as \v\/c — > 1. Thus we have to rely 
on the same construction as in Section B.2I. 



To order e the effective second order equation is given by QS.32D , except 
that now m(v) = mo 7k(u). We work out the various terms and set e = 1. 
Then the motion on the critical manifold of the Lorentz-Dirac equation is 
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governed by 



mo 7 k(v)v 



v . 

e(E + c~ x v x B 

„2 



67TC 3 



m 



7 



(v ■ V Q )(E + c^v x B) + (— ) 2 c" 1 ((£ x B) 



m ' 

2 



+c~ 1 (i; • E)E + c -1 ^ • 5)5 + ( - E 2 - Br 
+c' 2 (v ■ Ef + c- 2 (v ■ B) 2 + 2c~ 1 v ■ (E x B^c^v) 1 (.7.7) 



While singular perturbation theory provides a systematic method, Equa- 
tion (|7.7|) can also be derived formally. In (7J3) we regard k(v)v = 
e (E + c~ x v x B) as unperturbed equation, differentiate it once, and sub- 
stitute v inside the square brackets of (|7.2|) . Resubstituting v from the 
unperturbed equation results in Equation (|7.7|) . This argument is carried 
out more easily in the covariant form of the Lorentz-Dirac equation. The 
unperturbed part is 



m v» = ( e /c)F> w (z)v v 
and differentiating with respect to the eigentime, 

(m c/e)i^ = v x d x F^(z)v l/ + F^(z)v u 
Substituting Q and Q in Q yields 



(7.* 



(7.9) 



C " 67TC 3 



m c 



-v x d x F' lu v 



m c 



(7.10) 



Written in three-vectors Equation ( [7.1 0| ) coincides with ( |7.7p together with 
the equation for the energy balance. 

Of course, the justification of Equation ( |7.7| ) comes only from the struc- 



ture of the solution flow to (72). Higher order corrections, although rather 
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unimportant in our context, would have to be computed by the method ex- 
plained in Section |6TT| . 

Equation CffTlOD appears for the first time in the second volume of the 



Landau-Lifshitz Course in Theoretical Physics. It seems to be appropriate 
to call then Equation (|7.10|) the Landau-Lifshitz equation. The error in going 
from (|7. 1| ) to (|7.10| ) is of the same order as the one in the derivation of the 



Lorentz-Dirac equation itself. Thus we regard the Landau-Lifshitz equation 
as the effective equation governing the motion of a charged particle in the 
adiabatic limit. 

7.2 Some applications 

(i) Zero magnetic field. For zero magnetic field the Landau-Lifshitz equation 
simplifies to 

2 

m jK{v)v = eE+ [— j(v ■ V q )E (7.11) 

D7TC mo 

+ (_£_) 2 (( 1 , . E )E - 7 2 E 2 v + 7 V 2 (<; • Efv)} . 

ITLqC 

Of interest is a central potential. We set q = r, \r\ = r, f = r/\r\,(f) ex (q) = 
<p(r) which implies E = —<p'r. Then (|7.11 ) becomes 



6 2 C 

rriQ^ k,{v)v = —e<f)'f-\ [ — 7 (— (v ■ f)<f>"r 

one 6 mo 

— (v - (v ■ r)r)(j}' + ( ) 4>' 2 ((v ■ r)f - 7 2 v 

r itlqc 

+7 2 c" 2 (v -ffv)] . (7.12) 
The angular momentum L = r x mo^jv satisfies 

L = ^3 [ - — - 0' - {—Y 7 2 (1 - c~\v f ) 2 )0' 2 ] L . (7.13) 
ottc 6 m r ttiqC 

Thus the orientation of L is conserved and the motion lies in the plane 
orthogonal to L. No further reduction seems to be possible and one would 
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have to rely on a numerical integration. Only for the harmonic oscillator, 
4>{r) = ~ m WQr 2 , a closed form solution can be achieved. 

A somewhat more tractable case is to assume that ex varies only along 
the 1-axis. Setting v — (v, 0, 0), q = (x, 0, 0), and E = (—0', 0, 0), Equation 



(|7.11|) becomes 



6 2 6 

m o 7 3 v = — e<p'(x) — jd>"(x)v. (7-14) 

bite 6 m 



The radiation reaction is proportional to —<f)"(x)v, which we can be regarded 
as a spatially varying friction coefficient proportional to <p"(x). For a con- 
vex potential, <ft" > 0, like an oscillator potential, this friction coefficient is 
strictly positive and the resulting motion is damped until the minimum of 4> 
is reached. In general however, <fi" will not have a definite sign, like the double 
well potential 4>(x) ~ (x 2 — l) 2 , or the washboard potential <fr(x) ~ — cosx. 
At locations where <p"(x) < one has antifriction and the mechanical energy 
increases. This gain is always dominated by losses as can be seen from the 
energy balance 

d 6 2 c 

— [mo 7 + e0 + — -3 — 7 <j>'v\ 

at birc 6 m 

= ~ (-) 2 <P' 2 " — (/» -) 2 7^"- • (7-15) 
birc 6 mo mo bnc A ttlq 

The last term in ( |7.15| ) does not have a definite sign. But its prefactor is 
down by one order in e and therefore it is outweighed by — <f)' 2 . 

Equation ( |7.14p has one peculiar feature. If <p(x) = —a x, a > 0, over 
some some interval [a_,a + ], then 0" = over that interval and the friction 
term vanishes. The particle entering at a_ is uniformly accelerated to the 
right until it reaches a+. From Larmor's formula we know that the energy 
radiated per unit time equals (e 2 /67rc 3 )(e/mo) 2 aQ. This energy must come 
entirely from the near field without a mechanical contribution. The same 
behavior is found for the Lorentz-Dirac equation. If, locally, E = const and 
B = 0, then the Hamiltonian part is solved by the hyperbolic motion, i.e. 
a constantly accelerated relativistic particle. For this solution the radiation 
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reaction vanishes which means that locally the critical manifold happens to 
be independent of e. The radiated energy originates from the near field only. 

(ii) Zero electrostatic field and constant magnetic field. We set B = (0, 0, B) 
with constant B. Then ( |7.7| ) simplifies to 

2 

m olK (v)v = -(v x B) + (—) 2 \(v ■ B)B- 1 2 B 2 v (7.16) 

c 67rcr v m c 

+7 2 c~ 2 (u- B) 2 «] . 
We multiply by and obtain 

2 

m o7 i, = £(« x B) + (— ) 2 [(v ■ B)B - B 2 v] . (7.17) 

C 07TC v m c 

The motion parallel to -B decouples with = 0. We set t> 3 = and v = 
(u, 0), n.^ = (— -U2)Wi). Then the motion in the plane orthogonal to B is 
governed by 

7 it = uj c (u ± — (3u c u) , (7-18) 



with cyclotron frequency u c = eB/moc and (3 = e 2 /67rc 3 mo. (|7.18|) holds 



over the entire velocity range. For an electron (3uj c = 8.8 x 10 _18 B [Gauss]. 
Thus even for very strong fields the friction is small compared to the inertial 
terms. 

(|7.18|) can be integrated as 



| 7 = -^ c 2 ( 7 2 -l) (7.19) 

with solution 

Jt = [ 7o + 1 + (70 - l)e- 2/3 ^][7o + 1 - (7o - lK 2 ^]- 1 , (7.20) 

which tells us how u(t) 2 shrinks to zero. To determine the angular depen- 
dence we introduce polar coordinates as u = u(cosp>, simp). Then 

du „ dip 1 ,„ 
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Thus u(ip) shrinks exponentially, 

u ((p) = u(0) . (7.22) 

Since f3uo c = 8.8 x 10~ l8 B [Gauss] for an electron, even for strong fields the 
change of u in one revolution is tiny. 

To obtain the evolution of the position q = (r, 0), \r\ = r, we use that for 
zero radiation reaction, /3 — 0, 

r = — 7. (7.23) 



UJ, 



c 



By ( |7.22| ) this relation remains approximately valid for non-zero (3. Inserting 
u{t) from ( |7.20|) we obtain 



r(t) = r e-^ 2 '[l + (( 7o - 1)/2)(1 - e" 2 ^)]-! ( 7 .24) 

with r the initial radius and u(0)/c = (70 — l) 1 / 2 /7 the initial speed which 
are related through ( |7.23| ). In the ultra-relativistic regime, 70 ^> 1, and for 
times such that Pu^t 1, (|7.24|) simplifies to 

1 + 7o/3w# 

and the initial decay is according to the power law t~ l rather than exponen- 
tial. 

For an electron (3uj1 = 1.6 x 10~ 6 (B [Gauss]) 2 /sec. Therefore if we choose 
a field strength B = 10 3 Gauss and an initial radius of r = 10 cm, which 
corresponds to the ultra-relativistic case of 7 = 6 x 10 4 , then the radius 
shrinks within 0.9 sec to r(t) = 1 /im by which time the electron has made 
2 x 10 14 revolutions. 

(iii) The Penning trap. An electron can be trapped for a very long time in the 
combination of a homogeneous magnetic field and an electrostatic quadrupole 
potential, which has come to be known as a Penning trap. Its design has been 
optimized towards high precision measurements of the gyromagnetic (/-factor 
of the electron. Our interest here is that the motion in the plane orthogonal 
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to the magnetic field consists of two coupled modes, which means that the 
damping cannot be guessed by pure energy considerations using Larmor's 
formula. One really needs the full power of the Landau-Lifshitz equation. 
An ideal Penning trap has the electrostatic quadrupole potential 

e<j>{x) = ^muj 2 z (~xl - l - x\ + x 2 3 ) , (7.26) 

which satisfies A0 = 0, superimposed with the uniform magnetic field 

B = (0,0,B). (7.27) 

The quadrupole field provides an axial restoring force whereas the magnetic 
field is responsible for the radial restoring force, which however could be 
outweighed by the inverted part of the harmonic electrostatic potential. 

We insert E = — V0 and B in the Landau-Lifshitz equation. The terms 
proportional to {v ■ V q )E, E x B, (v ■ B)B, and B 2 v are linear in v, resp. 
q. The remaining terms are either cubic or quintic and will be neglected. 
This is justified provided 

^ < 1 (7.28) 

c 

and 

(m u 2 Je) r max < B, i.e. r max < c{u c /u 2 z ) , (7.29) 

if r max denotes the maximal distance from the trap center. With these as- 
sumptions the Landau-Lifshitz equation decouples into an in-plane motion 
and an axial motion governed by 

it = ^r + u^ - P^-^u^u+^u^ 2 ^] , (7.30) 
z = —uj 2 z — (3uj 2 z. (7-31) 

Here q = (r,z), v = (u, z), (x u x 2 )^ = (-x 2 , xj. 

The axial motion is just a damped harmonic oscillator with frequency uj z 
and friction coefficient 

lz = [3uj 2 z . (7.32) 
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The in-plane motion can be written in matrix form as 

^ = (A + (3V)i/j (7.33) 

with if) = (r,u) and An = 0, A 12 = 1, A 2 i = A 22 = iou z a y , Vu = 

0, V12 = 0, V 2 i = iu) c u) 2 z a y) V 22 = (cu z —u! 2 )l, where a y is the Pauli spin matrix 
with eigenvectors x±, o y x± — ±X±- The unperturbed motion is governed 
by the 4x4 matrix A. It has the eigenvectors ip+,± = (±i(l/uj + )x^, Xt) 
with eigenvalues ±iu + and i>-,± — (^(l/u^x^, Xt) with eigenvalues ±iu_, 
where 

u ± = l - (u c ± ^M^A ) . (7.34) 

The mode with frequency uj + is called cyclotron mode and the one with oj_ 
magnetron mode. Experimentally u c ^> u z and therefore u + <C The 
orbit is then an epicycle with rapid cyclotron and slow magnetron motion. 

The adjoint matrix A* has eigenvectors orthogonal to the ip's. They 
are given by <£>+,± = {^i{uj 2 /uj + )xT-> Xt) with eigenvalues ±iu + and <f-,± = 
(-(uj 2 z /uj-)x t , Xt) with eigenvalue =fiu-. 

Since (3 is small, the eigenfrequencies of A + /3V can be computed in 
first order perturbation. The cyclotron modes attains a negative real part 
corresponding to the friction coefficient 

7 + = ^4 ("5) 

D7TC 77&o C<J + — C<J_ 

and the magnetron mode attains a positive real part corresponding to the 
antifriction coefficient 

7- = • 7 ' 36 

OTTC 771q LU- — LO + 

As the electron radiates it lowers its potential energy by increasing the mag- 
netron radius. 

Experimentally B = 6 x 10 4 Gauss and the voltage drop across the trap 
is 10V. This corresponds to u z = 4 x 10 s Hz, uj + = 1.1 x 10 12 Hz, cj_ = 
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7.4 x 10 4 Hz. The conditions ( [Tgg ), ([T2§|) are easily satisfied. For the life- 
times one obtains (l/7 2 ) = 5 x 10 s sec, (I/7+) = 8 x 10~ 2 sec, and — (I/7-) = 
2 x 10 23 sec. Thus the magnetron motion is stable, as observed by keeping 
a single electron trapped over weeks. The cyclotron motion decays within 
fractions of a second. The axial motion is in fact damped by coupling to the 
external circuit and decays also within a second. 

The variation with the magnetic field is more clearly discussed in terms 
of the dimensionless ratio {uo c /uo z ) = A. Then 

co± = cj z ^(\± v / A 2 ~^2) , 

7± = ±(3u 2 z {\± v / A 2_ ^2) 3 /8v / A 2_ ^2. (7.37) 

For large A, uj + = A, w_ = A -1 , whereas 7+ = A 2 ,7_ = A~ 4 . As A -> y/2, 
we have uj + = uj_ = u z /\/2. However the friction coefficients diverge as 
(A — \/2)~ l l 2 . Let us call B c the critical field at which the mechanical motion 
becomes unstable. For B > B c , one has still periodic motion with frequency 
/y/2, but the coming instability is disclosed through the vanishing lifetime. 
In the mentioned experiment A = 2.7 x 10 3 and for fixed uj z the critical field 
strength would be B c = 30 Gauss. 

Notes and References 

ad [/].• The name Lorentz-Dirac is standard but historically inaccurate. Some 
authors, e.g. Rohrlich (1997), therefore propose Abraham-Lorentz-Dirac in- 
stead. The radiation reaction term was first derived by Abraham (1905), 
compare with Sections [| and |6] von Laue (1909) realized its covariant form. 
In the Pauli Handbuch article on relativity the equation is stated as in ( [7.1|) . 
The contribution of Dirac is explained in Section |2.2j . 

ad \7. | : The literature on the critical manifold of the Lorentz-Dirac equation 
is listed in ad \6.3j . The Landau-Lifshitz equation appears in all editions of 
their Course in Theoretical Physics. They provide no hint on the geometrical 
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picture of the solution flow nor on the errors involved in their approximation. 
It is rather surprising that the contribution of Landau and Lifshitz is ignored 
in essentially all discussions of radiation reaction, one notable exception be- 
ing Teitelbom et al (1980). For that reason the Landau-Lifshitz equation was 
rederived independently in Spohn (1998). There have been other attempts to 
replace the Lorentz-Dirac equation by a second order equation, Mo, Papas 
(1971), Bonnor (1974), Parrot (1987), Ford, O'Connell (1991,1993). Based 
on Ford, O'Connell (1991), Jackson (1999) uses the substitution in the case 
of a radiation damped harmonic oscillator and discusses applications. In the 
general case only Landau and Lifshitz obtain the correct center manifold 
equation. 



ad \7. Uniform acceleration is discussed in Rohrlich (1990). Constant mag- 
netic field is important for synchroton sources. Since the electron is kept 
on its circular orbit, Larmor's formula is precise enough. Landau and Lif- 
shitz (1959) give a brief discussion. The power law for the ultra-relativistic 
case is noted in Spohn (1999a). Shen (1972a,1978) discusses at which field 
strengths quantum corrections will become important. His results are only 
partially reliable, since he does not start from the Landau-Lifshitz equation. 
The Penning trap is reviewed by Brown and Gabrielse (1986), which includes 
a discussion of the classical orbits and their life-times. They state the re- 



sults ( |7.35|) , ( |7.36| ) as based on a quantum resonance computation. Since 



the result is classical, it must follow from the Landau-Lifshitz equation, 
Spohn (1999a). In the classical framework, more general trap potentials can 
be handled through numerical integration routines for ordinary differential 
equations. 
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8 Spinning Charges 



If an electron is modelled as a classical lump of highly concentrated charge, 
then merely by the interaction with its own radiation field the charge dis- 
tribution will start to rotate. A proper mechanical description must include 
then the angular velocity of the internal rotation and an equation for the 
torque. The argument seems to leave little choice and in this chapter we will 
progress a few steps in the direction of including the classical spin, which 
leads to unexplored and interesting territory. 

If we take the quantum mechanical description as starting point, however, 
as we should do, then the situation is more ambivalent. To be a little bit more 
specific we consider a quantum particle subject to slowly varying external 
forces, which is the standard semiclassical limit. The center of the wave 
packet evolves then according to an effective classical evolution equation of 
the form (|6.34 ). Of course, the energy-momentum relation E S (P) has to 
computed now from the quantum hamiltonian. Only if the model is fully 
relativistic, we can be sure a priori that E S (P) = (P 2 + m 2 ) 1 / 2 . Let us 
assume that in addition the particle carries a spin \. The corresponding 
spinor, f f 6 C 2 , is governed by 

»4^w ( s - 38 ) 

at 

with Hs{t) the time-dependent spin hamiltonian. In a relativistic theory we 
have 

= £ r ■ Kf " 1 + ^ B -<-M> " (§ - !) 2 < 8 - 39 > 

(«(t) ■ B„(r(t)))«(t) - (f - T X-)c- 1 «(t) x £„(r(t))] . 

Here cr = (o^, a 2 , 03) are the Pauli spin matrices with of = 1 and commu- 
tation relations [oi, Oj] = 2ieijk&k, hi, k = 1,2,3. g is the gyromagnetic ratio 
of the quantum particle, g = 2 for an electron. The spin passively adjusts 
itself to the fields along the semiclassical orbit t 1— > (r(t),u(t)) traced out 
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by the particle. To leading order there is no back reaction onto the transla- 
tional degrees of freedom. As we will see this is not the case for a classical 
spin and both degrees of freedom are coupled, a property which is shared by 
relativistically covariant Lagrangians for a particle with spin. We conclude 
that, in contrast to the translational degrees of freedom, a model including 
the classical spin serves only within limits as a phenomenological description 
for a quantum spin. 

There is another, physically more basic objection. Hs(t) is linear in cr. 
Thus defining the average spin s t =< fy, cri/j t > we see that in flOj) h 
drops out and s t satisfies a classical spin equation. In fact, ( |3.39| ) becomes 
the BMT equation for s t , an equation originally obtained on purely classical 
grounds. However, such an approximation by a classical angular momentum 
is valid only in the large spin number limit. In the semiclassical the spin 
degree if freedom remains fully quantum. Of course, in the standard polar- 
ization experiments, as for example the high precision measurements of the 
gyromagnetic ratio, interference is not probed and the classical picture serves 
well. 

8.1 Abraham model with spin 

Abraham models the charge as a nonrelativistic rigid body. Clearly, a com- 
plete mechanical description must specify both its center of mass and its 
angular velocity, which we denote by u)(t) G M 3 . The spinning charge gener- 
ates the current 

j(x, t) = (v(t) + u(t) X(x- q(t)))p(x - q(t)) . (8.40) 

Therefore the source term in Maxwell equations is modified as 

d t B(x,t) = -V x E(x, t), 

d t E{x,t) = VxB(x,t)-(v(t)+u>(t)x(x-q(t)))p(x-q(t)), 
V-E(x,t) = p(x-q{t), V-B(x,t) = 0, (8.41) 

which satisfies charge conservation, since p is radial. 
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The mass distribution of the rigid body is assumed to have the same form 
factor as the charge distribution, p m (x) = m\ y f{x 2 ). The bare moment of 
inertia is then 

lb = ~ J d 3 xp m {x)x 2 (8.42) 

with corresponding angular momentum S = I^u}. The electric dipole mo- 
ment of the charge distribution p vanishes by symmetry. The magnetic dipole 
moment of the current (|8.40|) is given by 



p, — ^ J d 3 xx x (y + (u) x x))p(x) 



= l^I hU j = g h -^S (8.43) 
2 nib 2rrib 

with bare gyromagnetic ratio — 1. As in the case of the bare mass, lb and 
(7b will be renormalized through the self-interaction. 

The Lorentz force equation comes now in two parts, one for the linear 
and one for the angular momentum. To be consistent we stick to the nonrel- 
ativistic form and have 

— m h v(t) = / d 3 xp(x - q(t))[E(x,t) + (v(t) + u(t) x (x - q(t))) x B(x,t)] , 

^bW(i) = J d 3 xp(x - q(t))(x - q(t)) x [E(x,t) (8.44) 

+ (v(t)+u(t) x (x-q(t))) x B(x,t)] . (8.45) 

If in addition there are external forces acting on the charge, then E and 
B in (|8.44j) , flS.45|) would have to be replaced by E + E ex and B + B cx , 
respectively. 

The Abraham model of Section |L4|is obtained by formally setting u>(t) = 
0. Note that this is not consistent with the Lorentz torque equation fl8.45|) , 
since u) ^ 0, in general, even for uj = 0. 

The Abraham model with spin has the conserved energy 

£ = ~ m h v 2 + l - I h u 2 + ~ J d 3 x(E 2 + B 2 ) (8.46) 
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and the conserved linear momentum 

V = m h v + J d 3 xE x B . (8.47) 

In addition the total angular momentum 

J = qx m h v + I h uj + J d 3 xx x (E x B) (8.48) 

is conserved. Of course, also the Abraham model without spin is invariant 
under rotations and must therefore have a correspondingly conserved quan- 
tity. Only it does not have the standard form of the total angular momentum, 
which from a somewhat different perspective indicates that inner rotations 
must be included. 

As good tradition already, we assume that the external forces are slowly 
varying and want to derive in this adiabatic limit an effective equation of 
motion for the particle including its spin. As a first step of this program we 
have to determine the charge solitons. We set 

q(t) = vt , w(t) = uj , (8.49) 
E(x,t) = E(x-vt) , B(x,t) = B(x-vt) (8.50) 

and have to determine the solutions of 

-v ■ VB = -V x E , 

-v ■ VE = V x B - p(v + uj x x) , (8.51) 
VE = p , V- B = 0, 

0= / d 3 xp(x)[E(x) + (v + u x x) x B(x)), (8.52) 



= J d 3 xp(x)x x [E(x) + (v + uj x x) x B{x)\ . (8.53) 



We will solve (|8.51| ) to ( |8.53| ) in Fourier space. 



The inhomogeneous Maxwell equations (|8.51|) are solved by 

E = E 1 + E 2 , B = B 1 + B 2 (8.54) 
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with 



^i(fc) = -i[k 2 - (k ■ vf]-\k - (k ■ v)v)p(k) , (8.55) 
E 2 (k) = -[k 2 - (k ■ v) 2 }- 1 ^ x k)(v • V fc )p(fc) , (8.56) 



Bi(fc) = i[k 2 - (k ■ v) 2 ]-\k x v)p(k) , (8.57) 
B 2 {k) = -[k 2 - (Jfe • u) 2 ] _1 (fe x (w x V fc ))p(fc) . (8.58) 

Note that Ei, B\ are odd and E 2 , B 2 are even in k. 

Using that odd terms vanish, in the Lorentz force equation we have 

- J d 3 kp*[k 2 - (v ■ fc) 2 ] _1 (w x k)(v • V fc )p 

- J d 3 kp*[k 2 - (v ■ k) 2 ]- x v x (Jfe x (u; x V fc ))p 
+ J d 3 k((u; x V fc )p*)[fc 2 - (v ■ k) 2 ]- 1 x (Jfe x v)p 

= _y ^[^-(vfc) 2 ]- 1 !^- 1 ^ 

((w x fe)(« • k) + v x (Jfe x (u> x Jfe)) — ((a? x Jfe) • v)k) = 

(8.59) 

for every v and u>, where we took into account that p is radial. 

The Lorentz torque equation requires more work. Using again that odd 
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terms vanish, we have 




(k x (fe — (fe • — x (u x (k x «))) 




(8.60) 



For the divergence of B 2 we find 



V fc • B 2 = 2[k 2 - (k ■ v) 2 }- 2 k 2 (uj ■ V fe - (v ■ u)(v ■ V fe ))p (8.61) 



d 3 k\V k p\ 2 2[k 2 - (k ■ v) 2 ]- 2 {k x w)(w • fe - (v ■ u)(v ■ k)) = . (8.62) 



Taking into account that p is radial, the torque vanishes only if either u> || v 
or oj±.v. If v = 0, the torque vanishes always. For u) oblique to v Equations 
( |8.51| ) to ( |8.53| ) have no solution. 

Physically the charge distribution is rigid, but the electromagnetic fields 
are Lorentz contracted along v. This mismatch yields a nonvanishing torque 
unless oj || v , resp. iV-Lv. Clearly, this mismatch is an artifact of the semi- 
relativistic Abraham model. For a relativistic extended charge distribution 
there is a charged soliton for every v and uj. 

In the adiabatic limit, there must be two disjoint effective equations of 
motion. If lo || v initially, then it will remain so approximately and there is 
a closed equation for q, v and uj ■ v. Similarly, if <jj1.v initially, this property 
is almost preserved in time, v and oj — (v ■ uj)oj vary slowly on the same 
time scale. We expect that through the emission of radiation an uj initially 



and therefore zero Lorentz torque results in the equation 
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oblique to v will rapidly relax to either being parallel or orthogonal to v, 
depending on the initial conditions. This admittedly rather sketchy picture 
raises some interesting dynamical questions, in particular how precisely the 
spinning particle succeeds in slowly turning its axis of rotation. We have not 
persued this issue, since it leads away from the quantum spin. 



8.2 Relativistic dynamics of charged particle with spin 

To be supplied. 



Notes and References 

ad\8[. BMT is an acronym for Bargmann, Michel, Telegdi (1955). The BMT 
equation is explained in Jackson (1999). Bailey, Picasso (1970) is an infor- 
mative article on how the BMT equation is used in the analysis of the high 
precision measurements of the electron and muon ^-factor. The BMT equa- 
tion with g = 2 is the semiclassical limit of the Dirac equation, Rubinow, 
Keller (1963), Bolte, Keppeler (1999), Spohn (1999c). 



ad \8.1\ Kiessling (1999) observes that the usual form of the total angular 
momentum is conserved only if the inner rotation of the charged particles is 
included. The charge soliton solutions with spin are determined in Spohn 
(1999b). 



ad \8.3[ . Just as for translational degrees of freedom, one way to guess the cor- 
rect effective spin dynamics is to impose Lorentz invariance. In addition, one 
could require that the equations of motion come from a Lagrangian action. 
In full generality, including an electric dipole moment, this program was car- 
ried out by Bhabha, Corben (1941). Alternative approaches are compared 
in Corben (1961), Nyborg (1962). A concise and useful summary is given by 
Barut (1964) who discusses also how the BMT equation fits into the general 
scheme. In the relativistic extended charge model of Nodvik (1964) spin is 
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included and the effective equation of motion is derived formally. Recent 
work on this model is Appel, Kiessling (1999). 
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9 Many Charges 



There is little effort in extending the Abraham model to several particles. 
We label their positions and velocities as qj(t), Vj(t),j — 1, . . . ,N. The j-th 
particle has the bare mass m^j and the charge ej, where for simplicity all 
particles have the same form factor /, i.e. Pj(x) = ejf(x 2 ). The motion 
of each particle is governed by the Lorentz force equation, as before, and 
the current in the Maxwell equations becomes now the sum over the single 
particle currents. Therefore the equations of motion read 

d t B(x,t) = -VxE(x,t), 

N 

d t E(x,t) = VxB(x,t)-J2pA x -q 3 ( t ))v J (t), (9.1) 

3=1 

N 

V-E(x,t) = ^ Pj (x-q j (t)),V-B(x,t) = 0, 

3=1 

d f 

-(m him (t)) = d 3 xp(x-q l (t))(E(x,t) + v l (t)xB(x,t)), (9.2) 

i = l,...,N, with 7i = (1 - v 2 )' 1 ' 2 . 

There are no external forces. Thus the force acting on a given particle is 
due to the other particles as mediated through the Maxwell field. In addition, 
there is the self-force which we have discussed already at length. Physically 
we trust our model only if particles are very far apart on the scale set by R p . 
If two particles are at a distance of a few R P) then there are strong forces 
which depend on the details of the phenomenological and unknown charge 
distribution. Thus we assume that initially 

\<£-<$\=0(e- 1 R p ), i^j. (9.3) 

We emphasize that the scale parameter e enters only through the initial 
conditions. e~ l is the typical distance of particles measured in units of R p . 
For the initial fields it is natural to again impose the condition of no slip. 
Then they are a linear superposition of charge soliton fields corresponding to 
ql vl i = l,...,N. 
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We expect ( |9.3| ) to remain valid at least over a certain macroscopic time 
span and we want to understand whether in this limiting regime there is 
a closed dynamics for the particles by themselves. Since the particles are 
far apart, it takes a time of order e~ 1 R p /c = e~ 1 t p for light to travel inbe- 
tween and the force on a given particle depends on the other particles at 
a macroscopically retarded time. This means that the effective equations 
of motion are closed, but nonlocal with a structure to be explained in the 
following section. In many circumstances the velocities can be regarded as 
small, |i>j|/c <C 1, and retardation effects are neglible. To lowest order this 
yields then the static Coulomb interaction. Somewhat unexpected even the 
first order correction has still the form of an effective Lagrangian. 



9.1 Retarded interaction 

We insert the solution of the inhomogeneous Maxwell equation ( |9.1| ) into the 
Lorentz force equation ( p.2|) . The forces are additive and the force on particle 
i naturally splits into self-force (j = i) and a mutual force (j 7^ i). For the 
self-force we use the Taylor expansion of Section ^. Thereby the mass is 
renormalized and the next order is the radiation reaction. For the mutual 
force we recall that in Section |6.4| we showed already that to leading order 



the field generated by charge j is the Lienard-Wiechert field. Thus, ignoring 
radiation reaction, we obtain as retarded equations of motion 

N 

mi{vi)vi = ^ e j (E retj (q i ,t) +«,x B reti (^,t)) , t > 0. (9.4) 

j=i 

Here rrii is the effective mass of particle i as defined in ( |6.1| ). E Tet j(x,t) 
equals ( |1.21| ) with q replaced by and t ret replaced by t Tet j which is defined 
by 

tretj = t-\x- ^-(tretj)! • (9.5) 

For x = q i the retarded time is 0(e~H p ). Similarly B ret j(x,t) equals ( |1.22| ) 
with q replaced by q^ and t ret replaced by t ret j. Note that the equations of 
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motion in ( |2.14| ) have the same structure. 

To solve (|9.4j) one needs the trajectories for the whole past. Our assump- 
tion of no initial slip is equivalent to 

q^ + qV + tv* , i = l,...,N , t<0 (9.6) 

which must be added to ( ^4|) . 

Using ( |9.4| ) we can estimate the size of the various contributions. The 
far field contributions to E TCt j and B TCt j are of 0(e 2 ) and the near field 
contributions are 0(e)vj. Thus i)i = 0(e 2 ) and the first order correction 
from the near field is 0(e 3 ). In the next order we see the radiation reaction 
which is proportional to v 2 and thus 0(e A ). In ( |9.4|) we would have to add 
the Lorentz-Dirac term of ( |6.1| ) for particle i. As a consequence ( |9.4| ) will 
pick up runaway solutions and we have to restrict to the critical manifold. 
Thereby the friction force becomes retarded. This is not surprising. The 
friction comes from acceleration which is due to the retarded motion of all 
other charges. 

The issues raised here remain largely unexplored, at present. One would 
like to know how well ( ^.4| ) approximates the true dynamics and over what 
time scale. 

9.2 Limit of small velocities 

We impose that the initial velocities are small. The natural scale turns out 
to be 

1^1 = 0(^0) (9.7) 

and, of course, we have to show that this order is maintained, at least for a 
certain time span. To preserve the relation q = v we have to adjust the time 
scale as e~ 1//2 relative to the macroscopic time scale. The accumulated force 
is then of order yfe which just balances the velocity in (|9.7f) . Therefore we 
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arrive at the following scale transformation 

t = e- 3/2 t', q 3 = e-'q'j , Vj = y/i^ , (9.8) 
x = e~ 1 x', E = e 3 ' 2 E' , B = e 3 ^ 2 B' , 

where the primed quantities are considered to be of 0(1). The field ampli- 
tudes are scaled by e 3 ! 2 so to preserve the field energy. 

There is little risk of confusion in omitting the primes. We set 

q){t) = eq 3 (e^H) , v]{t) = e^v^H) . (9.9) 

Then the rescaled Maxwell-Lorentz equations are 

y/ed t B{x,t) = -VxE(x,t), 

N 

y/idtE(x,t) = WxB(x,t)-J2y/iv £ j (t)y/ip jte (x-q%t)),(9.10) 

3=1 

N 

V • E(x, t) = VzPlAx ~ q £ j(t)) , V • B(x, t)=0, 

3=1 



sf t (m hi (l ~ evm 2 r 1/2 vm) (9.H) 
= J d 3 xy/Ip £ , i (x-qt{t))(E(x,t) + y/Ivt(t)xB(x,t)). 

On the new scale the velocity of light tends to infinity as c/ yfe and the 
charge distribution has total charge y/e, finite electrostatic energy m e , and 
shrinks to a 5-function as y/ep £: j(x) = y/e e~ 3 pj(ex). Recall that the scale 
parameter e is just a convenient way to order the magnitudes of the various 
contributions. 

Before entering into more specific computations, it is useful to first sort 
out what should be expected. We follow our practice from before and denote 
positions and velocities of the comparison dynamics by Vj, Uj,j — 1, . . . , N, 
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i.e. q £ j(t) = rj(t),Vj(t) = Uj(t). Since the velocities are small, the kinetic 
energy takes its nonrelativistic limit 

T oK) = ~ (m hj + i m cj )u 2 , (9.12) 



up to a constant, compare with ( |3.19|) . Note that the mass of the particle 



is renormalized through the interaction with the field. For small velocities 
magnetic fields are small and retardation effects can be neglected. Thus the 
potential energy of the effective dynamics should be purely Coulombic and 
be given by 

1 N 

U Q {r 1 ,... : r N ) = ~Y,TT^—\- ^ 

2 47r|r, ; — r -l 



So what is the next order? For the kinetic energy we merely expand in 
( |3.19p with the result 



Ti(uj) = e(-m hj + — m^)u$ • (9.14) 

The next order correction to the Coulomb forces requires more explicit con- 
siderations, which will be explained in the following section. There are cor- 
rections due to retardation and to the magnetic field, which combine into a 
velocity dependent potential as 

1 N ee- 

Ui(r u ui, . . . ,r N , u N ) = -e- ^ 4 i r ' 3 _ \ { v i- v i + i v i ' ' v j)) 



3 I 



(9.15) 

with ?ij = (ri-r^/ln-rjl. 

In principle we could continue the expansion. It is of interest to see at 
what scale radiation effects will be important. They are proportional to v 2 . 
On the microscopic scale v = e 2 , as argued before. In rescaled velocities and 
accumulated over the time span e~ 3 ^ 2 this results in a loss of energy of the 
order e 3 ^ 2 . Thus the next order correction to the comparison dynamics is 
dissipative and of order e 3//2 . 
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We recall that \vj\/c = 0(y/e). Thus we may set e = 1 at the expense of 
reintroducing the velocity of light, c. Then up to an error of order (\vj\/c) 3 ^ 2 
the effective dynamics of the N charges is conservative and is governed by 
the Lagrangian 

N 4 11 2 

^Darwin = ((>bj + -m ej )-U 2 + (-77l b j + —m ej )c~ 2 u1j) 

3=1 

1 N 1 

-^Darwin is known as the Darwin Lagrangian and widely used in plasma physics. 



9.3 The Darwin Lagrangian 



As can be seen from (|9.1]j , (|9.2| ) the forces are additive. Thus it suffices 
to consider two particles only. As initial conditions we choose the linear 
superposition of the two charge solitons corresponding to the initial data 
q®, v® , i = 1,2. We solve Maxwell equations and insert in the Lorentz force. 
As already explained, in the self-interaction the contribution from the initial 
fields vanishes for t > t p . In the mutual interaction the initial fields take a 
time of order \fe to reach the other particle and their contribution vanishes 
for t > \fe\q\ — q 2 \- Thus for larger times we are allowed to insert in (|9.2|) 
the retarded fields only, which yields 

e^("lbl7l«l(*)) = ^rct,llW + ^ret,12 00, (9.16) 

ej t {m h2l2 v £ 2 {t)) = F ret>21 (t) + F rct , 22 (t) , (9.17) 

where 

t 

F ret , tJ (t) = J ds J d 3 fcp;(efc)p,(efc)e ifc - ( ^ (t) -^ (s)) 

- e 1/2 (|fc| -1 sm(\k\(t - a)/y/e))ik - e(cos(|fc|(i - s)/Vi))^( s ) 
+e 3/2 (|fc|- 1 sin(|fc|(t- s)/yfe)v\{t) x (ik X ^(s))) , (9.18) 
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i, 3 = 1,2- 

For the self-interaction we set ek = k' , e~ 3 l 2 t = t' . Then 

oo 

F^Mt) = e 3/2 y'dry'rf 3 A;|p 1 (fc)|V fc -^W-' ? ^ t - e3/2r »/ e 



( - y/edkl" 1 sin |fe|r)ife - £:(cos \k\r)vl(t - e 3/2 r) 

+e 3/2 {\k\- 1 sin |fc|r)«?(t) x (ik x v\{t - £ 3/2 r))) .(9.19) 
We Taylor expand as 



e- l {q\{t)-q\{t-eV 2 T)) = e^rv - l -e 2 r 2 v 

V\{t-8 3 ' 2 T) 



v - e 3/2 rv 



(9.20) 



Then, up to errors of order e 5 ^ 2 , 

oo 

F Iet , u {t) = JdrJd 3 k\p 1 {k)\ 2 {e[-{\k\- 1 sm\k\r)^T 2 (k-v)k 
o 

+ (cos \k\r)Tv] +e 2 [( - (|fc| _1 sin \k\r)^r 2 (k ■ v)k 

+ (cos \k\r)rv) ( - ]^ 2 (k ■ v) 2 ) 

-(cos \k\r)^T 3 (k -v)(k- v)v + (|fc| _1 sin \k\r) 

(r 2 (k ■ v)v x (k x v) + ^r 2 (fc -v)(vx(kx «)))] }(,9.21) 



which, upon integration, agrees with (|9.12| ), (|9.14j ). 

For the mutual interaction we leave the /c-integration and set e x l 2 t = t' . 
Then 



ret. 12 



oo 

(t) = ^ J dr J d 3 kpl(ek)p 2 (ek)e tk < q ^ t) - q ^ t - v ~ £T)) 



- £ 1/2 (|fc| _1 sin \k\r)ik - e(cos \k\r)v £ 2 {t - yfer) 
+e 3/2 {\k\- 1 sin\k\T)v £ 1 (t) x (ik x v e 2 (t - y/er))) . (9.22) 
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We Taylor expand as 



<?!(*) - Ql(t ~ y/er) = r + y/erv 2 - - st 2 v 2 , 

«?(*) = vi , v £ 2 {t - y/er) = v 2 - y/erv 2 (9.23) 



with r = q\(t) — q £ 2 {t). Then, up to errors of order e 5//2 , 

oo 

J dr J rf 3 A;p*(efc)p 2 (efc)e ifc - r | - e^ 1 sm\k\r)ik 



ret, 12 



+e 2 [(Ik]- 1 sin \k\r) {-\^ 2 {k- v 2 )k + l - r 2 (k ■ v 2 ) 2 ik 
+v 1 x (ik x v 2 )) + (cos \k\r)(rv 2 - ir(k ■ v 2 )v 2 )] j 

= (e ie2 /47r) ( - eV^rl" 1 + e 2 [(- V r (i> 2 • V r ) - ^V r (v 2 ■ V r ) 2 ) |r| 
-(« 2 - v 2 (u 2 • V.))^!- 1 + ( Vl x (V r x v 2 ))|r|- 1 ]) . (9.24) 
We define the potential part of the Darwin-Lagrangian for two particles 



as 



1 e v-i ■ Vo £ (v-i • r)(v 2 • r)\ 

L P = (e 1 e 2 4n)(-—+-^—^+- yi >y2 > ). 9.25 
V \r\ 2 \r\ 2 |r| d / 

Then 

-OO = ~e(j t (V V1 L P ) - V r L P ) + 0(£ 5/2 ) . (9.26) 



ret, 12 1 



Inserting ( |9.21| ) and ( |9.26| ) into the Lorentz force equation ( |9.16| ), ( ^.17| ), we 
conclude that, upon neglecting contributions of order e 5 / 2 , the dynamics of 
the charges is governed by the Darwin Lagrangian (|9.16|) . To control our 
Taylor expansion one has to resort to the contraction argument of Section 



5T3| . It becomes now considerably more involved. We note that in the effective 
equations of motion there is no mechanism which would preclude head on 
collisions. Thus the dynamics governed by L Darwin can hold only until the 
first collision. For the subsequent motion one has to go back to the full 
microscopic evolution. 
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Notes and References 

ad \9.Q and The Darwin Lagrangian is discussed in Jackson (1999). In 
Kunze, Spohn (1999b) the errors relative to the motion governed by the 
Darwin Lagrangian are estimated. 
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